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ACT  ( Continue 

LORA  is  a model  developed  to  predict  the  performance  of  long-range  active  sonar  sys- 
tems. The  supporting  computer  program  uses  3-8  epu  seconds  per  run  on  the  UNiVAC 
1110  computer  and  requires  38,000  words  of  core  storage.  The  performance  prediction 
model  comprises  models  that  estimate  propagation  loss,  reverberation  level,  and  probability 
of  detection.  The  propagation  loss  models  arc  ( I ) empirical  equations  (modified  AMOS)  for 
near-surface  ducts  and  (2)  ray  theory  (Pedersen  and  Gordon)  for  direct-path,  bottom-bounce. 
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SUMMARY 


OBJECTIVE 

Develop  a computer  model  for  predicting  the  performance  of  active  sonar  sys- 
tems for  long  horizontal  ranges,  that  is,  for  multiple  convergence  zones  and  bottom 
bounces. 


COMPUTER  PROGRAM  CAPABILITIES 

The  model  and  supporting  computer  program  have  been  developed  and  embrace 
the  following  capabilities: 

1 . Active  sonar  performance  predictions  are  available  for  up  to  five  convergence 
zones  or  bottom  bounces. 

2.  The  target  and  active-sonar  source  may  be  situated  at  any  ocean  depth. 

3.  Optional  passive  coherent  and  incoherent  propagation  losses  are  available  for 
up  to  30  convergence  zones  or  bottom  bounces. 

4.  Curvilinear  techniques  are  used  with  the  sound-speed  profile  to  reduce  false 
caustics. 

5.  Multi-ping  reverberation  is  available,  as  well  as  doppler  shifting  of  the  signal 
caused  by  target  motion. 

6.  Five  probability-of-detection  models  for  different  assumptions  of  signal 
degradation  and  type  of  detector  are  available. 

7.  The  computer  program  is  able  to  perform  many  successive  runs  without 
failure. 

8.  Computer  costs  are  low,  approximately  $2  for  a complete  set  of  perform- 
ance predictions  using  an  average  commercial  computer. 

9.  Core  storage  required  is  about  38,000  words. 


RECOMMENDATIONS 

The  development  of  an  adequate  subsurface  duct  model  is  needed,  perhaps  using 
normal  mode  theory.  Propagation  loss  in  the  duct  should  be  related  to  duct  width, 
duct  strength,  source  and  target  positions,  and  the  frequency  of  the  propagating  sound. 

Volume  reverberation  calculations  could  be  improved  by  averaging  propagation 
losses  from  transmitter  to  scatterers  over  the  whole  scattering  volume.  Nonspecular 
scattering  could  be  included. 

The  input  package  of  the  computer  program  could  be  modified  to  include  a 
curve  of  recognition  differential  versus  signal-to-noise  ratio  to  render  the  program 
more  sensitive  to  certain  active  sonars. 
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INTRODUCTION 


In  order  to  assess  the  effectiveness  of  long-range  weapon  systems,  an  active 
sonar  model  known  as  LORA  has  been  developed.  LORA  is,  essentially,  an  extension 
of  FAST  NISSM,1  the  utility  version  of  the  Navy  Interim  Surface  Ship  Model 
(NISSM).2  LORA  was  the  next  logical  iteration  in  a series  of  models  developed  to 
predict  active  sonar  performance.  Its  supporting  computer  program  is  a high-speed, 
user-oriented  program,  with  a flexible  input  routine.  Program  execution  is  rapid  and 
inexpensive,  costing  less  than  $2  (4  cpu  seconds  plus  input-output  charges)  for  a full 
set  of  near-surface  (direct-path),  multi-bottom-bounce,  and  multi-convergence-zone 
performance  predictions. 

The  LORA  model  affords  analysts  the  opportunity  to  evaluate  the  long-range 
detection  capabilities  of  operational  and  future  active  sonar  systems.  Specifically, 
LORA  can  be  used  to 

• Compare  the  performance  of  candidate  sonar  systems  under  identical 
environmental  conditions. 

• Evalute  the  performance  of  a given  sonar  system  under  various  environmental 
conditions. 

• Analyze  the  sensitivity  of  performance  to  changes  in  selected  environmental 
and  system  parameters. 

• Estimate  the  improvements  in  performance  resulting  from  proposed  changes 
in  system  components. 

• Establish  a set  of  performance  reference  data,  which  would  make  possible  the 
comparison  of  performance  measurements  made  in  different  environments. 


The  body  of  this  report  contains  the  total  model  description.  First,  the  active 
sonar  equation  is  presented  to  illustrate  the  procedure  used  to  determine  signal 
excess.  Next,  models  are  given  for  the  various  kinds  of  losses  which  degrade  transmit 
signal  and  target  echo.  Then  noise  and  reverberation  models  are  given.  F'inally,  several 
models  for  probability-of-detection  are  presented  for  various  assumptions  of  signal 
distortion  and  detector  characteristics. 

The  appendixes  contain  information  about  the  LORA  computer  program. 
Appendix  A discusses  the  logical  structure  of  the  program;  Appendix  B defines  the 
input  parameters;  Appendix  ('  describes  the  input  format;  and  Appendix  D contains 
sample  runs. 
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MODEL  DESCRIPTION 


The  model  considered  here  is  capable  of  predicting  the  performance  ot  mono- 
static, active  sonar  systems,  hull-mounted  or  towed.  Some  of  the  salient  features  ot 
the  model  include: 

• No  limitation  to  source  and  target  depths  within  the  ocean. 

• Automatic  curve  fitting  for  the  sound-speed  profile  to  reduce  the  effect  of 
false  caustics  caused  by  the  digitizing  of  the  sound-speed  data. 

• Multiple  convergence-zone  predictions  (maximum  of  five). 

• Multiple  bottom-bounce  predictions  (maximum  ol  five). 

• Multiple-ping  time  histories  tor  reverberation  (maximum  of  nine  pings). 

• Choice  of  five  probability-of-detection  models. 

• Optional  passive  coherent  and  incoherent  propagation  losses  for  a maximum 
of  30  bottom  bounces  or  convergence  zones. 

Model  inputs  are  basic  descriptive  variables  defining  the  ocean  medium,  the 
transmitting  and  receiving  system,  and  the  signal  processor.  Predictions  are  restricted 
to  active  sonar  systems  which  operate  at  frequencies  between  0.5  and  35  kHz.  Two 
restrictions  have  been  imposed  on  environmental  representations  to  facilitate  the 
application  of  high-speed  ray-tracing  techniques.  Eirst.  horizontal  homogeneity  ot  the 
water  column  is  assumed;  that  is,  a single  sound-speed  profile  is  specified  along  any 
range  track  of  interest.  Second,  both  surface  and  bottom  boundaries  are  assumed 
locally  flat  with  a constant  bottom  depth,  although  corrections  are  made  to  account 
for  earth  curvature. 

Model  output  is  probability-of-detection  for  those  ranges  from  the  active  sonar 
which  can  be  irradiated  at  a specific  target  depth.  Intermediate  quantities  such  as 
propagation  loss,  reverberation  level,  signal  level,  and  signal  excess  are  supplementan 
calculations  and  are  useful  outputs  in  their  own  right.  Model  implementation  is 
effected  by  a FORTRAN  computer  program  which  requires  about  38.000  3ft-bit 
words  on  the  UNI  VAC  1110. 

Probability-of-detection  is  used  as  the  measure  of  sonar  system  performance 
and  is  calculated  as  a function  of  horizontal  range  between  source  (own  ship)  and 
target  A useful  measure  of  acceptable  performance  can  be  defined  as  those  range- 
intervals  where  the  probability-of-detection  exceeds  some  predetermined  probability 

level. 

Single-ping  detection  probability  is  a function  of  signal  excess,  which  is  the 
amount  by  which  the  target  echo  level  exceeds  the  masking  background  level. 


The  signal  differential,  LSJX,  often  referred  to  as  the  signal-to-noise  level,  is  the 
difference  within  a designated  frequency  band  between  the  echo  level,  E,  and  the 
masking  background  level,  M\  that  is 

Lsn  = E-M  (dB)  ( 1 ) 

The  minimum  recognizable  value  of  echo  level,  EM  , is  that  portion  of  the  masking 
background,  M,  such  that 

Em  = M + Rd  (dB)  (2) 

where 

Rp  is  the  recognition  differential 

The  recognition  differential  provides  a measure  of  sonar  operator  performance.  It  is  a 
psycho-acoustic  measure  of  a sonar  operator’s  ability  to  recognize  a target  echo  50% 
of  the  time  given  a specified  probability  of  false  alarm.  In  terms  of  EM,  any  echo  level 
can  be  expressed  as 

£ - £m  + St;x  (3) 


where 

St  x is  referred  to  as  signal  excess 

If  a signal  is  transmitted,  then  for  a given  target  strength,  ST.  the  echo  level  is  given  by 


E - L{)  + Sy  - //•> 


(4) 


where 


and 


L{)  is  the  signal  level  1 yard  from  the  source 


/A,  is  the  two-way  propagation  loss  plus  the  degradation  accrued  from  the 
transmit  and  receive  beam  patterns 


Signal  excess.  Stx,  is  obtained  by  combining  Eq.  2.  3,  and  4: 


ShX  = L0  +St-H2  -M-Rd 


(5) 


This  is  the  active  sonar  equation  for  echo  ranging. 


PROPAGATION  LOSS 


f 


Sound  pressure  is  the  root-mean-square  increase  in  pressure  over  hydrostatic 
pressure  due  to  the  propagation  of  harmonic  acoustic  energy  through  the  position  of 
measurement.  Propagation  loss  is  a measure  of  the  reduction  of  far-field  sound  pres- 
sure from  that  of  a reference  sound  pressure  at  a specified  distance  from  the  sound 
source.  The  expression  for  propagation  loss  is 

//  = -10  log  (p2lp^ ) in  dB  (6) 


where 

p is  the  sound  pressure  at  some  range  of  interest 

and 

p0  is  the  reference  sound  pressure  level  at  a reference  range 

This  model  assumes  that  four  physical  processes  in  which  the  propagating  acoustic 
energy  interacts  with  the  ocean  environment  contribute  to  a reduction  (in  some 
refractive  cases  an  increase)  in  sound  pressure  and  therefore  an  increase  in  propaga- 
tion loss.  These  processes  are  refraction-spreading,  absorption,  surface  reflection,  and 
bottom  reflection  of  the  acoustic  pressure  waves.  These  four  processes  are  assumed  to 
contribute  independently  to  the  total  propagation  loss: 

Hr  = Hr  + H,  + !fs  + Hr  in  dB  (7) 


where 

is  refraction-spreading  loss 
//,  is  absorption  loss 
//s.  is  surface-reflection  loss 

and 

Hfl  is  bottom-reflection  loss 

1 ach  of  these  four  processes  has  the  following  characteristics  in  common.  Let  />,  be 
the  sound  pressure  at  the  start  of  the  process  and  p2  be  the  sound  pressure  at  the 
enu  of  the  process,  then  the  ratio  p \ lp2  is  independent  of  pl . This  characteristic 
results  directly  in  the  associativity  of  the  terms  on  the  right  side  of  Lq.  7.  since  pi 
can  be  identified  with  p0  and  p2  with  p in  Lq.  6.  Thus,  the  components  of  propaga- 
tion loss  can  be  calculated  independently  and  assembled  in  any  order. 

Several  techniques  for  determining  propagation  loss  are  used  in  this  model, 
depending  on  the  sound  speed  and  the  source  and  target  depths.  The  empirical 
Acoustic.  Meteorological,  and  Oceanographic  Survey  (AMOS)  equations3  are  used 
when  surface-duct  calculations  are  warranted.  Simple  subsurface-duet  equations  are 
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used  if  both  source  and  target  are  in  the  duct.  Otherwise,  direct-path  calculations  are 
performed  using  ray-tracing  techniques  based  on  a sound-speed  profile  fitted  with 
continuous-gradient  curvilinear  segments.  Bottom-bounce,  convergence-zone, 
surface-bounce,  and  reliable-acoustic-path  calculations  are  performed  by  using  curvi- 
linear ray-tracing  techniques.  For  active  sonar  predictions,  only  the  path  yielding  the 
minimum  propagation  loss  at  each  range  is  used  for  determining  propagation  loss 
versus  range.  Average  losses  for  boundary  reflections  are  obtained  as  a function  of 
frequency  and  angle  of  incidence.  The  number  of  surface  and  bottom  reflections  are 
counted,  and  boundary  loss  is  calculated  and  added  into  the  total  loss  fora  path.  The 
result  is  the  total  propagation  loss  caused  by  the  environment,  which  is  represented 
by  Fq.  7.  Beam  pattern  responses  (transmit  and  receive)  also  degrade  the  intensity 
for  a ray  path.  Beam  pattern  responses  in  dB.  //<£.  are  added  to  Hr  of  Eq.  7 to  give 
the  propagation  loss  used  in  comparing  paths  competing  to  give  minimum  propaga- 
tion loss  to  the  target. 

Table  1 gives  a summary  of  the  components  needed  to  compute  propagation 
loss  for  several  modes  .,f  acoustic  propagation. 


Table  1.  Applicability  of  propagation-loss  terms  to  various  propagation  modes. 


Propagation  Mode  Propagation  Loss  Equation 


Surface  Duct  (SD) 

Subsurface  Duct  (SSD) 

Direct  Path  (DP) 

Bottom  Bounce  (BB) 
Convergence  Zone  (CZ) 
Surface  Bounce  (SB) 

Reliable  Acoustic  Path  (RAP) 
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Passive  one-way  propagation-loss  outputs  are  also  available.  The  options  are: 

1 . Coherent  summation  over  all  the  paths  contributing  to  acoustic  pressure  at 
the  target. 

2.  Incoherent  summation  over  all  the  paths  contributing  to  acoustic  pressure  at 
the  target. 

3.  Fading  coherent  summation  of  all  the  paths  contributing  to  acoustic  pressure 
at  the  target.  (This  option  is  simply  a linear  decrease  of  coherence  with  range  until 
incoherence  is  obtained  at  50  kyd.  It  is  admitted  that  this  is  not  a realistic  fading 
coherence  model.) 

4.  Passive  simulation  in  which  the  receive  beam  pattern  is  included  in  the  pres- 
sure calculation  for  each  ray  path  before  summation. 

5.  Communications  simulation  in  which  the  transmit  beam  pattern  response 
(evaluated  for  the  ray  angle  at  the  transmitter)  and  the  receive  beam  pattern 


response  (evaluated  for  the  ray  angle  at  the  target)  are  included  in  the  pressure  calcu- 
lation for  each  ray  path  before  summation.  This  model  is  not  complete  in  that  time 
histories  of  the  signals,  with  appropriate  summation  over  the  phase-distorted,  time- 
displaced  signal  waveforms,  are  not  calculated. 

These  models  have  nothing  to  do  with  active  sonar  modeling  but  are  useful 
byproducts  of  the  calculations  needed  for  active  sonar  performance  predictions.  It 
must  be  noted  that  one-way  passive  propagation  losses  are  not  in  general  one-half  of 
the  two-way  active  propagation  losses.  This  is  because  summation  over  several  paths 
connecting  source  and  target  can  produce,  say,  a 10  dB  increase  over  the  single-path, 
minimum  propagation  loss  used  for  active  sonar  performance  prediction.  When  propa- 
gation loss  is  calculated  for  active  sonars,  it  is  assumed  that  the  pulse  length  is  short 
enough  in  time  to  be  resolved  in  every  multipath  occurrence.  The  assumption  for 
passive  sonars  is  that  the  signal  is  CW  and  constant  over  the  maximum  time  interval 
between  any  two  paths  in  the  summation. 


Absorption  Loss 

The  expression  for  absorption  loss  is 
Ha  = otS 

where 

S is  the  curvilinear  path  length  (kyd) 

and 

a is  the  absorption  coefficient  (dB/'kyd) 

Values  of  a are  obtained  from  an  equation  developed  by  Hall  and  Watson4  that  com- 
bines the  low-frequency  predictions  of  Thorp5  with  the  high-frequency  predictions 
of  Schulkin  and  Marsh.6  The  resulting  expression  is 

1.776 /KS  , 1 - 6.54  X \Q-^p  fo.65053/Tf~  0.026847/ 

( X = + — - I + 

32768  +/3  1 + 32768/3  \ /2 +/f  fT 

where 

/is  frequency  in  kHz 
P is  pressure  in  atmospheres 

and 


6 


fT  is  the  relaxation  frequency  in  kHz 


' 


i 

[ 


Pressure  is  a function  of  depth  in  feet,  z,  so  that 
P = 1.04  + 0.03 1244z  atm 
Tite  relaxation  frequency  is  determined  by 

fT  = 21.9  X 10<3^+‘02)/(5r+2297)kHz 


where 

T is  temperature  in  Fahrenheit 

Figure  1 shows  a as  a function  of  frequency  for  several  average  ocean  tempera- 
tures computed  by  Eq.  8.  A curve  showing  Thorp’s  deep  ocean  prediction  is  presented 
for  comparison. 


Frequency,  kHz 


Figure  1.  Absorption  coefficient  vs  frequency. 


Propagation  Loss  in  the  Near-Surface  Ducts 

The  near-surface  empirical  propagation-loss  equations  were  obtained  by  Marsh 
and  Schulkin  from  the  AMOS  transmission-loss  measurements  for  frequencies 
between  2 and  25  kHz.3  f or  completeness  the  equations  applicable  to  surface-duct 
propagation  loss  will  be  presented  here,  even  though  they  have  been  reproduced  in 
numerous  reports.  Certain  limits  to  their  applicability  have  been  imposed,  and  new 
subsurface  duct  equations  are  used.  The  equations  for  the  scattering  loss  coefficient 
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will  be  shown  first;  then  expressions  for  scaling  range  and  depth;  then  four 
propagation-loss  equations  from  the  AMOS  study,  as  well  as  a table  which  describes 
their  application  with  respect  to  the  scaled  ranges  and  depths;  then  the  normal  mode 
cutoff  term  for  lower  frequencies,  which  is  to  be  added  to  the  chosen  propagation- 
loss  expression;  and  lastly,  the  subsurface-duct  equations. 


Surface-Scattering-Loss  Coefficient.  For  energy  propagating  in  a surface  duct, 
the  surface-scattering  loss  per  bounce,  P.  is  obtained  from  an  analytic  expression 
developed  by  Schulkin  and  Marsh:7 


P=  10  log  1 1+C/7//4. 14)4 1 . /7/<4.2691 

= 1.59(/7t)1/2  , Jh> 4.2691 


(9) 


where 

/is  frequency  in  kHz 

and 

h is  waveheight  in  feet 

The  surface-scattering-loss  coefficient,^,  is  obtained  by  dividing  P by  the 
bounce  distance,  RB.  that  is 

= VfRH  dB/kyd  (10) 

where 

Rb  = (Cj  -cp1/2/  1500g0  kyd 

where 

Cy  and  are  the  sound  speeds  in  ft/sec  at  the  layer  depth  and  the  surface, 
respectively 

and 

g()  is  the  sound-speed  gradient  in  sec'1  in  the  surface  duct 
The  surface  reflection  loss  is 

Hs  = as{R-R')  \orR>R' 

= 0 for  R^R' 

where 

R is  the  horizontal  range  between  source  and  target,  in  kyd 
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where 


ZL  is  the  layer  depth  in  feet 

and 

r,  is  a scaled  range  defined  in  the  next  section 

A curve  of  surface  loss  versus  the  product  of  frequency  and  wave  height  is 
illustrated  in  Fig.  2. 


figure  2.  Surface-reflection-loss  curve  (from  Ref.  7). 


AMOS  Scaling  Relationships  and  Propagation-Loss  Equations.  The  AMOS  equa- 
tions3 may  be  expressed  in  terms  of  the  following  variables: 

R horizontal  range  between  source  and  target  (kyd) 

/-! : surface  layer  depth  (feet) 

Z.y.  source  depth  (feet) 

Zf:  target  depth  (feet) 

/:  frequency  (kHz) 

The  variables  R,  Zr  ■ Zx . and  ZT  are  used  to  define  scaled  variables  which  in  turn  are 
used  in  four  surface-duct  propagation-loss  equations.  The  selection  of  the  appropriate 
propagation-loss  equation  is  determined  by  scaled  zones  and  the  scaled  source  and 


target  depths.  Three  propagation  zones  are  distinguished,  a zone  of  direct  radiation,  a 
zone  of  single-surface-bounce  scattering,  and  a zone  of  multiple-surface-bounce  scat- 
tering, as  defined  in  Table  2. 


Table  2.  Selection  of  the  proper  AMOS  propagation-loss  equation. 


Insonification  Zone 
(Sealed  Range) 

Source-Target  Positions 
(Scaled  Depths) 

Propagation-Loss 
Lquation  Used 

0 ^ r 

zx  ^ ! ami  Zj  K 1 

n = iix 

z \ and/or  Zj-^>  1 * 

//  = Mini//,  ,//.,)** 

r,  <r  + 1/2 

ZX'  value* 

//  = Mini//,. //.,)** 

r,  + 1/2  < r 

c^..  z^-any  value* 

//  = Min  i//,y/4)** 

*Thc  AMOS  equations  are  used  in  this  model  only  if//  1000  feet 

and  if 


/.x  and/r<  32y/7J 

Thus,  the  source  or  target  may  be  located  at  30\/2i  . the  traditional  best 
depth  to  avoid  detection,  but  not  much  deeper  for  the  AMOS  equations  to  be 
used. 

♦•The  minimum  of  the  quantities  contained  within  the  parentheses  is 
retained. 

The  scaled  variables  are 

r = rvt; 

zx  = VZX/Z, 
zT  = v/Z7-/Z/ 

and 

r,  = ( 1 -zx)/4  + ( 1 -zT)/ 4,  zx<\,  rr<  1 
= ( 1 -zx  )/4  + ( Vz^\  )/5 , zx  < 1 , zT>  1 
= (V^d)/5  +(l-sr)/4,sjr>l.r7.<l 

The  four  propagation  loss  equations  are 

//,  = 20  log  ( R/R() ) + aR  + (z/r,  )G(zx z:r)  + 00 


//,  = 20  log  </?/*„)  + aR  + (25-  V \ZX-Z,\  - V\Z^Z[\  + ?/?)(// 25)1/3  + 60 
^3  = 20  log  </?//?,,>  + a/?  + 2(r-r,)F(zx,z7.)  + | 1 -2<r— /-j  > | G(:x.zf  ) + 60 

//4  = 10  log  {R/R^)+  {ot+ots)R  + FiZyZj) 

- a^Zl  (r+0.5 ) + 1 0 log  | \"rZ,  (r+0.5  )1  + 60 


where 


R{)  is  a reference  range  of  1 yard 

a is  the  absorption  coefficient  given  by  Hq.  8 

a is  the  surface  reflection-loss  coefficient  given  by  Hq.  0 and  10 

G(zx,zt)  = 0.1X102Jl*jr*rlC/'/25),/3,|r^-«7.|<l 


20(//25)l/3, 


U’rtl  > ' 


F(zx,zr)  = 0.4C/10"*  + 10'r  + 1 0 'X~‘T  ) 


c,  = \. 


= 0.5/i/3,/>8 

One  of  the  equations  //, , //2,  //3.  or  //4  will  be  used  for  1 1 in  Table  I . depending  on 
the  relationships  between  the  scaled  parameters  r.  r{ , zx,  and  zT.  Table  2 shows  the 
selection  of  the  proper  equation. 

Modification  of  the  AMOS  Equations  for  Lower  Frequencies.  The  AMOS 
surface-duct  propagation-loss  model  discussed  above  is  empirical  and  not  valid  for  fre- 
quencies less  than  2.2  kHz.  Predictions  based  on  the  AMOS  model  for  frequencies 
below  2.2  kHz  typically  do  not  show  good  agreement  with  Fleet  operational  data. 
This  disagreement  with  low-frequency  data  can  be  explained  in  part  by  the  low- 
frequency  cutoff  phenomenon  as  described  by  Arase8  and  ( lay.9  In  the  paper  by 
Arase,  experimental  attenuation  rate  is  studied  as  a function  of  frequency.  For  fre- 
quencies less  than  cutoff,  the  attenuation  rate  increases  as  a function  of  decreasing 
frequency.  Arase  found  good  agreement  between  experimental  attenuation  rate  and 
theoretical  first-mode  attenuation  rates  based  on  a normal-mode  propagation  model. 


To  extend  the  applicability  of  the  AMOS  surface-duct  propagation  model  to 
lower  frequencies,  an  additive  cutoff  loss  term  has  been  incorporated.  It  is  based  on 
an  approximation  to  the  normal-mode  surface  duct  model  of  Pedersen  and 
Gordon10,11  The  normal-mode  cutoff  term  is 

Hco  = 8686r,  R 

where 

R is  the  horizontal  range  in  kyd 

and 

t{  is  the  first-mode  damping  factor,  which  is  given  by 
r,  =(nJg^)'P  lm(MXl )/Cs 

where 

/is  frequency  in  kHz 

g0  is  the  sound-speed  gradient  in  the  duct  in  sec-1 
//u(A/A'1 ) is  the  imaginary  component  of  the  first  eigenvalue 

Propagation  Loss  in  the  Subsurface  Duct.  A subsurface  duct  is  formed  when  the 
sound-speed  profile  has  a concavity  to  the  right,  as  illustrated  in  Fig.  3a.  The  dotted 
lines  in  Fig.  3a  extend  from  the  upper  to  the  lower  boundaries  of  the  several  ducts 
pictured.  If  both  the  source  and  the  target  depths  are  between  these  boundary  depths, 
then  the  subsurface-duct  equations  are  used,  except  for  two  special  cases:  ( 1 ) when 
the  subsurface  duct  is  part  of  a surface  duct  (Fig.  3b),  and  (2)  when  the  axis  of  the 
duct  is  the  deep  minimum  (deep  axis)  of  the  sound-speed  profile  (Fig.  3c).  In  case  ( I ) 
the  AMOS  surface-duct  equations  are  used,  and  in  case  (2)  ray  tracing  is  used. 

Since  the  equation  for  subsurface-duct  propagation  loss  derived  from  the  AMOS 
data  gives  losses  consistently  greater  than  spherical  spreading  loss,  and  since  a duct  in 
\ the  ocean  is  expected  to  trap  sound  energy  and  give  losses  on  the  order  of  cylindrical 

' spreading,  the  subsurface-duct  equation  given  in  Ref.  3 was  modified.  The  new  equa- 

tions are 

lln  = 20  log  R + aR  + II.  +60.  R < R . 

(II) 

= 1 0 log  RRy  + a V?  + //,  + 60.  R > Ry 

where 

Hp  is  propagation  loss  in  the  subsurface  duct 
R is  range  in  kyd 


• •. . 


a is  the  absorption  coefficient  in  dB/kyd 
Hy  is  the  depth  loss  factor  (see  below) 

R is  the  maximum  distance  a ray  starting  on  the  duct  axis  can  travel  without 
vertexing  (derived  below) 

The  depth  loss  factor,  Hz.  is  a function  of  source  depth.  Zx , in  feet,  target  depth.  Z? . 
in  feet,  duct  axis  depth,  ZD.  in  feet,  duct  width.  Zw.  in  feet,  and  frequency./,  in  kHz. 
and  is  given  by 

Hy  = O.^/S^IO^1^'  + 10^I/Z*  + 10  IZt~ZdIIzW)  (12) 

where 

is  set  to  unity  for/<8  kHz 

The  quantity  Rv  in  Eq.  1 1 is  derived  in  the  following  manner.  Figure  4 depicts 
two  paths  which  have  vertices  (become  horizontal)  at  the  edges  of  the  duct,  an  up- 
ward path  and  a downward  path.  In  this  case  the  downward  path  is  used  to  deter- 
mine Rv,  since  it  produces  the  larger  value.  Ray  theory  for  a constant-gradient  sound- 
speed  profile  segment  yields  the  following  equation  (or  Rv  : 


R = 

’ hrfti 


/c2  -c2 

M I) 


where 


CM  is  the  sound  speed  at  the  boundaries  of  the  duct 
CD  is  the  sound  speed  at  the  duct  axis 


C -C 

yh  = -7^- — is  the  sound-speed  gradient  in  the  lower  branch  of  the  duct 
Propagation  by  Ray  Acoustics 

The  Wave  Equation.  The  propagation  of  acoustic  energy  through  the  ocean 
medium  is  approximated  by  the  wave  equation: 

_ 2 1 5 2p 

v p = IT  . 2 

C~  <5  r 


figure  4.  Ray  vertexing  ranges  in  a subsurface  duct. 


where 

p is  the  excess  pressure  at  some  constant  phase  on  the  wave  front 

and 

C is  the  sound  speed:  both  p and  C are  functions  of  the  spatial  and  time 
dimensions 

The  Eikonal  Equation.  Wave  fronts  of  harmonic  waves  which  are  defined  by 


W(x,y,z)  = C’0  U-t0 ) 


(14) 


can  represent  a time-independent  solution  to  the  wave  equation  called  the  eikonal 
equation,  which  is 


(f )!  ♦ (f*)‘  * (r)2  - (¥ 


The  eikonal  equation  can  be  derived  by  relating  the  pressure  wave  front  to  W and  by 
directly  substituting  into  Eq.  13  or  by  using  Eq.  14  to  apply  Huygens'  principle  for 
constructing  successive  wave  fronts  normal  to  preceding  wave  fronts.1" 

The  eikonal  equation  is  a special  case  of  the  wave  equation,  valid  only  if  the 
wavelength  of  the  propagating  sound  is  much  less  than  the  radius  of  curvature  of  any 
portion  of  the  ray  path. 


i 
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Applicability  for  Sonar  Systems.  All  current  active  sonar  systems  operate  at  fre- 
quencies large  enough  (wavelengths  small  enough)  for  the  eikonai  equation  to  be  valid. 
Passive  sonar  systems,  however,  are  designed  for  detection  at  low  frequencies.  At  500 
Hz  the  wavelength  is  about  10  feet,  which  does  not  preclude  use  of  the  eikonai  equa- 
tion. At  100  Hz  the  wavelength  is  about  50  feet,  which  suggests  that  the  applicability 
of  the  eikonai  equation  and  the  ray  acoustics  which  spring  from  it  may  not  be  valid 
for  realistic  sound-speed  changes  in  the  ocean. 

Definition  of  a Ray  Path.  Normals  to  the  wave  surface  (constant  phase)  in  space 
at  a given  time  can  be  constructed  using  Huygens’  principle  or  can  be  computed  using 
the  eikonai  equation,  Eq.  1 5.  The  continuous  path  of  a normal  to  the  wave  surface 
obtained  by  varying  time  is  called  a ray  path.  It  is  along  ray  paths  that  acoustic 
energy  theoretically  propagates. 

Since  sound  speed  is  assumed  to  be  a function  of  depth  only,  ray  paths  are  re- 
stricted to  the  vertical  plane  through  the  source  and  target  positions.  Thus,  ray  tracing 
is  accomplished  in  two  dimensions  which  are  (horizontal I range,  R.  and  depth.  Z. 

The  origin  is  at  the  ocean  surface  vertically  above  the  source,  and  Z is  positive  down- 
ward, as  in  Fig.  5.  Positive  angles  are  measured  downward  from  the  horizontal. 


figure  5.  Coordinate  geometry  for  ray  tracing. 
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Propagation  Loss  in  Ray  Theory.  Propagation  loss  is  defined  as  the  ratio  of  the 
source  strength  to  the  sound-field  intensity  in  decibels.12 


U = 10  log  j 


( 16) 


where 

f is  the  source  strength  in  energy  units  per  unit  solid  angle  per  second 

and 

/ is  the  intensity  in  units  of  energy  crossing  an  elemental  area  JA  perpendicular 
to  the  ray  path  and  bounded  by  an  elemental  solid  angle 

The  ratio  Fjl  is  easily  obtained  using  the  concept  of  solid  angle: 12 

/ c/n 

F (IA 


with  the  absolute  value  sign  used,  since  neither  F nor  / can  ever  be  negative.  With 
cos  <f>x  =CX/C\  . f.q.  16  becomes 


//=  10  log 


+ 60  in  dB  re  i yd 


(17) 


where 

II  is  propagation  loss  in  decibels  re  0 dB  at  1 yard  slant  range  from  the  source 
along  the  ray  path 

R is  range  in  kyd 

<t>T  is  the  angle  of  the  ray  at  the  ending  (target)  depth 
Cj,  is  the  sound  speed  at  the  source  depth 
Cj,  is  the  sound  speed  at  the  vertex  depth 

and 

60  represents  the  conversion  in  units  from  yards  to  kiloyards  in  the  range  and 
range-derivative  factors 

4>x  is  the  starting  angle  of  the  ray  at  the  source  depth 

is  the  constant-depth  derivative  of  the  range  with  respect  to  the  starting 
angle 
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Except  for  the  units,  Eq.  1 7 is  general,  independent  of  algorithms  for  determining  R 


and 


for  which  expressions  will  be  given  in  Eq.  20,  21,  23,  and  24. 


Curvilinear  Sound-Speed  Profile.  The  adjectives  “straight-line”  and  “curvilinear,” 
used  as  descriptors  of  sound-speed  profile,  signify  the  functional  representation  of  the 
sound  speed  from  the  surface  to  the  bottom  of  the  ocean  when  sound  speeds  are 
known  only  for  discrete  depths.  “Straight-line”  refers  to  a sound-speed  profile  con- 
structed of  straight-line  segments  connecting  the  points.  “Curvilinear”  refers  to  fitting 
the  sound-speed  profile  to  within  a given  tolerance  with  quadratic  equations  of  the 
form 


C = cl\\-Ka(Z-Za)2)-'  (18) 

where 

Z is  depth 

C , Za,  and  Ka  are  three  curve  parameters  evaluated  for  each  curve  segment 


The  curve  parameters  between  successive  segments  are  dynamically  adjusted  so  that 
successive  curves  and  their  derivatives  are  continuous  at  the  sound-speed  profile  points. 

Pedersen13  has  shown  that  a straight-line  sound-speed  profile  produces  undefined 
propagation  losses  for  certain  ranges  calculated  using  ray  theory.  This  problem  is  a 
direct  result  of  having  discontinuities  in  the  first  derivative  anywhere  on  the  sound- 
speed  profile.  Introducing  more  points  compounds  the  problem,  since  more  discon- 
tinuities in  the  first  derivative  are  added. 


Pedersen  and  Gordon14  introduced  the  curvilinear  approximation  to  the  sound- 
speed  profile,  Eq.  18,  and  required  that  the  curve-fitting  segments  and  their  deriva- 
tives be  everywhere  continuous  over  the  profile.  Numerical  examples  given  by 
Pedersen  and  Gordon  show  that  a curvilinear  approximation  to  the  sound-speed 
profile  eliminates  the  false  caustics  (false  regions  of  infinite  intensity)  characteristic 
of  the  straight-line  approximation. 

The  major  obstacle  to  implementing  the  curvilinear  method  is  developing  a com- 
puter program  which  can  reliably  fit  the  sound-speed  profile  within  an  adequate 
tolerance.  Some  computer  programs  were  written  and  are  still  being  used  which  re- 
quire trial  and  error  on  the  part  of  the  user  to  adjust  and  add  points  to  the  sound- 
speed  profile  in  order  to  get  an  adequate  fit.  This  kind  of  software  is  unacceptable  for 
support  of  an  active  sonar  performance  prediction  program  which  must  run  to  com- 
pletion without  idiosyncrasies.  Fortunately,  a set  of  subroutines  does  exist  which 
performs  the  task  adequately.  These  subroutines  were  developed  by  Charles 
Bartberger15  of  NADC.  He  used  Pedersen’s  and  Gordon’s  equations  and  applied  vary- 
ing tolerance  versus  depth.  The  subroutines  used  in  the  LORA  program  are  essentially 


IX 


Bartberger’s,  with  changes  in  one  of  the  tolerance  equations  to  obtain  a closer  fit.  The 
new  tolerance  equations,  in  the  form  of  Bartberger’s  equations,  are 


T = Tj(yJ  + ft/sec 


where 

ym,  is  the  average  depth  of  all  the  points  being  fitted 

and 

8 is  the  smallest  of  the  absolute  straight-line  gradients  connecting  adjacent 
points  of  the  set,  with  the  functions 

Tj{yav)  = 2.9  - 2.5  exp  (-0.00000272^) 


and 


Tg(g)  ="T  I 1 + tanh  4.8(g- 0.35)1 

In  a numerical  evaluation  with  a typical  surface-duct  sound-speed  profile  which 
extends  to  10.000  teet,  the  curve-fitted  profile  differs  from  the  linear  profile  by  less 
than  0.2  ft/sec  at  the  duct  maximum,  less  than  0.5  ft/sec  down  to  1 .000  feet,  and 
less  than  0.01  ft/sec  between  4.000  and  10.000  teet.  The  standard  deviation  is  less 
than  0.2  ft/sec  over  the  whole  profile  (110  points  in  the  sample). 

Equation  1 8 can  manifest  three  types  of  curves,  depending  on  the  signs  of  the 
curve  parameters.  These  are  shown  in  Fig.  6.  The  sound-speed  profile  will  be  con- 
structed from  segments  which  appear  in  these  forms.  Curves  of  Types  1 and  11  occur 
most  frequently.  Type  111  seldom.  Even  though  a curve  of  Type  1 or  11  may  fit  a 
minimum  or  maximum  adequately,  an  extra  layer  is  introduced  at  that  point.  Refer- 
ence 15  gives  both  a complete  description  of  the  mathematical  technique  for  curve- 
titting  the  sound-speed  profile  and  a listing  of  the  computer  software. 

Earth  Curvature  Correction.  Bartberger,s  gives  a derivation  for  the  first-order 
spherical  earth  curvature  correction  to  the  sound-speed  profile.  With  slight  changes  in 
nomenclature  his  equations  are 

C]  = c/>  i + z;irc) 

zt  =z;(\  + z;  1 2R,) 


where 

<7  is  input  sound  speed 


Type  I:  C2  > 0,  K >0 
a a 


Type  II  C2>0,K  <0  Type  III:  (?  < 0.  K >0 

a a a a 


C 


figure  6.  Schematic  plots  of  the  three  types  of 
curves  used  to  fit  the  sound-speed  profile  with 
curvilinear  segments  (from  Ref.  I4>. 


Z{  is  corresponding  input  depth 

Cj  and  Z-  are  the  corrected  sound  speeds  and  depths 

and 

R is  the  radius  of  the  earth  in  corresponding  units 

Snell’s  Law.  This  model  assumes  that  the  sound  speed  varies  in  only  one  spatial 
dimension,  depth.  Snell’s  law  is  therefore 

cos  0 _ C ( j 9 

cos  <t>f)  C0 


where 

tf>  and  0()  are  angles  that  a ray  path  makes  with  respect  to  the  horizontal 

and 

C and  C0  are  sound  speeds  that  correspond  to  the  depths  at  which  the  angles 
are  measured 

Snell's  law  applied  to  the  sound-speed  profile  can  yield  the  following  information: 

1 . The  starting  angle  for  a ray  to  get  to  a given  depth. 

2.  The  angle  of  the  ray  at  any  depth,  given  the  starting  angle  and  depth. 

3.  Determination  of  whether  it  is  possible  for  a ray  to  vertex  at  a given  depth, 
given  the  starting  angle  and  depth. 

4.  The  depths  that  a ray  cannot  reach,  given  the  starting  angle  and  depth. 


i 


I 
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Snell’s  law  is  an  important  concept  that  has  been  used  extensively  in  this  model, 
both  for  the  selection  of  rays  to  be  traced  and  in  the  interpretation  of  the  outputs  in 
anticipating  modes  of  propagation,  shadow  zones,  and  the  like.  Frequently  users’ 
questions  about  the  interpretation  of  the  LORA  program  outputs  are  answered 
through  some  application  of  Snell’s  law. 

Curvilinear  Ray  Tracing.  Curvilinear  ray  tracing  means  calculating  the  locus  of 
points  in  range  and  depth  of  a ray  governed  by  a sound-speed  profile  which  is  defined 
by  curvilinear  segments  of  the  form  of  Hq.  18.  Four  important  quantities  are  obtained 
from  ray  tracing  between  specified  depths:  range,  R,  time,  T,  range  derivative  with 
respect  to  starting  angle,  D,  and  absorption  loss,  //( . 

A layer  is  defined  as  the  depth  region  over  which  one  given  curvilinear  segment 
is  valid.  An  interface  is  defined  as  the  depth  at  which  one  curvilinear  segment  makes 
transposition  to  another.  Figure  7 shows  a typical  layer.  In  ray  tracing,  two  depths 
are  specified  and  the  ray  parameters  R.  T.  D,  and  //,  are  calculated  and  summed  as 
the  tracing  proceeds  through  the  layers.  //,  will  be  described  in  the  next  section. 
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Sound  Speed 


Range 


I igure  7.  Layer  geometry. 


Let  C Kaj,  and  Zai  be  the  segment  parameters  for  the  curvilinear  sound-speed 
profile  layer  i between  interfaces  i and  / + 1 . Let  the  upper  and  lower  depths  tor 
layer  i be  Z,(  and  Z,.,  respectively.  (Notice  that  Zu  = Z(  and  Z2i  + Z(+|  .1  Let  Cj,  be  the 
vertex  sound  speed^the  maximum  sound  speed  the  ray  can  achieve,  the  sound  speed 
at  which  the  ray  becomes  horizontal.  Let  </>j  be  the  angle  which  the  ray  makes  with 
the  horizontal  (downward  is  positive).  Pedersen  and  Gordon14  detine  nine  auxiliary 


quantities: 


A'  ** 

A = + C — — 

Ai  ~ S i 


A = cj(c2-c2) 


a;  = (c~  + c 4 yc2 

RU  = Hzv  rzJ 

Rn  ~ A^zn  -Zaj) 

B\i  = cot  ^zu~zJ 

% = cot  oyz^-zj 

T\i  = tan  Hz\  rzat 

Ty  = tan  <t>i(Z2i-Za^ 

The  ray  parameters  /?..  /j.  and  //  for  a ray  passing  completely  through  layer  i 
can  now  be  expressed  in  terms  of  the  above  nine  quantities,  l or  a Type  I curvilinear 
segment  (see  Fig.  6)  the  range  in  layer  i is 

R,  = j cos'1  \Dt  | RVRV\  + tan  0,.  tan  0(+,  I ( 20 


and  for  Types  II  and  111 


j_.  rKl+lan^ 

' " Ai  11  L 1^1  +tan0(+,_ 


The  travel  time  in  layer  i is 


7 ; = 2T  Wi'Tu  + T*> 

V 

The  derivative  of  range  with  respect  to  vertex  sound  speed  is 
(SR.  \ 

C>\8c)  = = 


The  derivative  with  constant  Z indicates  ray  tracing  is  between  two  constant  depths. 
Using  Snell’s  law 


(i& --*(!& 


where 


<(>x  is  the  starting  angle  of  the  ray  at  the  transmitter 


W\V/ 


tan  </> j.  tan  <pY 

-c^  \ = ~c~  B> 


If  a ray  forms  an  apex  in  layer/.  Ru  is  replaced  by  i£>  r1/2 . tan  0(.  is  set  to  zero, 
and  Eq.  20  and  23  are  still  applicable.  If  a ray  forms  a nadir  in  the  layer.  /?2,  is  re- 
placed by  |Z}|'/2  and  tan  <pi+i  is  set  to  zero.  Pedersen  and  Gordon14  list  problems  and 
solutions  for  cases  when  a vertex  sound  speed  is  at  an  interface.  The  LORA  program 
selects  vertex  sound  speeds  in  such  a way  that  these  problems  will  never  occur;  how- 
ever, the  ray  tracing  algorithm  is  written  to  handle  them  according  to  the  specifica- 
tions of  Pedersen  and  Gordon. 

The  total  range.  R.  time.  T.  and  range  derivative.  D.  are  obtained  by  summing 
the  appropriate  quantities  calculated  for  each  layer  over  the  total  number  of  layers 
traversed  from  one  end  of  the  ray  path  to  the  other.  Thus,  some  layers  contribute 
more  than  once  in  the  summation.  Using  Eq.  20-23  and  converting  from  teet  to 
kiloyards 

R = 2 R(./3000  kyd  \ 


T = L 7'  see 


D = = (tan 0V/3()O()  (;.)  L kyd 


23 


i 


R and  D are  used  directly  in  Eq.  17  for  the  calculation  of  propagation  loss.  T is  to  be 
used  to  obtain  the  reverberation  level  at  the  time  at  which  target  echo  arrives  at  the 


Absorption  Loss  Along  the  Ray.  The  path  length  in  the  z'th  layer  is  approxi- 
mated by  the  slant  range 


S’.  = [Rr+(Zi+ , -Z,.)2l1/2/3000  kyd 


where 


Rj  is  the  range  in  feet  given  by  Eq.  20  and  21 
Z is  the  upper  interface  depth  in  feet 

and 

ZJ+|  is  the  lower  interface  depth  in  feet 
The  total  absorption  loss  is 


IlA  = S v,. 


where 


u(.  is  the  absorption  loss  in  dB/kyd  calculated  for  the  mid-depth  at  each  layer  of 
the  sound-speed  profile  after  curve-fitting 


£ indicates  that  the  summation  is  over  all  the  layers  a given  ray  path  enters 
1 (some  layers  may  have  more  than  one  contribution) 

This  value  for  //(  is  used  in  Table  1 in  all  propagation  modes  except  the  duct  modes. 

The  absorption  loss  is  calculated  by  using  Eq.  8.  The  pressure  is  calculated  for 
the  mid-depth  of  each  layer,  Zm , and  the  relaxation  frequency  is  calculated  for  the 
temperature,  Tm , at  Zm.  The  temperature  is  given  by 

Tm  = /I,  - (0.000434  Zm)2/A{  \ 


A,  = 122.83-4.6676  V 5132- Cm  +0.0 165  Zm 


L 


where 

Zm  is  the  mid-depth  of  the  layer 

and 

C’m  is  the  corresponding  sound  speed 

The  form  of  Eq.  25  results  from  an  approximation  to  Wilson’s  equation,16 
inverted  to  obtain  temperature  by  the  regula  falsi  method.  It  is  understood  that 
Wilson’s  equation  is  properly  a least-squares  polynomial  fit  to  the  regression  variables 
temperature,  pressure,  and  salinity,  and  that  inversion  to  obtain  temperature  as  a func- 
tion of  sound  speed,  pressure,  and  salinity  produces  a larger  standard  deviation  than  if 
a new  regression  were  performed  on  the  original  data.  However,  absorption  loss  is  not 
severely  influenced  by  temperature,  and  the  errors  introduced  by  inverting  Wilson’s 
equation  and  by  using  Eq.  25  are  negligible  for  the  total  propagation  loss  along  a ray 
path. 


Equation  25  is  a trial  and  error  fit  as  to  form,  but  some  of  the  coefficients  have 
been  adjusted  to  produce  a minimum  standard  deviation  in  temperature  prediction. 

H.  W.  Frye’s  regression  analysis  program17  was  used  to  fit  polynomials  to  a large  set 
of  temperatures  obtained  by  inverting  Wilson’s  equation  for  specific  sound  speeds, 
pressures  (depths),  and  salinities.  Realistic  values  for  ocean  parameters  produced  a 
standard  deviation  of  less  than  0.1°F.  McGirr  and  Hall18  have  shown  that  an  error  of 
even  l°F  would  not  seriously  debilitate  the  prediction  for  the  absorption  coefficient. 

Surface  Loss.  Loss  caused  by  a ray  reflecting  from  the  ocean  surface  is  a model 
input.  The  user  of  the  LORA  computer  program  can.  as  an  option,  use  the  value  for 
surface  loss  per  bounce  calculated  by  Eq. c).  The  total  loss  caused  by  surface  reflections 
is  the  number  of  surface  bounces  a ray  path  undergoes  times  the  surface  loss  per 
bounce.  This  loss  is  added  to  the  other  components  of  propagation  loss  for  one  of  the 
modes  of  propagation  described  in  Table  1 . 

Bottom  Loss.  Loss  caused  by  a ray  bouncing  off  the  ocean  bottom  is  a model 
input.  In  the  computer  program  it  is  calculated  by  linearly  interpolating  a user- 
supplied  piecewise  linear  array  of  bottom  loss  values  versus  grazing  angle.  The  bottom 
loss  per  bounce  is  multiplied  by  the  number  of  bottom  bounces  to  obtain  total  bot- 
tom loss,  which  is  used  in  the  equation  for  bottom  bounce  propagation  in  Table  1 . 

Ray  Paths.  Vertex  sound  speed,  the  maximum  sound  speed  a ray  can  obtain,  is 
given  by  Snell’s  law.  At  this  sound  speed,  the  ray  becomes  horizontal  and  thereafter 
traces  a mirror-image  path  about  the  vertical.  With  the  ending  angle  0 = 0.  Eq.  Id 
gives  the  expression  for  vertex  sound  speed: 


25 


where 


-1.  . 


<j>0  and  C0  are  ray  angle  and  sound  speed,  respectively,  at  the  same  depth 

The  sound-speed  profile,  C(Z).  is  a single-valued  function  of  depth.  The  inverse 
of  C(Z),  depth  as  a function  of  sound  speed,  is  multivalued.  Thus,  given  Cv,  a number 
of  values  for  the  vertex  depth,  Z,.,  may  be  obtained.  Given  a sound-speed  profile  and 
the  depth  and  angle  at  which  a ray  starts,  one  upper  and  one  lower  vertex  depth  can 
be  found  by  scanning  the  profile  upward  and  downward  from  the  starting  depth  and 
finding  the  depths  (upward-downward)  at  which  the  ray  achieves  the  sound  speed  CJ . 
The  ray  is  confined  between  these  upper  and  lower  vertex  depths,  and  mirror-image 
paths  are  produced  as  the  ray  travels  on.  In  some  cases  Cj,  is  greater  than  any  sound 
speed  found  in  the  upward-downward  scan  of  the  sound-speed  profile;  this  indicates 
that  the  ray  will  strike  the  surface  or  bottom  of  the  ocean.  In  this  model  these 
boundaries  provide  mirror  images  of  the  path  with  range  and  give  the  same  replication 
of  path  as  if  the  ray  had  gone  through  a refractive  turning  point.  Once  a ray  is  traced 
from  its  upper  vertex  depth  (or  the  surface)  to  its  lower  vertex  depth  (or  the  bottom), 
one-half  of  a cycle  of  the  ray  is  obtained  (in  range  versus  depth).  The  second  half  of 
the  cycle  is  obtained  by  annexing  the  mirror  image  to  the  first  half.  As  many  cycles 
as  desired  can  be  appended  to  extend  the  horizontal  range  of  the  ray.  This  model 
uses  a maximum  of  30  cycles  (up  to  30  cycles  appended  to  points  on  the  basic  ray 
path  in  the  first  cycle)  for  both  active  sonar  two-way  propagation  loss  and  passive 
one-way  propagation  loss.  The  program  user  can  specify  the  maximum  number  of 
refractive  cycles  or  bottom  bounce  cycles  to  be  used. 

This  model  is  concerned  with  tracing  rays  between  several  depths,  computing 
propagation  loss  for  the  resultant  ranges,  and  using  these  losses  to  determine  the  sig- 
nal return  from  a target  or  the  reverberation  from  surface,  bottom,  or  volume 
acoustic  scatterers.  Thus,  rays  are  traced  from  transmitter  depth  to  target  depth,  as 
well  as  to  the  surface  and  bottom.  When  a ray  is  traced  from  one  depth  to  another, 
five  ray  parameters  are  calculated  for  later  use  in  the  calculation  of  propagation  loss 
and  for  interpolation:  absorption  loss.  //,, , range.  R.  time.  T.  derivative  of  range  with 
respect  to  starting  angle,  D.  and  boundary  loss,  II.  For  convenience  in  subscripting.  A 
will  be  used  for  HA  in  the  following  discussion.  The  parameters  A.  R.  T.  D.  and  /i  are 
associative  from  segment  to  segment  on  a ray  path,  so  that  several  segments  can  be 
annexed  together  by  adding  the  corresponding  components.  The  parameters  for  a path 
segment  can  be  thought  of  as  forming  a vector  space 

F=  (A.  R.  T.  D.  II) 


where  the  linear  combination  of  the  components  is  the  only  property  ofjnterest  in 
this  study  (scalar  product  and  norm  have  no  physical  significajice  here).  1’  will  be 
used  as  a shorthhand  notation  for  its  components,  and  when  V is  subscripted,  the 
subscript  applies  to  each  of  the  components.  An  example  of  V for  a ray  traced  from 
depth  a to  depth  b is 

^ R.te  Tah.  Dah.  Bah) 


2b 


p 

i 


One  of  the  many  possible  range  components,  Rah,  is  depicted  in  Fig.  8. 

Because  of  the  cyclic  nature  of  the  ray  path,  there  are  an  infinite  number  of 
vectors  Vab  which  can  connect  two  depths.  Thus,  following  the  dotted  ray  path  in 
Fig.  8,  one  can  see  that  the  path  will  cross  depths  a and  h infinitely.  A ray  path  con- 
necting two  depths  which  does  not  contain  cyclic  components  or  mirror  image  com- 
ponents about  some  vertical  (range  = constant)  is  defined  as  a ray  path  segment. 

In  computing  propagation  loss  from  transmitter  to  target,  it  is  only  necessary  to 
trace  three  segments:  (1)  upper  vertex  depth  or  surface-to-transmitter  depth, 

(2)  transmitter  depth  to  target  depth,  and  (3)  target  depth  to  lower  vertex  or  bottom 
depth.*  These  segments  are  depicted  in  Fig.  9.  with  one  of  the  vector  components.  R. 
the  horizontal  range,  also  shown.  The  vectors  resulting  from  calculating  the  segments 
are 
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I ijiure  8.  Range  component  and  range  cycle. 


* These  segments  apply  when  transmitter  depth  is  less  than  the  target  depth.  It  the  transmitter  depth  equals 
the  target  depth,  then  segment  2 is  null.  If  the  transmitter  and  target  depths  are  reversed,  then  “transmitter”  and 
“target”  designations  in  the  description  of  the  three  segments  should  be  reversed. 
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l igure  9.  Path  segments. 


where  the  subscripts  imply  the  following: 

SX  -*  surface  to  transmitter 

XT  -*  transmitter  to  target  I 

and 

TB  -*■  target  to  bottom 

* 

The  B components  are 

♦ 

Bs v = 1/2  X surface-loss-per-bounce  or  zero  if  the  ray  does  not  strike 
the  surface 

bxt  = 0 j 

and 

B,h  = 1/2  X bottom-loss-per-hounce  or  zero  if  the  ray  does  not  strike  j 

the  bottom 

Once  the  path  segments  are  determined  for  a given  starting  angle  at  the  trans- 
mitter, they  can  be  combined  through  addition  to  give  four  basic  path  types.  These 
basic  path  types  are  illustrated  in  Tig.  10. 


28 


Depth 


If 

i 


Figure  10.  l our  basic  path  types  connecting 
transmitter  and  target  depths. 
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The  vector  representation  for  the  basic  path  types  is 


AT 


v = ~>V  + V 

' 2 ~rSX  y XT 


V = V + TV 
'XT  ~yTB 


v = ~>y  + V + ~>v 

y4  “ ' SX  r XT  ~v  TB 


(27) 


These  four  basic  path  types  are  the  only  ways  in  which  a ray  with  angle  ±0vat  the 
transmitter  can  travel  from  the  source  to  the  target  depth  without  completing  a cycle. 
The  vector  describing  a cycle  is 


Vc  = 2(KSv  + K^  + Kra) 


(28) 


All  possible  paths  and  their  parameters  from  the  source  to  the  target  may  now  be 
described  by 


T/i 


= vk  +nvc,  k = 1, 2.3 .4;  n = 1,2.3, 


(29) 


I 

t 

I 


where 

Vk  = is  the  vector  whose  components  are  absorption  loss,  range,  time,  range 
derivative,  and  boundary  loss* 

Equations  27-29  define  every  possible  path  connecting  transmitter  and  target 
depths,  given  that  the  complete  ray  path  has  a maximum  sound  speed  C„. 

Selection  of  Starting  Angles.  The  ray  paths  of  the  previous  section  all  had  the 
same  absolute  starting  angle  at  the  source.  (px . Given  the  maximum  sound  speed  C,,, 
<px  is  given  by  Snell’s  law,  Eq.  18,  as 

<px  = ± cos*1  ( CxtCx, ) 

This  model  considers  two  categories  of  rays,  those  which  ultimately  strike  the 
bottom  and  those  which  do  not.  The  former  category  is  called  the  bottom  bounce 
(BB)  propagation  mode,  and  the  latter,  the  convergence  zone  ( CZ )**  propagation 
mode.  For  BB  rays  (7,  exceeds  the  bottom  sound  speed.  C . For  CZ  rays  Cj.  is  less 
than  C„.  For  each  propagation  mode,  BB  and  CZ.  40  C’.'s  are  selected,  giving  a total 
of  80  unique  0v’s.  By  using  both  positive  and  negative  values  for^v.  160  ray  paths 
are  obtained. 


The  40  BB  C„’s  are  selected  by  using  the  equation 


cos  !</>„,-  ! ! -o  x p ( / / 4 1 1 ( 0,„  - 0,,’i 


where 

cj1  is  the  / th  vertexing  sound  speed  to  be  used  in  ray  tracing 

C.  is  the  sound  speed  at  the  source  depth 

<t>m  = 1 .553343033,  the  radian  value  for  89  degrees 

/ = 0, 1 39 

and 

<t>H  = cos*1  (CX/CB),  Cx  <CB 

= 0,  cx  > cB 

is  the  angle  in  radians  at  the  source  depth  of  the  ra\  which  just  grazes 
the  bottom 

•The  boundary-loss  component  is  one-half  the  actual  loss,  since  the  seements  containing  nonzero  bound- 
ary  losses  are  always  doubled  in  Eq.  27. 

**This  categorization  is  somewhat  arbitrary  in  that  refracted  and  surface  reflected  ray  paths  are  also 
included. 
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In  addition,  each  of  the  following  sound  speeds  is  included  in  the  set  if  it  is  greater 
than  both  Cx  and  C'B : 

C5±6 


where 

C'v  is  the  surface  sound  speed 

is  any  relative  maximum  in  the  sound-speed  profile 

and 

8 is  a small  constant  increment 

Selection  of  the  40  values  of  C\.  for  the  C'Z  propagation  mode  is  as  follows:  f irst, 
all  interface  sound  speeds  are  found  with  values  between  Cx  and  Cg.  Then  8 is  sub- 
tracted from  each  and  the  result  stored.  If  the  interface  sound  speed  is  a relative 
maximum,  8 is  added  and  the  result  is  included  in  the  set.  The  surface  sound  speed 
is  treated  in  the  same  way  as  a relative  maximum.  The  source  and  target  sound  speeds 
plus  8 are  included  in  the  set  Then  all  these  sound  speeds  are  ordered  in  increasing 
magnitude.  The  spacing  between  successive  sound  speeds  in  the  ordered  set  is 
examined,  and  more  sound  speeds,  approximately  evenly  spaced,  are  added  to  the  set 
to  fill  up  an  array  of  40  vertex  sound  speeds. 

The  BB  and  CZ  vertex  sound-speed  arrays  are  used  as  drivers  to  step  through  the 
starting  angles  at  the  source.  They  have  been  chosen  to  reduce  caustic  effects  and 
accurately  outline  the  shadow  zones.  These  sets  of  vertex  sound  speeds  also  produce 
roughly  evenly  spaced  range  values  between  successive  rays  for  average  sound-speed 
profiles  with  the  following  qualifications:  in  range  regions  with  highly  fluctuating 
propagation  losses,  the  rays  tend  to  be  spaced  very  close  together  (say.  0.1  kyd).  and 
in  range  regions  with  slowly  varying  propagation  losses,  the  rays  are  spaced  farther 
apart  (say,  2 kyd).  The  angular  coverage  is  ±84  degrees;  thus,  the  ocean  is  effectively 
filled  with  sound. 

Ray  Path  Interpolation.  Ray-tracing  methods  involve  specifying  a ray's  starting 
depth  and  starting  angle  (or  its  equivalent,  vertex  sound  speed)  and  following  the  ray 
through  the  medium  to  see  where  it  goes.  Range  and  range  derivative  are  obtained, 
from  which  propagation  loss  can  be  calculated.  Range  and  range  derivative  are  outputs 
of  the  ray-tracing  algorithm.  But  for  exceptional  cases,  a specific  range  for  a given 
depth  cannot  be  produced  except  by  iterative  techniques  or  interpolation.  The  tech- 
nique for  producing  propagation  loss  for  specific  ranges  becomes  a part  of  the  model. 
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For  interpolating  propagation  loss  between  adjacent  rays  of  the  same  class,  this 
model  uses  an  interpolation  technique  whose  basic  assumption  is  that 


//  = 20  log  {RIR..)  + IL 


where 


//  is  propagation  loss  at  range  R 


//0  is  a reference  propagation  loss  at  a reference  range  Ru 


Two  adjacent  rays  at  a given  depth  produce  the  propagation  losses  //  and  Hh  at  the 
ranges  R0  and  Rh . Equation  30  can  then  be  solved  for  Ru  and  //()  to  produce  the 
following  equation  for  If  : 


H / 20  /20 

//  = 20  log  |<  I-/-')  10  a + F 10  b \ 


where 


Xa-*b  ) 

The  computer  program  interpolates  in  1-kyd  increments.  Thus,  interpolation  is 
made  for  propagation  loss  for  all  integer  range  values  between  Ra  and  Rh.  Correspond- 
ing times  are  obtained  by  linear  interpolation.  These  times  are  used  later  for  interpo- 
lating reverberation  arrays. 

Interpolation  between  successive  rays  means  that  the  two  rays  have  successive 
vertexing  sound  speeds  in  one  of  the  arrays:  these  arrays  are  arranged  either  in 
increasing  magnitude  (for  convergence-zone  loss)  or  decreasing  magnitude  (for 
bottom-bounce  loss).  Also,  the  ray  paths  must  be  of  the  same  class:  both  ray  paths 
must  have  the  same  value  of  k and  n for  the  Vkn  defined  in  Eq.  29.  In  Fig.  1 I rays  a 
and  b would  be  interpolated,  but  ray  c would  not  be  interpolated  with  either  ray  a or 
ray  b,  even  though  ray  c has  the  same  vertex  depth  as  ray  a and  the  next  successive 
vertex  depth  compared  to  ray  b.  The  key  to  the  interpolation  lies  in  the  path  type  k 
and  cycle  number  n subscripts.  See  Fig.  10  and  Fq.  29  to  obtain  relationships  in 
Fig.  I I . 


There  is  one  very  important  exception  to  the  interpolation  rule  of  only  interpo- 
lating between  successive  rays  of  the  same  class:  the  FORA  program  does  not  calcu- 
late rays  which  vertex  exactly  at  the  source  or  target  depths;  the  rays  vertexing  near 
these  depths  have  vertexing  sound  speeds  slightly  larger  than  the  sound  speeds  at  the 
source  and  target  depths.  If  the  source  or  target  depth  falls  in  a region  of  the  sound- 
speed  profile  where  the  sound  speed  changes  very'  slowly  with  depth,  a ray  vertexing 
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Figure  1 1.  Example  of  ray  interpolation. 

near  the  source  or  target  depth  may  travel  for  several  kiloyards  before  it  recrosses  that 
depth.  The  distance  between  crossings,  as  depicted  by  the  ‘‘hole’’  in  Fig.  12.  is  omitted 
from  interpolation  if  the  rule  that  “only  rays  of  the  same  class  may  be  interpolated”  is 
followed.  This  is  an  unacceptable  limitation  in  the  interpolation  scheme,  since  chang- 
ing the  precision  of  the  calculations  can  change  the  size  of  the  “hole.”  For  a ray  that 
vertexes  very  near  the  source  and  target  depths,  interpolation  is  accomplished 
between  the  ray  path  of  the  first  crossing  and  the  same  path  at  the  second  crossing, 
even  though  the  ray  paths  constitute  two  different  classes.  This  interpolation  scheme 
is  reasonable,  since  A.  R.  T,  D,  and  /i  are  continuous  in  the  affected  region. 

Another  special  case  arises  when  source  and  target  are  at  the  same  depth.  The 
first  ray,  even  though  it  has  a very  small  starting  angle,  often  produces  a range  value 
of  more  than  1 kyd  when  it  crosses  the  target  depth.  This  case  is  dealt  with  easily  by 
setting  Ra  to  0.001  kyd  and  I1A  to  0 dB.  the  reference  values,  and  using  Hq.  3 I for 
interpolation.  A linear  interpolator  for  propagation  loss  versus  range  would  be  inac- 
curate. but  Fq.  31  produces  accurate  spherical  spreading  values. 


I ipurc  12.  t he  “hole"  in  the  interpolation  scheme. 
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Equation  3 1 is  valid  for  interpolation  for  short  ranges,  but  what  about  long 
ranges?  Interpolation  between  successive  rays  spaced  a few  kiloyards  apart  at  ranges 
exceeding  10  kyd  is  so  close  to  linear  interpolation  that  there  is  no  noticeable 
difference. 


Beam  Pattern  Correction  to  Propagation  Loss 

A beam  pattern  is  the  relative  response  in  decibels  of  a transmit  or  receive  array 
of  hydrophones  as  a function  of  some  angle.  Since  ray  paths  in  this  model  are  all 
contained  in  the  vertical  plane,  the  response  need  only  be  known  as  a function  of 
vertical  angle  in  the  vertical  plane  through  own-ship  and  target  positions.  The  angle 
of  maximum  response  defines  the  main  axis  of  a beam  pattern.  The  angle  that  this 
axis  subtends  with  the  horizontal  is  called  the  depression-elevation  (D/E)  angle.  The 
maximum  response  is  0 dB.  with  all  other  responses  measured  relative  to  it.  The 
abscissa!  angular  values  used  in  defining  the  response  function  are  positive  when 
measured  downward  (clockwise)  with  respect  to  the  main  axis.  Since  the  maximum 
of  the  response  function  is  0 dB.  responses  at  all  other  angles  must  be  less  than  0 dB. 

Beam  patterns  as  inputs  to  the  computer  program  are  represented  as  piecewise 
linear  functions  of  angle,  and  the  response  values  are  in  decibels  down  from  the 
main  axis  response.  Thus,  all  response  values  are  entered  as  positive  numbers,  and 
the  program  makes  the  change  of  sign  at  the  time  of  interpolation.  Intermediate 
responses  as  a function  of  angle  are  obtained  by  linear  interpolation.  The  remainder 
of  the  discussion  here  will  treat  responses  as  negative,  according  to  the  usual 
definition. 

The  computer  program  uses  the  theoretical  response  of  a long  line  as  a default 
function  for  beam  pattern.  Both  the  transmit  and  receive  beam  patterns  can  be  repre- 
sented by  this  function,  and  each  can  have  its  own  respective  beamwidth.  The 
response  function  for  a continuous  line  of  length  L is1^ 


/i(0)  = 20  log 


sin  ( sin  l0-0n  1 1 
~j~sin  |0-0nl 


where 

X is  the  wavelength  of  the  acoustic  energy 
0 is  the  look-angle  for  calculating  response 


(32) 


0O  is  the  depression  angle  of  the  main  axis  and  is  obtained  by  mechanically 
rotating  the  line  array,  not  electronically  steering  it 
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Equation  32  can  be  represented  in  terms  of  the  absolute  value  of  the  main-lobe  half- 
power beamwidth  2a  as 


B(a)  = 20  log  — = 20  log 

V? 

The  above  equation  is  in  the  form 

1 _ sin  x (33) 

vT  -v 


.7 TL  . , 

sin  ( — sin  a) 


77 /. 


sin  a 


where 


v 


7 TL  . 

-r—  sin  a 


(34) 


Solving  Eq.  33  for.v  numerically  gives  .v  = 1.39155738.  If  this  value  for.v  is  used  in 
Eq.  34,  the  beam  pattern  response  function  of  Eq.  32  can  be  written  as 


But>)  = 20  log 


sin 


1 39 

(-F—  sin  I0-&.I 
sin  a ° 


1,39 
sin  a 


sin  10-0.,  | 


where 

2a  is  the  beamwidth  at  3.0103  dB  shown  on  each  side  of  the  main  lobe 

The  active  sonar  equation.  Eq.  5.  includes  the  term  IF , which  is  two-way 
propagation  loss  from  the  source  to  the  target  and  back  to  the  receiver  plus  signal 
degradation  caused  by  the  transmit  and  receive  beam  patterns.  It  is  convenient  at  the 
time  of  interpolating  for  propagation  loss  to  use  F (given  by  Eq.  31)  to  interpolate 
linearly  for  the  angle  at  the  source-receiver.  The  beam  pattern  response  is  obtained 
by  interpolating  the  input  transmit  and  receive  beam  patterns  (or  by  calculating  the 
sini  v )/.v  response  function)  using  the  angle  at  the  source.  These  beam  pattern  values 
are  added  to  two-way  propagation  loss  and  the  result  is  called  //,.  The  total  trans- 
mit and  receive  beam  pattern  degradation  is  //|(.  which  is  used  in  Table  1 . 

Except  for  surface-duct  propagation  loss  outputs,  the  LORA  program  never 
displays  true  two-way  propagation  loss  values:  all  displayed  values  contain  beam  pat- 
tern corrections.  If  the  user  desires  to  see  uncorrected  propagation  loss  outputs,  he 
should  input  an  omnidirectional  beam  pattern,  that  is,  zero  response  for  all  angles. 
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MASKING  BACKGROUND 


' 


I- 


The  total  masking  background  level.  A/,  may  be  decomposed  into  two  basic- 
components,  total  noise  level,  Ly,  and  total  reverberation  level,  LR.  Both  noise  level 
and  reverberation  level  are  expressed  in  decibels  (re  I /rbar)  and.  when  added  inco- 
herently, yield  masking  level: 


10 

M = 1 0 log  ( 10  N 


/.„/  1U 

10  R ) 


The  quantities  Ly  and  LR  are  composed  of  primary  sources  of  noise  and  reverberation. 


Noise 


The  total  noise  level  is  assumed  to  include  both  self-noise  and  ambient  noise. 
Self-noise,  L.  is  primarily  composed  of  electrical  receiver  noise,  other  electrical  inter- 
ference. and  acoustic  noise  generated  by  own-ship  motion  and  machinery.  Ambient 
noise  is  produced  by  rain  and  wind  agitating  the  ocean  surface,  distant  shipping,  and 
biological  activity.  Data  for  these  separate  noise  components  are  difficult  to  obtain; 
measurements  of  the  equivalent  isotropic  noise  spectrum  level.  A'.,  are  more  readily 
available.  20-21-22 

Assuming  the  noise  to  be  isotropic  and  noise  spectrum  level  to  be  constant  over 
the  receive  band,  the  equivalent  total  plane-wave  noise  level  in  the  active  band  is 


Ls  = \s  - Dl  + 10  log  t Bwk  MdBi 


(35) 


where 

A’v  is  the  noise  spectrum  level  < d B re  I /ahar  in  a 1-Hz  bandwidth)  for  the  active 
receiving  frequency  at  the  receive  array  (with  self-noise  added  incoherentlv  ) 

1)1  is  the  directivity  index  of  the  receive  array 

and 

liWR  is  the  receive  bandwidth  (II/) 

The  user  may  enter  his  own  value  for  directivity  index.  Otherwise,  he  has  the 
option  of  using  an  approximate  expression  obtained  using  horizontal  and  vertical 
half-power  beamwidths.  AO  and  A0 

0(0.0)  = 1.  if — r A0<  0<  i A0.  and  - Jp  A0  < 0 < -V  A0 

= 0,  otherwise 


3b 


t 


witli 


Dl  = 10  log 
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= 10  log 


n* 

B> 

r_jzL__ 

[A0  sin  (^)J 


b(0,<p ) sin  0 d<j)  dQ 


(361 


where 

Dl  is  directivity  index,  given  in  Uriek* 
h(0,<p)  is  the  beam  pattern  (intensity  ratio) 

6 is  the  azimuthal  angle 
0 is  the  vertical  angle 

A 0 is  the  horizontal  half-power  beamwidth  in  radians 
and 

A 0 is  the  vertical  half-power  beamwidth  in  radians 

Notice  that  for  omnidirectional  transmission  in  the  vertical  plane,  that  is.  A 0 = tt. 

Eq.  36  gives 

"^7T 

D!  = 10  log  (^  ) (vertical  omnidirectionality) 

which  reduces  to  zero  if  the  horizontal  beamwidth  is  2ir.  With  horizontal  beamwidth 
and  vertical  beamwidth  in  degrees,  the  directivity  index  of  Eq.  36  is 


r\  t | p,  | 360  I 

8 Lao  sin  (0 .0087266463  A0)  J 


The  directivity  index  is  only  used  to  calculate  the  bandlevel  noise,  Ls . at  the  out- 
put of  the  beamformer.  Often,  in  system  studies,  Ls  is  known  and  need  not  be  calcu- 
lated by  using  Eq.  35,  but  instead  can  be  entered  directly  as  a program  input. 


'Ref.  19,  p.  48. 
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Reverberation 

Acoustic  energy  generated  by  the  active  sonar  transmitter  propagates  through 
the  ocean  medium,  reflects  and  reradiates  from  a target,  and  propagates  back  to  the 
sonar  receiver.  This  return  energy  is  the  signal.  Some  of  the  energy  put  into  the  water 
by  the  active  sonar  scatters  from  the  ocean  boundaries,  the  surface  and  the  bottom, 
and  from  volume  seatterers  such  as  marine  organisms,  bubbles,  and  suspended  partic- 
ulate matter.  The  return  energy  from  seatterers  other  than  the  target  is  called  reverb- 
eration. Reverberation  masks  the  signal,  just  as  does  noise,  but  differs  from  noise  in 
its  spectral  characteristics.  Noise  has  many  frequency  components  over  a broadband. 
Reverberation  is  narrowband,  with  the  same  frequency  as  the  transmit  frequency 
except,  perhaps,  for  a small  frequency  shift  (doppler)  due  to  transmitter  motion 
relative  to  the  environment.  This  model  assumes  that  the  transmitter  is  stationary 
relative  to  the  environment  and  that  the  doppler  correction  is  made  to  the  signal 
because  of  target  motion  along  the  range  track  between  the  source  and  the  target. 


There  are  three  components  of  reverberation  level  tall  levels  given  in  decibels! 
surface  reverberation  level,  Ls,  bottom  reverberation  level.  La.  and  volume  reverbera- 
tion level,  . The  total  reverberation  level.  J.R , is  expressed  in  terms  of  these  compo- 
nents as 


lR 


1 0 log  ( 1 0 


/ ..  m 


/.„/ to 
It)  H + 
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For  each  scattering  mechanism,  the  LORA  model  uses  only  hackscattering;  that  is, 
the  propagation  path  to  the  scattering  center  is  identical  to  the  path  from  the  scatter- 
ing center  back  to  the  receiver. 

The  bandwidth  of  energy  returned  from  seatterers  may  be  broadened  due  to  the 
dispersive  effects  of  the  scattering  mechanisms.  Also,  the  surface  (wave  motion)  may 
cause  a slight  frequency  shift,  which  is  of  interest  in  high-resolution  sonars.  These 
effects,  however,  are  not  considered  here. 

Both  surface  and  bottom  reverberation  are  strongly  dependent  upon  grazing 
angle  and  boundary  roughness.  Volume  reverberation  is  mainly  characterized  by  the 
biological  organisms  composing  the  scattering  layer.  I he  scattering  layer  extends  in 
depth  to  about  2500  feet  and  is  subject  to  large  diurnal  fluctuations. 


Regardless  of  the  scattering  mechanism  encountered,  unwanted  energy  returned 
to  the  source  is  strongly  affected  by  the  acoustic  propagation  paths.  Both  surface 
reverberation  and  volume  reverberation  can  include  direct,  surface-duct,  bottom- 
bounce.  and  convergence-zone  paths.  Bottom  reverberation  can  include  direct  and 
surface-bounce  paths. 

The  scattering  strengths  attributed  to  the  various  layers  of  seatterers  may  be 
represented  approximately  by  a single  number,  the  scattering  strength  integrated  over 
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the  water  column,  that  is.  the  column  scattering  strength.  Such  an  approximate 
representation  allows  the  volume  element  of  scatterers  to  be  replaced  by  an  effective 
scattering  “area.”  This  concept  is  consistent  with  most  active  sonar  systems  operating 
in  deep  waters,  since  most  currently  operational  sonars  employ  pulse  lengths  exceed- 
ing 100  p's.  Pulse  lengths  of  this  order  correspond  to  resolution  lengths  on  the  order 
of  2500  feet.  Consequently,  ensonifieation  of  the  entire  scattering  column  is  usually 
ensured. 

A general  expression  for  reverberation  level  from  one  scattering  source  may  be 
written  in  the  form 

Lf.  = Lt)  + /iv(0v  ) + BK(<px)  ~-Hk  + Sk(<t>k)  + 10  log  (Ak/A())  in  dB  (38) 

where 

is  the  source  level 

By  is  the  transmit  beam  pattern  response 

Br  is  the  receive  beam  pattern  response 

0Y  is  the  ray  angle  at  the  transmitter  (and  also  at  the  receiver) 

Hk  is  the  one-way  propagation  loss  for  a ray  propagating  between  the  transmitter 
and  A th  scattering  element 

Sy  is  the  backscattering  strength  of  the  Ath  scattering  element  and  is  10  log  of 
the  ratio  of  intensity  in  the  direction  of  backscatter  at  1 yard  from  the 
center  of  a scattering  element  of  area  A()  to  the  intensity  of  a plane  wave 
incident  to  the  element  with  angle  <pk 

/l/r  is  the  effective  scattering  area  of  the  Ath  scattering  surface  or  the  cross- 
sectional  areas  associated  with  column  volume  scattering 

and 

/1(|  is  a reference  area  which  has  the  same  spatial  units  as  ;hc  element  producing 
Sk : A()  is  I yd’  in  the  LORA  model 

Watson  and  McCirr1  describe  a technique  for  approximating  the  effective  total 
scattering  areas  for  the  surface  and  bottom  backseatterers  as  well  as  the  equivalent 
scattering  area  for  the  backscattering  water  column.  Tabic  5 summarizes  these 
approximations.  If  the  units  of  the  scattering  areas  in  Table  5 arc  adjusted  to  square 
yards,  then  the  scattering  areas  are  numerically  equal  to  .1A  . !, , ot  I q.  38.  I he  param- 
eters in  the  expressions  for  scattering  area  are:  Cs,  Cfl.and  ( ',  .sound  speeds  at  the 
ocean  surface,  ocean  bottom,  and  center  of  the  volume  scattering  column.  c\.  0lr 
and  0,  . angles  down  from  the  horizontal  for  rays  incident  at  the  surface,  bottom  and 
volume  scatterers.  r.  pulse  length;  R.  horizontal  range.  Ad.  horizontal  beam  width. 
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Table  3.  Approximations  for  effective  scattering  areas  for 
surface,  bottom,  and  volume  backscatters. 


Reverberation 

Scattering 

Scattering 

Source 

Strength 

Area 

Surface 

(M2CsT)MA0;cos(Xi 

Bottom 

SB 

OI2CbT)RA9/cos  0/; 

Volume 

Sy 

(l/2C/.mA0/cos0, 

Surface  Backscattering  Strength.  For  surface  backscattering  strength.  Ss . the 
LORA  model  uses  the  Chapman-Harris  expressions  for  diffuse  backscattering.'3 
merged  with  Eckarf s high-angle  expression  for  facet  reflection  from  the  ocean  sur- 
face,24,25 and  interpolated  values  representing  Richter’s  data26  for  low  grazing 
angles.  The  empirical  Chapman-Harris  expressions  apply  for  grazing  angles  from  10- 
20  degrees  to  40-60  degrees,  depending  on  wind  speed.  The  Chapman-Harris  expres- 
sions for  backscattering  strength,  sCH . are 


SCH  = 3. 3/3  log (<V30>- 42.1  log  13  + 2.6  ) 
& = 41.6  (^/^r0-58  1 
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where 

<i>s  is  the  grazing  angle  in  degrees  at  the  surface 

Vw  is  the  wind  speed  in  knots 

and 

/is  the  frequency  in  kHz 

Eckart24  derived  theoretical  equations  for  surface  backscattering  strength  for 
sound  impinging  on  the  ocean  surface  at  angles  steep  enough  to  preclude  sell- 
shadowing at  the  surface.  His  general  expression  requires  knowledge  of  the  spatial 
spectrum  of  the  surface  waves:  the  ocean  surface  must  be  described  by  two  compo- 
nents of  the  surface  wave-number  vector,  depending  on  the  direction  and  wavelength 
of  the  sinusoidal  components  of  the  surface  waves.  W ith  the  simplification  that  the 
wavelength  of  the  incident  sound  is  much  shorter  than  the  ocean  surface  wavelengths 
reflecting  the  sound,  the  following  expressions  describe  surface  backscattering  strength  . 

Sr  = -1 0 log  (8jt(*2 ) - 2. 1 7 a"2  cot2  <t>s  (40) 


where 


a2  = 0.003  + 0.00294  Vw  (mean-square  waveslope) 


(41 1 
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is  the  grazing  angle  (off  the  horizontal)  at  the  surface 


and 

Vw  is  the  wind  speed  in  knots 

Chapman  and  Scott’5  reported  data  for  high  angles  of  incidence  and  recommended 
the  use  of  Eckart’s  theory,  which  they  expressed  in  the  form  of  Eq.  40.  They  also 
recommended  the  Cox  and  Munk27  empirical  expression  for  random-direction 
mean-square  waveslope  given  by  Eq.  41,  modified  for  units.  Incidentally,  Cox  and 
Munk  also  give  with-the-wind  (a^, ) and  crosswind  (a". ) mean-square  waveslopes  of 

a2  = 0.00316  W 
w 

and 

a2  = 0.003  + 0.00102  W 

c 

where 

W is  wind  speed  in  m/s 

The  LORA  model  does  not  consider  wind  direction,  but  may  in  the  future  if  a satis- 
factory low-grazing-angle  theory  supported  by  data  over  the  full  active  sonar  fre- 
quency range  becomes  available. 

Gardner28  finds  significant  backscattering  at  low  grazing  angles  at  the  ocean  sur- 
face for  5-kllz  sound.  His  reported  data  and  theoretical  derivation,  which  uses  diffuse 
backscattering  from  a self-shadowing  surface,  indicate  that  the  surface  backscattering 
strengths  should  be  20  to  40  dB  higher,  for  grazing  angles  less  than  !0  degrees,  than 
those  predicted  by  the  Chapman-Harris  expressions.  Eq.  30.  Almost  every  investigator 
in  the  field  who  noticed  the  high  backscattering  strengths  at  low  grazing  angles 
attributed  the  effect  to  near-surface  volume  scatterers  such  as  biological  organisms 
; or  resonating  air  bubbles.  Gardner’s  work  shows  the  theoretical  possibility  of  high 

' backscattering  strengths  at  low  grazing  angles.  Unfortunately  he  validates  his  predic- 

tions at  only  one  frequency,  5 kHz.  His  reported  backscattering  strengths  are  I 2 dB 
higher  than  Richter’s  data-6 at  the  same  frequency.  Richter's  data  spans  the  frequen- 
cies 0.8-1  2.8  kHz  and  the  grazing  angles  2-38  degrees.  Richter  interpreted  his  high 
backscattering  strengths  at  low  grazing  angles  to  be  caused  by  near-surface  volume 
scatterers.  The  trends  in  his  data  at  the  various  frequencies  would  seem  to  support 
his  view.  Since  Richter's  data  spans  the  frequencies,  grazing  angles,  and  wind  speeds 
of  interest,  and  since  Richter’s  data  indicates  substantially  higher  values  than  those 
predicted  by  Eq.  3‘>  for  low  grazing  angles,  Richter’s  data  will  be  used  for  a lower 
limit  on  surface  backscattering  strength  for  grazing  angles  of  0-20  degrees.  This 
choice  is  a compromise  in  Gardner's  direction,  for  he  says  the  low-angle  backscatter- 
ing strengths  are  even  higher.  Since  Richter  did  not  summarize  his  data  as  some 
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function  of  wind  speed,  frequency,  and  grazing  angle,  a reduced  data  set  based  on  his 
reported  data  set  will  be  interpolated.  The  reduced  data  set  is  given  in  I able  4.  In  the 
LORA  computer  program  the  data  is  linearly  interpolated,  with  the  appropriate  end 
value  used  as  the  default  value  if  any  of  the  parameters  of  wind  speed,  frequency,  or 
grazing  angle  is  out  ol  range. 


Table  4.  Reduced  Richter  data  tor  surface 
backscattering  strength. 


Wind  Speed 
5 knots 

Wind  Speed 
20  knots 

Geometric 

frequency 

(kHz) 

Grazing  Angles  deg 

Grazing  Angle,  deg 

6 

15  38 

6 

15  38 

1.15 

-57 

-57  -46 

-57 

-48  -39 

2.3 

-58 

-63  -45 

-58 

-50  -39 

4.5 

-52 

-57  -51 

-52 

-45  -42 

9.0 

-48 

-55  -5 1 

-48 

-46  -30 

The  model  for  surface  backscattering  strength  is  a combination  of  the  Richter 
data  the  Chapman-Harris  expressions,  and  the  Chapman-Scott  expressions.  Surface 
backscattering  strength  is  calculated  by  each  method  for  the  grazing  angle  ol  interest 
and  the  maximum  value  is  retained.  The  results  for  some  selected  frequencies  and 
wind  speeds  are  .mown  in  Fig.  13. 


Bottom  Scattering  Strength.  The  expression  for  backscattering  from  the  bottom 
includes  a term  for  diffuse  scattering  at  small  angles,  in  accordance  with  Lambert  s 
law.  This  term  was  derived  by  Mackenzie-9  and  merged  by  Watson  with  a large- 
angle  facet  reflection  term  based  on  Schmidt's  data  from  the  Norwegian  Sea.  u 
resulting  expression  for  scattering  strength  in  decibels  for  I yd-  ol  bottom  measured 
at  I yard  in  the  direction  of  backscattering  is 


SH  = 1 0 log  \nB  sin2  <t>B  + exp  (-1 00  cot2  4>B)]  < 4 - > 

where 

fj"  = )(/  « and  is  the  bottom  backscattering  coefficient 
U is  the  bottom  backscattering  coefficient  in  dB 

n 

and 

<$>}j  is  the  grazing  angle  at  the  bottom 

Figure  14  shows  Lambert’s  law  scattering  for  VR  = -27  dB.  facet  reflection  lor 
the  Norwegian  Sea,  and  the  merging  of  the  two  curves  resulting  in  Lq.  4_. 
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Volume  Backscattering  Strength.  The  primary  mechanism  responsible  tor  vol- 
ume reverberation  consists  of  migrating  biological  scatterers  extending  from  the 
surface  to  depths  of  about  2500  feet.  Often  these  scatterers  occur  in  identifiable 
layers  which  exhibit  a diurnal  vertical  migration.  The  scattering  strengths  ot  these 
layers  depend  on  geographic  location  and  season  as  well  as  time  ot  day.  Groups  ot 
scatterers  have  been  known  to  collect  and  then  disperse  within  a tew  minutes,  giving 
rise  to  transient  scattering  strengths  of  significant  magnitude.32  Consequently,  meas- 
urements of  volume  scattering  strength  as  a function  of  depth,  while  providing  a 
means  of  studying  the  behavior  of  migrating  scatterers,  do  not  yield  values  considered 
stable  enough  for  most  performance  prediction  efforts. 

The  LORA  model  uses  the  concept  of  backscattering  from  a water  column  con- 
taining uniformly  distributed  volume  scatterers.  I he  water  column  is  defined  as  a 
column  with  1 yd2  cross-sectional  area  extending  from  the  ocean  surface  to  the  bot- 
tom of  the  deepest  scattering  layer.  If  the  backscattering  coefficient.  p,  tor  1yd  is 
known  as  a function  of  depth,  then  the  column  backscattering  strength.  .S’,  can  be 
obtained  by  integrating 


where 

Zc  is  the  depth  of  the  deepest  scattering  layer  in  feet 

The  user  must  supply  the  value  for 5,  for  input  to  the  LORA  computer  program. 
He  must  find  data  either  for  column  strength  or  for  depth-dependent  volume  scatter- 
ing strength  and  perform  his  own  numerical  integration.  C omprehensive  discussions 
of  column  strength  can  be  found  in  Ret.  33  and  34. 

Reverberation  Ray  Paths.  Rays  are  traced  from  the  source  depth  to  bottom, 
surface,  and  volume  backscatterers.  Many  different  classes  ol  ray  paths  are  defined. 
Ray  paths  within  the  same  class  can  be  interpolated;  rays  in  different  classes  cannot 
be  interpolated. 

Figure  15  shows  17  classes  of  ray  paths  which  may  contribute  to  the  reverbera- 
tion level  at  the  time  the  echo  returns  to  the  receiver  from  a target  in  Zone  1 (this  is 
the  approximate  region  where  there  exists  some  ray  path  to  the  target  which  b s not 
completed  a cycle).  All  of  the  paths  in  Fig.  I 5 are  specular.  Also,  the  direction  of 
backscattering  is  along  the  path  of  the  incident  ray.  Paths  4 and  in  I ig.  1 5 appear 
to  be  redundant  but  are  not.  Path  4.  if  continued,  would  strike  the  bottom,  like 
Path  7.  Path  5.  if  continued,  would  not  strike  the  bottom,  like  Path  l>.  The  surtace- 
bottom-surface  path  which  terminates  at  a bottom  backseat terer  is  not  included 
because  it  contributes  less  to  reverberation  level  than  Path  3.  which  itself  usually 
makes  a negligible  contribution  to  the  reverberation  level  in  a deep  ocean. 
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Long-range  reverberation  paths  are  shown  in  Fig.  16.  One  to  four  cycles  are 
added  to  Paths  8,  9,  14.  1 5,  16,  and  1 7 ot  Fig.  1 5 in  order  to  obtain  multizone 
reverberation  (Zones  2 through  5).  Only  those  paths  which  do  not  strike  the  bottom 
are  used  tor  multizone  reverberation,  since  the  high  bottom  loss  associated  with 
active  sonar  frequencies  reduces  multiple-bottom-bounce  reverberation  levels  to 
insignificance.  Shallow-water  studies,  however,  warrant  using  contributions  to  reverb- 
eration from  multi-bottom-bounce  paths.  The  appropriate  paths  can  be  incorporated 
into  the  LORA  program  if  needed. 

Reverberation  level  is  obtained  by  summing  the  contributions  from  all  backscat- 
terers.  There  are  41  different  classes  of  reverberation  paths,  1 7 shown  in  Fie.  15  and 
24  indicated  in  Fig.  16.  A maximum  of  40  rays  are  traced  for  each  path,  and  corre- 
sponding two-way  travel  times  are  calculated.  Propagation  loss  is  computed  for  each 
ray.  and  transmit  and  receive  beam  pattern  corrections  are  added  to  the  loss.  The 
appropriate  backseat tering  area  in  Table  3 is  calculated,  as  well  as  one  of  the  appro- 
priate scattering  strengths:  bottom,  surface,  or  volume.  These  values  and  the  source 
level  are  used  in  F.q.  38  to  obtain  reverberation  in  dB  re  1 pB.  which  theoretically  is 
measured  at  the  output  of  the  beamformer.  These  reverberation  levels  and  their 
corresponding  times  are  stored  in  two-dimensional  arrays  of  41  ray  classes,  with  40 
rays  in  each  class.  Minimum  and  maximum  times  of  each  group  of  40  rays  are  stored 
lor  later  examination  to  determine  if  interpolation  is  appropriate  for  a given  two-way 
travel  time  trom  source  to  target  as  a function  of  range.  If  the  minimum  and  maximum 
times  bracket  the  target  echo  time,  the  array  of  reverberation  values  is  interpolated. 

It  a reverberation  array  is  multivalued  in  time,  the  array  is  interpolated  for  each 
instance  in  which  reverberation  time  crosses  the  target  echo  time.  All  reverberation 
components  contributing  power  to  the  receiver  at  the  time  of  the  target  echo  return 
are  added  incoherently  (the  phases  ot  the  levels  are  random),  so  that 

Lf.  = 10  log  ( I'  lO^'710)  (43) 

where 


and 

lv  are  reverberation  levels  which  belong  to  the  various  paths  shown  in  Fig.  15 
and  16.  grouped  according  .o  bottom,  surface,  or  volume  backscattering 
mechanism 

I he  Lf,  are  summed  by  Fq.  37  to  give  total  reverberation. 

Surface  and  Subsurface  Duct  Reverberation.  Because  predicted  propagation 
losses  are  suspect  when  ray  tracing  is  used  for  relatively  shallow  ducts,  the  duct  re- 
verberation calculations  use  propagation  losses  predicted  by  the  empirical  AMOS 
equations.  It  a surface  duct  is  present  and  the  source  depth  is  both  'ess  than  1000 
teet  and  less  than  32^/Z,  teet.  where  Zj  is  the  layer  depth  in  feet,  then  the  AMOS 
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I ijiuro  16.  Reverberation  paths  lor  nth-zone 
backseat  terers. 
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equations  are  used.  Propagation  loss  from  source  to  surface  is  equivalent  to  propaga- 
tion loss  from  source  to  a target  deptlr  of  zero.  Propagation  loss  from  source  to  vol- 
ume backscatterers  is  obtained  by  setting  target  depth  to  one-half  the  layer  depth. 

The  reverberation  equation,  Hq.  38,  is  applicable  for  reverberation  calculations 
based  on  the  AMOS  equations,  with  certain  assumptions  for  angles  of  incidence.  The 
transmit  and  receive  beam  pattern  responses  needed  in  Hq.  38  are  computed  by  using 

</>v  = -y  <t>! , source  in  the  surface  duct 
= <t>! , source  below  the  surface  duct 


where 


*i. 


cos 


and  where 

C’v  is  the  sound  speed  at  the  source 
and 

CJ  is  the  sound  speed  maximum  in  the  layer 

The  surface  backscattering  strength  is  computed  for  the  angle  at  the  surface 
given  by 


rv 


d-  = COS'1  (rr-COS  0V  ) 
C/. 


where 

Cs  is  the  surface  sound  speed 
and 

C'x  and  <px  are  defined  above 

The  surface  backscattering  area  needed  in  Hq.  38  is  easily  calculated  from 
Table  3.  since  all  the  parameters  are  known. 

The  volume  scattering  strength  for  a uniform  scattering  column  is  functionally 
related  to  the  depth  of  the  scattering  column.  Zc.  by 

5,  = l\  * 10  log  Z(. 
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where 


{/,  is  a constant  containing  scattering  strength  per  unit  volume  and  a constant 
for  conversion  of  units 

The  assumption  that  the  backscatterers  are  uniform  with  depth  gives  the  scattering 
strength  in  the  duct,  with  layer  depth.  Zy  : 

SDy  = Uy  + lOlogZ,  , ZL<ZC 

= Uv  + IO!ogZc.  , Z,  >ZC 

Thus,  the  volume  scattering  strength  of  the  duct  is  related  to  the  scattering  strength 
of  the  column,  Sy,  by 

Spy  = Sy  + 10  log  (Zl  /Zc).  Z,  <Z( 

= Sy,  ZL  >ZC 

Sy  is  input  parameter,  and  Snv  is  automatically  calculated  in  the  program.  The  scat- 
tering area  associated  with  the  backscattering  column  in  the  duct  is  easily  computed 
from  Table  3. 

Subsurface  ducts  provide  only  volume  reverberation  and  only  if  the  source  is 
located  within  the  boundaries  of  the  duct.  The  columnar  volume  backscattering 
strength  for  the  subsurface  duct  is 

S.soi’  ~ Sy  + 1 0 log  <Z(),/Z(.) 


where 

Zu,  is  the  width  of  the  duct  or  the  portion  that  overlaps  the  scattering  column 
whose  maximum  depth  is  Z( 

Multiple-Ping  Reverberation.  Reverberation  level  is  a function  of  time  and  is  the 
incoherent  sum  of  squared  pressure  levels  produced  by  all  of  the  backscattering 
mechanisms  described  above.  It  a sonar  system  broadcasts  more  than  one  ping,  later 
pings  will  have  their  target  echo  signals  masked  by  the  additional  incoherent  reverb- 
eration levels  of  earlier  pings.  The  LORA  program  allow-,  the  user  to  specify  total 
number  of  pings.  X,,.  and  the  time  between  pings.  Tf,.  To  the  reverberation  level  at 
time  T.  the  two-way  time  from  source  to  target,  are  added  (incoherently,  according 

to  the  form  of  I q.  43 1 the  reverberation  levels  at  times  T + Tp,  T + 2Tp 

T + {X p—  I )Tp.  These  reverberation  levels  are  obtained  by  interpolating  all  members  of 
the  stored  reverberation  table  which  have  starting  and  ending  times  bracketing  the 
extended  time  of  an  earlier  ping.  I'he  user  can  specify  one  to  nine  pings,  with  single- 
ping reverberation  as  the  default  option. 
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Doppler  Gain.  In  many  active  sonar  tactical  situations,  target  doppler  can  he  the 
determining  factor  in  target  detectability.  Doppler  gain  is  the  increase  of  signal-to- 
reverberation  ratio  caused  by  the  shifting  of  the  target  echo  frequency  band  away 
from  the  reverberation  frequency  band.  Doppler  gain  occurs  when  the  target  has  a 
velocity  component  relative  to  the  ocean  along  the  range  track  between  the  target 
and  the  transmitter-receiver.  A doppler  shift  in  frequency,  Af,  occurs  between  the 
target  echo  and  reverberation  center  frequencies.  If  the  receiver  is  tuned  to  the  target 
return  frequency. /+  Af,  then  the  reverberation  level  is  reduced  in  the  doppler 
receive  band  as  compared  to  the  reverberation  level  in  the  receive  band  if  no  shift 
had  occurred.  It  is  difficult  to  obtain  an  accurate  estimate  of  doppler  reduction  of 
reverberation,  since  pulse  shape  (in  the  frequency  domain),  target  and  reverberation 
echo  distributions,  and  sonar  operator  skill  in  finding  the  proper  receive  frequency 
and  bandwidth  are  not  model  inputs.  However,  some  assumptions  are  possible  which 
can  give  a rough  but  tidy  model  for  estimating  doppler  gain. 

The  first  assumption  is  that  if  the  source  emits  a tone  pulse  at  frequency  f,  then 
the  return  intensities  from  backscatterers  or  target  are  gaussian  distributed.  The  sec- 
ond assumption  is  that  all  distributions  have  the  same  standard  deviation  equal  to 
one-half  the  receive  frequency  bandwidth,  BWK.  Thus,  it  is  assumed  that  the  sonar 
operator  can  match  the  bandwidth  of  the  target  echo.  The  third  assumption  is  that 
the  sonar  operator  can  set  the  receive  frequency  precisely  at  /+  Af.  the  doppler- 
shifted  target-echo  frequency.  The  fourth  assumption  is  that  all  significant  energy  to 
be  used  for  target  detection  falls  in  the  new  band./+  A ft  1 / 2 BWR 

Figure  1 7 shows  the  assumed  standardized  normal  distributions  for  the  target- 
echo  intensity  and  reverberation  intensity.  Since  all  energy  for  detection  is  in  the 
bandwidth  Rwr  . only  the  shaded  area  under  the  reverberation  distribution  curve  with- 
in that  band  will  mask  the  signal.  It  must  be  noted  that  since  each  of  the  distributions 
is  normalized,  the  actual  intensity  distributions  are  obtained  by  a multiplicative  con- 
stant or  an  additive  level  in  decibels.  The  difference  in  area  of  the  two  curves  in  Fig. 

17  is  used  to  determine  a bias  in  decibels  which  is  subtracted  from  the  reverberation 
level,  either  from  the  separate  reverberation  components  of  Eq.  37  or  from  the  total 
reverberation  level. 

The  standardized  normal  distribution  is 


<I>(  v) 


V 


f-f 

_7T  *z  _c 


,'\p  ( — — ) dy 


where 

v is  in  units  of  one  standard  deviation 
rite  area  under  the  target  return  curve  within  the  receive  band.  BWR  . is 
Ar  = <l>(  I ) - «!>(-! ) = 0.68 
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area  to  area  under  heavy  line!. 


The  area  under  the  normalized  reverberation  intensity  distribution  within  the  band 
BWR  centered  on /+  A/ is 


/I*  = <l>|(  A />  0.5  Bwr  )/ 0.5  Bwr  | - lJ>((  A/  - 0.5  BWR  )/0.5  BWR  i 
= T( 2Af/BWR  + !)-«!><  2A//fl  - 1) 


Tlie  doppler  gain  is  Gn 
(Jn  - 1 0 log 


10  log  {At!Ar  ) or 


0 OX 

‘H  2Aj/BWR  + D-4-CA f!BWR 


where 

A/  = 0.7Tr/Uz12 

Vc  <s  target  velocity  component  (along  the  range  track  between  the  source  and 
the  target  and  relative  to  the  ocean  as  the  frame  of  reference) 


/ is  the  transmitter  frequency  in  kHz 


Iht  total  reverberation  level,  LR,  from  stationary  boundary  and  volume  back- 
scatterers  is  corrected  tor  the  doppler  shift  away  front  the  receive  band  by 

Lr  ('t> 


where 

L’h  is  the  corrected  reverberation  level 


SIGNAL  PROCESSING  MODELS  FOR  PROBABILITY  OF  DETECTION 

1 he  LORA  model  contains  five  optional  models  for  predicting  the  probability 
ot  detection  ot  a signal  masked  by  white  gaussian  noise  and  reverberation: 

1 . Linear  correlator  and  a known  receive  signal. 

-■  Lmcar  correlator  and  a constant-amplitude,  phase-distorted  receive  signal. 

3.  Quadrature  correlator  and  a constant-amplitude,  phase-distorted  receive 
signal. 

4.  Quadrature  correlator  and  a random-amplitude,  random-phase-distorted 
receive  signal  (Rayleigh  channel l. 

5.  Threshold  detector  tor  total  energy  in  the  receive  band  and  a random-phase 
gaussian-amplitude  receive  signal. 

I he  first  tour  signal  processing  models  are  presented  by  Schumacher  and 
leeter.  L.  K.  Arndt  ot  NUC  made  modifications  to  those  models  to  replace  the 
time-bandwidth  products  with  a functional  dependence  on  recognition  differential 
I he  filth  signal  processing  model  is  presented  by  Watson  and  McCirr.1  All  five  models 
will  be  presented  after  the  following  sections,  which  describe  the  probability  of  false 
alarm  and  the  modification  to  recognition  differential  to  account  for  reverberation- 
limited  masking  noise  and  signal-clipping  losses. 

Probability  of  False  Alarm 

I he  probability  ot  talse  alarm.  Pfa,  is  the  probability  that  some  threshold  is 
exceeded  in  some  resolution  cell  of  a sonar  plan  position  indicator  (or  some  equiva- 
lent display)  so  that  an  operator  is  convinced  that  a target  is  present  when  it  is  not 
Sinee  /’o  applies  only  to  one  resolution  cell,  the  total  probability  of  false  alarm  for 
the  indicator  is  given  by 


PFA  = 


(45  l 


where 


N 's  the  total  number  of  resolution  cells  in  the  dopl.n 


i 


If  the  display  is  divided  into  1 -degree  sectors  and  20  range  bins,  a 360-degree  indi- 
cator would  have  7200  resolution  cells.  Thus,  in  practice  PJa  assumes  a small  magni- 
tude, with  typical  values  ranging  between  I O'2  and  10" 12 . Usually  P,a  tails  in  between 
10  4 and  10-6  (L.  K.  Arndt,  discussion).  The  model' input  is  Pfa  and  not  Py4 . In  the 
sonar  setting,  the  operator  can  change  his  Pja  by  increasing  or  decreasing  the  gain  on 
the  indicator.  Certain  subjective  influences  on  the  operator,  such  as  a requirement 
that  no  false  alarms  be  reported,  can  produce  wide  variance  in  estimates  of  Pfa  How- 
ever, the  probability-of-detection  models  are  relatively  insensitive  to  P(a.  In  all  the 
detection  models  presented  here,  the  probability  of  detection  is  much  more  sensitive 
to  the  recognition  differential. 

Recognition  Differential 

The  Case  of  Reverberation-Limited  Masking  Background.  The  recognition  dif- 
ferential is  the  ratio  in  decibels  of  the  signal  level  to  the  masking  level  required  for  a 
U.5  probability  of  detection.  The  recognition  differential  is  a psychoacoustic  function 
of  the  sonar  operator’s  mental  outlook  and  the  characteristics  ot  the  masking  noise. 

At  the  same  false  alarm  rate,  operators  will  tend  to  detect  a signal  less  often  in  a 
reverberation-limited  background  than  in  an  ambient-noise-limited  background.  The 
performance  of  the  linear  correlator  has  been  investigated  for  the  presence  of  a known 
signal  in  gaussian  white  noise  and  in  reverberation-limited  noise.36-37  The  results  indi- 
cate that  reverberation-limited  noise  increases  the  recognition  differential  2 dB  over 
an  ambient-noise-limited  background  for  both  CW  and  FM  transmit  signals.  Arndt  has 
produced  an  expression  to  give  a smooth  transition  for  the  correction  term  Ar  to  be 
added  to  the  recognition  differential  (see  Fig.  18): 


where 


erfc(.v) 


(47) 


< 


and 

Rxr  is  the  ratio  of  reverberation  to  ambient  noise  in  dB 

The  input  recognition  differential  supplied  by  the  user  of  the  computer  program  is 
the  recognition  differential  for  the  ambient-noise-limited  performance  model. 

Clipping  Correction  for  Recognition  Differential.  For  multielement  receiving 
arrays,  signal  digitization  can  produce  effective  loss,  since  the  discrete  time  sampling 
and  nonlinear  quantization  of  the  signal  level  produce  a nonoptimum  signal  at  the 
replica  correlator  as  compared  to  a signal  received  by  an  analog  linear  receiver.  For 
the  assumption  of  a very  low  signal-to-noise  ratio.  Horton3,  ami  Schultheiss  s have 
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Reverberation-to-Noise  Ratio  (dB) 


I igurc  18.  C orrection  curve  tor  recognition  differen- 
tial, showing  the  transition  from  a noise- 
limited  to  a reverberation-limited 
situation  (L.  K.  Arndt,  discussion ). 


determined  upper  bounds  on  clipping  loss  for  a replica  correlator.  Horton  specifies 
clipping  loss  as  a 2-dB  response  deterioration  or  as  a 1-dB  signal-to-noise  loss 
I xpcrimentul  results  by  Teeter’6  for  pulsed  CW  and  linear  FM  signals  indicate  the 
following  clipping  losses  to  be  subtracted  from  the  signal  level: 

Clipping  Loss,  Ambient-Noise-Limited:  2 dB 

Clipping  Loss.  Reverberation-Limited:  I dB 

Instead  of  the  correction  being  applied  to  the  signal  level,  it  is  here  applied  to  the 
recognition  differential.  Arndt’s  smooth  transition  curve  (see  I ig  I'M  can  be  used  to 
obtain  the  correction  term,  which  is 

R\k 

A = I + erfc  ( — ^ — ) for  clipping  <4N) 

= 0 for  no  clipping 


where 

RSI<  is  the  ratio  of  reverberation  to  noise  in  dB 
and 

erfc(.v)  is  defined  by  Lq.  47 


I 
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Reverberation-to-Noise  Ratio  IdB) 


I igure  19.  Clipping  loss  in  the  presence  of 
reverberation. 


The  Corrected  Recognition  Differential.  The  corrected  recognition  differential. 
R[)c-  is  obtained  by  combining  Lq.  4(>  and  48: 


V rd  + \ + \ 

r\k 

Rnc  ~ rd  + - _ - erfe  ) • 110  clipping 
Rnr 

= /?D  + 3 - erfc  (— — ) , clipping 


<4‘>) 


Probability  of  Detection 

The  following  discussion  will  deal  with  several  models  for  probability  of  detec- 
tion for  various  assumptions  of  signal  degradation. 

Model  1 Performance  of  the  Linear  Correlator  With  a Known  Signal  in  White 
Gaussian  Noise.1'  The  linear  correlator  (Fig.  20)  is  (by  the  maximum-likelihood  and 
minimum-mean-squared-error  criteria)  the  optimum  receiver  structure  for  the  detec- 
tion of  a completely  known  signal  in  the  presence  of  additive  white  Gaussian  noise. 
n(t A39  The  received  signal,  sR  ft),  is  at  worst  an  attenuated,  phase-shifted  version  of 
the  transmitted  signal,  sJt): 

sR(t)  = kA(l)  sinlcu/  + <>(t)  + 0 | , 0 < / < T 


= 0.  otherwise 
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k is  an  arbitrary  constant 

A< i)  is  the  amplitude 

o>  is  the  transmit  frequency  + 2ir 

<t>( t)  is  the  phase 

6 is  an  arbitrary  constant 

and 

T is  the  signal  duration 

The  receiver  is  matched  to  sx  ( t).  This  is  equivalent  to  the  unlikely  situation  (especially 
for  echo  ranging)  where  the  propagation  channel  does  not  distort  the  signal  in  ampli- 
tude or  phase. 

The  performance  of  the  linear  correlator  under  optimum  conditions  is'9 

PD  = erfc  (x-d)  (50) 

where 

Pp  is  the  probability  of  detection 

.v  is  the  threshold  corresponding  to  a given  false  alarm  rate.  Pfa 

and 

d is  the  output  signal-to-noise  ratio  given  by 


/. 

(10 


5 b 


where 


= r*2 

J n K 


(t)  dt  = I T 


is  the  mean-squared  received  signal  power 
T is  the  signal  length 

n is  the  mean-squared  input  noise  power  in  a 1-Hz  bandwidth 
Z.  v is  the  input  signal-to-noise  ratio  in  dB 


Bwu  is  the  receiver  bandwidth  in  Hz 
A new  term,  the  processing  gain.  gp.  is  defined  here  as 


*/•  s - tbwr 


Then  tq.  50  becomes 


'sy/u\ 


PD  = crfc  I x-(gP  10"  l I 

The  probability  of  false  alarm.  Pfa.  is  the  probability  that  the  threshold  v is  exceeded 
when  there  is  no  signal  energy,  so 

Pn  = P,a  = ertc(.v)  (53) 

Thus,  given  Pfa.  the  threshold  can  be  determined  by  inverting  Lq.  53.  In  the  computer 
program,  inversion  is  done  by  interpolating  tabulated  values. 

The  recognition  differential  is  defined  as  the  value  ot  signal-to-noise  ratio  where 
Pn  = 0.5  for  a given  P,a.  Thus,  using  F.q.  5 2 

KJH) 

0.5  = erfc  ( v - gr  10  ' ) 


erfc(0)  = 0.5 


»Ji(i 

- (gr  10  ) 


with  the  result 


g,. 


10 


■V 
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Using  Eq.  5 2 and  54 


PD  = ert'c  [x(l  - 10 


iLSK  ' HDh'2'\ 


(54) 


where 

x is  obtained  by  inverting  Eq.  53. 

This  signal  processing  model  is  included  because  it  provides  an  upper  bound  on 
signal  processor  performance.  The  reverberation-limited  case  is  realized  by  substitut- 
ing/^. for  Rd  . 

Model  2 - Performance  of  the  Linear  Correlator  With  a Constant-Amplitude. 
Phase-Distorted  Signal  in  White  Gaussian  Noise.  The  assumption  is  made  here  that 
the  medium  causes  phase  distortion  ot  the  signal,  but  not  amplitude  distortion,  f or 
the  transmit  signal 


Schumacher  and  Teeter"5  have  shown  that  this  model  is  very  nearly  equivalent  to 
Model  1 , with  the  signal-to-noise  ratio  degraded  by  2 dB.  The  probability  of  detection 
is 

PD  = erfe  |.v(  1 - 1 (f^oc"2 , | 

where 

.v  is  the  threshold  obtained  by  inverting  Eq.  53 

Z.vv  is  the  signal-to-noise  ratio  in  dB 
and 

Rnc  is  the  corrected  recognition  differential  of  Eq.  49 

Model  3 — Performance  of  the  Quadrature  Correlator  With  a Constant-Ampli- 
tude. Phase-Distorted  Signal  in  White  Gaussian  Noise.  The  quadrature  correlator 
matched  filter  is  the  optimum  receiver  (by  the  maximum-likelihood  criterion)  for 
the  detection  of  a constant-amplitude,  uniformly  distributed  random-phase  signal. 
sR( t),  in  gaussian  white  noise,  n( t).  The  transmit  signal  is  described  by  Eq.  55.  and 
the  receive  signal  is  described  by  Eq.  56.  The  quadrature  correlator  is  schematically 
represented  in  Fig.  2 1 . 


sR(t)+mn 


I inure  21.  Quadrature  correlator  matched  filter. 

The  performance  of  the  quadrature  correlator  under  these  conditions  is39 

I 

C 00  2 i1 

Pn  = Q(d.x)  = / ; exp  * c — - lAzd)  dz  (57) 

Jx 

where 

(J  is  known  as  Marcum's  Q Function4'1 
/n  is  a modified  Bessel  function  of  the  first  kind4" 
v : v' -2  In  /' , is  the  threshold 


and 


d = 


and 


Tl%  thc  re^eive  signal  pulse  length 
BWR  1s  the  receive  signal  bandwidth 

LSN  1s  the  signal-to-noise  ratio  in  dB 


time-bandwidth  P rod  nit  h7'/y|,r”  "'{jsl  n ^ d"  57° differentiaI  to  remove  the 


^ - Q\%  lo'^'V'2  , 

By  the  definition  of  recognition  differential 


05  = . v^JhTT 


fa 


(58) 


Of  detection!^  '"°de'  "’PUtS'  S<‘  Can  l,e  foun<*  By  interpolation.  The  probability 


i) 


Ql(gplOiLs^Roc,'l° 


fa 


where 

Q is  defined  by  Eq.  57 

8P  ls  numerically  evaluated  using  Eq.  58 

Ls»  is  the  signal-to-noise  ratio  in  dB 

R°c  ‘S  the  COrrected  recognition  differential  of  Fq  40 
and  1 

P"  ' "K  Probability  of  falsy.  ..|,irtli 

ami  RSmPha^r:,^  m‘  ***''""' 

Ins  is  the  most  reasonable  of  all  the  models  m . a8"?  VVhl,e  Gal,ss,an  Noise. 

— — <■ »*..  T„e 
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the  signal  in  amplitude  and  phase.  The  transmit  signal  is  given  In  Kq.  55  and  the 
receive  signal  by 


{(r)  = k £ a{A(t)  sin  [cor  + <j)(tJ  + 0f  | [u(t-t^  ) - )] 


i+ 1 


where 

a.(i  = 1 .2 are  independent  samples  from  a Rayleigh  distribution 

and 

0-,  uU-tj),  and  t(  are  defined  in  Eq.  56 

The  performance  of  this  model  is  characterized  by  the  probability  of  detection  w 


/>  = /> 
'd  fa 


(iw-) 


(59) 


where 


d2  = 2 TBwr  10  SN 


/.  1 o 


Again  if  Eq.  51  and  the  definition  for  Rp  are  used,  TBWR  can  be  eliminated  for  RD : 

I 


0.5  = P. 


( •V'"1) 

\ 1 +gp  to  ,K  / 


la 


Thus,  the  processing  gain  is 


log 


£/> 


>/  1 p~K  DCIU) 


p log  2 


and 


log 


(44 


d2  - —r'F  io^'^10 

log  2 


for  use  in  Eq.  59. 
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Model  5 - Random-Phase,  Gaussian-Amplitude  Signal  in  White  Gaussian  Noise. 
From  the  Central  Limit  Theorem,  the  distribution  of  the  signal  level.  Z-sv,  tends  to  be 
gaussian  when  it  is  the  sunt  of  independent  variables  having  the  same  distribution  as 
the  others.  The  signal  level  in  the  presence  ot  noise  is 

Lsn  = SL  - 2 //  + BxUf>x  ) + Br(<I>r)-M 

which  is  assumed  to  be  gaussian,  since  the  source  level.  S;  : one-way  propagation  loss. 
//;  transmit,  Bx , and  receive,  BR , response  patterns;  and  masking  noise.  M,  are  random. 
The  signal  excess  is 

^t  :x  ~ Asv  ~ ^dc 


where 

R/)C  is  the  recognition  differential  corrected  tor  the  presence  ol  reverberation 
and  for  signal  clipping,  as  in  Eq.  49 

The  probability  of  detection  is 
PD  = <i> 

where 

a is  the  standard  deviation  of  the  signal  excess  in  decibels 

and 

4>  is  the  standardized  normal  distribution  function  given  by  Eq.  44 

Note  that  the  probability  of  false  alarm  plays  no  part  in  this  model  for  probability 
of  detection. 
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NOTABLE  PROGRAM  CAPABILITIES 


The  nominal  first-convergence-zone  prediction  capabilities  of  NUC  models 
ancestral  to  LORA  have  been  extended  to  allow  performance  predictions  for  a 
maximum  of  five  convergence  zones  and  bottom  bounces. 

Previous  NUC  active-sonar  models  apply  only  to  sonars  mounted  on  surface 
ships.  The  LORA  model  allows  any  source  and  target  depths. 

Propagation-loss  predictions  are  useful  byproducts  of  the  performance  predic- 
tions. Active  and  passive  propagation  losses  may  be  calculated  for  a maximum  of  30 
convergence  zones  or  bottom  bounces. 

The  sound-speed  profile  is  automatically  curve-fitted  by  computer  routines 
designed  by  C.  Bartberger.  who  used  the  mathematical  techniques  of  Pedersen  and 
Gordon.  Curve-fitting  the  sound-speed  profile  reduces  false  caustics  introduced  into 
the  propagation  loss  when  the  sound-speed  profile  has  discontinuous  derivatives. 

Multi-ping  reverberation  is  available;  that  is.  reverberation  caused  by  pings  pre- 
ceding the  current  ping  is  included  in  the  total  reverberation  level.  Single-ping  reverb- 
eration is  the  standard  output:  multi-ping  reverberation  for  a maximum  of  nine  pings 
is  optional. 

If  the  user  wishes,  reverberation  effects  can  be  reduced  because  of  the  doppler 
motion  of  the  target.  The  user  can  enter  doppler  gain  or  use  the  internal  calculations 
to  include  the  doppler  effect. 

Five  optional  probability-of-detection  models  are  available,  depending  on  the 
user’s  choice  of  detector  characteristics  and  his  assumptions  of  signal  degradation. 

The  computer  program  can  run  nearly  without  failure.  It  has  been  run  for 
thousands  of  cases.  An  example  of  its  capability  is  its  processing  of  240  sequential 
data  sets  in  the  same  program  execution;  that  is.  the  program  read  new  data  240 
times,  produced  performance  predictions,  and  terminated  normally.  The  total  cost 
for  that  set  of  runs  was  S60  (UNIVAC  1110  computer  charges),  of  which  S30  was 
for  paper.  The  LORA  model  and  supporting  computer  program  can  thus  give  the 
sonar  system  designer  the  ability  to  compare  the  performance  of  active-sonar  systems 
inexpensively. 

The  LORA  program  requires  38.000  words  of  36-bit  core  storage. 
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RECOMMENDATIONS 


Volume  reverberation  may  be  estimated  more  closely  by  using  a scattering 
strength  which  varies  with  depth.  The  added  raytracing  would  increase  program  run- 
ning costs  substantially.  In  addition  an  algorithm  would  have  to  be  developed  which 
smooths  reverberation  peaks  due  to  raytracing  caustics.  Such  an  algorithm  would 
depend  on  propagation  mode  (that  is,  direct  path,  convergence  zone,  etc.),  scattering 
geometry,  beam  pattern,  and  the  intuition  of  the  designer  of  the  algorithm. 

Nonspecular  scattering  could  also  be  added  to  the  reverberation  model.  This  is 
another  costly  calculation,  and  the  user  of  the  program  should  be  allowed  to  switch 
to  a cheaper  reverberation  model  if  he  desires. 

The  subsurface  duct  model  could  be  improved.  The  current  model  uses  cylindri- 
cal spreading  of  sound  in  the  duct.  At  low  wavelength-to-duct-width  ratios  cylindrical 
spreading  does  not  apply  since  some  acoustic  energy  leaks  out  of  the  duct.  Normal- 
mode theory  can  give  accurate  predictions  for  propagation  loss,  but  normal-mode 
models  are  typically  expensive  to  run.  To  keep  computer  costs  down  a subsurface 
duct  model  should  be  derived  which  approximates  normal-mode  theory  results. 

The  new  model  should  incorporate  such  parameters  as  duct  width,  source  and  receiver 
positions,  frequency  of  the  propagating  sound,  and  some  indicator  for  duct  strength. 

The  recognition  differential  model  used  here  is  not  representative  of  some  active 
sonars.  A helpful  addition  to  the  program  would  be  the  ability  to  enter  recognition 
differential  versus  signal-to-noise  ratio  in  order  to  make  the  program  more  sensitive 
to  certain  active  sonars. 
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Appendix  A 

LOGICAL  STRUCTURE  OF  THE  LORA  COMPUTER  PROGRAM 


The  LORA  computer  program  contains  three  major  computational  sections: 
environmental,  reverberation,  and  detection.  Figure  A-l  shows  the  computational 
structure  of  the  LORA  program.  All  inputs  are  entered  through  SUBROUTINE 
READ1N.  The  input  parameters  are  flagged  to  indicate  if  they  belong  to  environ- 
mental, reverberation,  or  detection  categories.  The  flag  "LOOP"  in  Fig.  A-l  is  set 
in  READ1N.  depending  upon  which  inputs  are  entered.  Table  B-l  in  Appendix  B 
gives  the  grouping  of  the  input  parameters,  environmental,  reverberation,  or  detec- 
tion. If  the  user  enters  only  reverberation  and  detection  parameters,  the  whole 
environmental  section  will  be  skipped  automatically,  without  user  specification, 
eliminating  redundant  calculations. 

OU  1 P and  OUTR  in  Fig.  A-l  are  input  parameters  set  by  the  user  if  he 
desires  special  outputs  for  propagation  loss  and  reverberation,  respectively.  The 
outputs  tor  the  detection  calculations  include  summaries  of  propagation  loss  and 
reverberation  which  are  sufficient  for  active  sonar  performance  prediction  work. 
However,  the  user  may  want  the  optional  passive  coherent  and  incoherent  propaga- 
tion loss  outputs  or  the  separated  components  of  reverberation  surface,  bottom, 
and  volume  as  well  as  the  relative  contribution  to  reverberation  from  previous 
pings.  Use  ol  OU  I R and  OU1P  also  allows  the  user  to  skip  the  detection  or 
reverberation  calculations  il  lie  is  interested  only  in  reverberation  or  propagation 
loss  outputs.  II  the  user  does  not  >et  OUTP  or  OUTR,  the  program  sets  them  to 
zero,  and  the  detection  outputs  a e the  standard  outputs. 

With  successive  runs,  it  may  he  that  environmental  parameters  are  changed 
which  do  not  influence  the  sound  speed  profile,  so  that  it  is  unnecessary  to  curve- 
lit  the  profile  lor  every  run.  The  dotted  line  in  Fig.  A-l  indicates  that  the  sound- 
speed  profile  calculations  are  considered  to  be  part  of  the  environmental 
calculations. 

The  environmental  calculations  in  l ig.  A-l  include 

• Searching  the  sound-speed  profile  for  surface  and  subsurface  ducts  and 
applying  the  AMOS  eom.tions  if  applicable. 

• C urve-fitting  the  sound-speed  profile  with  continuous-gradient  curve 
segments  after  the  profile  has  been  corrected  for  earth  curvature. 
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• Ray  tracing  to  obtain  path  segments  between  source  and  target,  source 
and  surface,  source  and  bottom,  and  source  and  volume  scattering  depths. 

• Assembling  the  path  segments  to  obtain  a multitude  of  different  paths 
for  the  purpose  of  computing  propagation  loss  and  reverberation 
components. 

• Interpolating  the  ray  paths  to  obtain  two-way  times  and  propagation 
losses  corresponding  to  integer  kiloyard  steps  in  horizontal  range  between 
source  and  target. 

• Obtaining  one-way  passive  coherent  and  incoherent  propagation  losses 
(optional). 

The  reverberation  section  performs  the  following  computations: 

• Interpolating  arrays  of  reverberation  components  (stored  during  ray  path 
computations)  corresponding  to  the  two-way  times  to  the  target  for 
integer  kiloyard  horizontal  ranges. 

• Grouping  (by  incoherent  summation)  the  effects  of  surface,  bottom,  and 
volume  reverberation. 

• Calculating  total  single-ping  reverberation  for  two-way  travel  times  to 
the  target. 

• Calculating  total  reverberation  from  time-delayed  previous  pings. 

• Calculating  total  reverberation  from  all  pings  for  the  current-ping  two- 
way  times  to  the  target. 

The  detection  section  calculates  and  prints  the  following  tabular  outputs  (see 
the  pages  in  Appendix  D labeled  “Active  Performance  Prediction”): 

• A heading  and  a brief  statement  of  the  propagation  mode. 

• A listing  of  the  values  of  the  single-value  inputs,  with  a trailing  letter  to 
indicate  status.* 

• A listing  of  the  values  of  certain  important  internally  calculated 
parameters. 


♦The  status  symbols  indicate  the  following  (see  p.  95  ff): 

U user  has  specified  this  input  value  on  this  tun. 

Blank  set  by  user  on  a previous  run;  unchanged  on  this  run. 

I typical  value  set  in  absence  of  user  specification. 

C calculated  internally  by  the  program  at  the  user’s  request. 

W wrong  value  specified  by  user.  (Out  of  range  of  allowable  values.  The  printed  value  is  either  a 
typical  or  an  internally  calculated  value.  I 
N not  an  input;  internally  calculated  parameter. 
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and  under  the  following  headings: 

• RANGE,  all  integer  kiloyard  ranges  to  the  target  for  which  the  current 
propagation  mode  is  applicable. 

• TIME,  corresponding  two-way  travel  times  for  the  rays  which  give  mini- 
mum loss  to  the  target. 

• LOSS,  minimum  two-way  propagation  loss  to  the  target,  including  beam 
pattern  corrections. 

• ANGLE,  starting  angles  for  the  rays  at  the  source  depth  which  give 
minimum  loss  to  the  target.* 

• SIGNAL,  signal  level  of  the  return  signal  in  the  receive  band  (as  measured 
at  the  output  of  the  beamformer). 

• REVERB,  reverberation  level  before  doppler  correction  (doppler  gain 
“GDOP”  is  listed  with  the  single-value  inputs). 

• REV  LIM,  a “YES”  or  “NO”  to  indicate  if  the  masking  noise  is  reverb- 
eration limited  (that  is,  doppler-corrected  reverberation  level  exceeds 
ambient  noise  level  in  the  receive  band  by  3 dB). 

• S/MN,  signal-to-masking-noise  ratio  in  the  receive  band. 

• EXCESS,  signal  excess. 

• PROB  DETECT,  probability  of  detection. 


*This  is  useful  information  for  aiming  the  beam  pattern  on  successive  runs.  No  values  for  ANlil  1 are 
given  in  the  surface-duct  or  subsurface-duct  modes  of  propagation. 
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Appendix  B 

INPUTS  TO  THE  LORA  COMPUTER  PROGRAM 


Table  B-l  contains  a list  ot  the  inputs  to  the  LORA  computer  program.  Each 
input  is  designated  by  a name  that  corresponds  to  a parameter  used  in  other  com- 
puter programs  at  NUC  These  names  are  listed  in  the  column  "NUC  Names.” 
Modifications  to  the  Nt  ( names  were  required  by  users  at  the  Na\al  Intelligence 
Support  Center  (NISO,  because  their  computer  word  contains  lour  characters  rather 
than  six  as  at  NUC  Thus,  new  and  sometimes  more  mnemonic  and  acronymistk 
names  were  developed.  These  occur  in  I able  B-l  in  the  column  under  "N1S(  N.ime^ 
Input  parameters  identified  bv  either  set  of'  parameter  names  are  acceptable  to  ilk 
program  without  user  specification.  Output  labels,  however,  are  identified  In  using 
either  ot  the  data  literals.  "NI  C or  "NISC  ( see  1 able  B 2 1 to  denote  the  appropri 
ate  set  of  parameter  names.  The  column  under  “Description”  in  Table  B-l  contains 
brief  definitions  of  the  input  parameters.  “Range”  is  the  range  of  values  a given 
parameter  is  allowed.  Any  input  value  outside  this  range  is  either  set  to  a " I epical 
Value"  or  calculated  in  the  program. 

Table  B-2  is  the  list  ot  data  literals  used  to  exert  specific  control  over  the 
program.  The  user  must  type  these  one  to  a card  (line,  record). 

Not  all  of  the  input  parameters  need  to  be  entered  by  the  user  In  fact  only 
one  input  parameter  is  required  to  be  entered  by  the  user,  the  sound-speed  profile. 
“SSP."  Typical  or  calculated  values  will  be  produced  for  every  other  input  param- 
eter that  is  not  entered.  Of  course  the  outputs  may  have  little  significance  to  the 
user  who  does  not  exercise  his  prerogative  to  enter  his  own  data. 


Table  B-l.  LORA  inputs. 


ENVIRONMENTAL  PARA  ML  FIRS 


3 1 

ft 

N1SC 

Name 

Description 

Units 

Range 

Typical  Value 

OUTP 

OUTP 

Mag  for  spccidi  output  of  propagation 
losses  (=0.  no  special  output;  =1,  print 
and  continue  to  reverberation  calcula- 
tions; =2,  print  and  go  back  to  read  new 
data) 

0 to  2 

X 

0 , 

ZX 

SU 

Source  depth 

ft 

1 to  40.000 

20 

ZTG 

TD 

Target  depth 

it 

1 to  40,000 

50 

i ri  y 

1 K1  Q 

Frequency  of  transmitter 

kHz 

0.099  to  25 

3.5 

LWA 

VVAV1 

Wave  height 

ft 

0 to  20 

4 

VWI 

WIND 

Wind  speed 

knots 

0.1  to  35 

15 

SALT 

SALT 

Salinity 

°/oo 

33  to  37 

35 

US 

SLPB 

Siunal  loss  per  bounce  at  the  surface 

dB 

-2  to  10 

3,  calculated  it  negative 

ZBM 

BD 

Bottom  depth  (if  different  than  the  larg- 
est depth  in  the  sound-speed  profile) 

ft 

1 to  40,000 

largest  depth  in  SSP 

KIT  I S 

NBBB 

Maximum  number  of  bottom  bounces 

1 to  30 

i 

ZONES 

ZONE 

Maximum  number  of  convergence  zones 

- 

1 to  30 

i 

PH  DX 

XDL 

Transmit  depression-elevation  angle  of 
the  main  axis  of  the  beam  pattern 

deg 

-85  to  85 

5 

Pll  UR 

RDE 

Receive  depression-elevation  angle  of  the 
main  axis  of  the  beam  pattern 

deg 

-85  to  85 

5 

DLLPX 

XV  BM 

Effective  vertical  beamwidth  transmit 

deg 

0.1  to  180 

10 

DILPK 

RVBM 

Effective  vertical  beamwidth  receive 

deg 

dB/yd* 

dB/yd“ 

0.1  to  180 

10,  set  to  180  if  BEAMR 
is  omnidirectional 

MUB 

BSS 

Bottom  backseat tering  strength 

-40  to  -5 

-27 

MUV 

VSS 

Volume  backscattering  strength  of  the 
water  column 

-80  to  -30 

-49 

ZU 

SLD 

Scattering  layer  depth  (water  column 
depth) 

ft 

0 to  10,000 

2000 

LAT 

DUMP 

LAI 

DUMP 

Latitude 

flag  for  diagnostic  printout  (=0,  no 
diagnostic  printout;  =-l , print  ray  trac- 
ing outputs;  =1,  print) 

deg 

-90  to  90 

0 

SSP 

SSP 

Sound-speed  profile  array  of  alternat- 
ing depth  and  sound-speed  values  or 
depth  and  temperature  values 

ft,  ft  /see 

m.  m/sec 

ft,  °F 
m,  °C 

4750  to 
5160** 
1440  to 
1590** 
30  to  85** 
0 to  30** 

no  typical  value 

SSP  must  be  input  for 
program  to  run 

Aim 

BRLP 

Bottom  reflection  loss  profile  array  o( 
alternating  grazing-angle  and  bottom- 
reflectiondoss  values 

deg,  d B 

not  checked 

zero  loss  at  am  grazing 
angle 

♦The  input  parameter  is  automatically  set  to  this  value  in  absence  of  user  specification.  K ontd) 

♦♦This  range  of  values  applies  only  at  the  ocean  surface,  that  is.  the  second  value  in  the 
input  array  al  ter  elements  have  been  ordered  with  respect  to  increasing  depth. 


Tabic  B-l.  (Contd) 


NUC 

NISC 

Name 

Name 

Description 

Units 

Range 

Typical  Value 

BEAMX 

XV  RP 

Transmit  vertical  response  pattern  (i.e.. 

deg,  dB 

not  checked 

omnidirectional  (i.e.,  0 

beam-pattern)  array  of  alternating 

dB  for  all  angles);  if 

angle  and  response  values,  where  response 

the  array  contains  one 

is  given  as  positive  dB  down  from  the 

element  only,  the 

maximum  response,  defined  as  0 dB  at 

sin  (x)/x  response  pat- 

an  angle  of  0 degrees 

tern  for  a continuous 
line  is  calculated 

BIAMR 

RVRP 

Receive  vertical  response  pattern 

Same  as  for  BEAMX 

REVERBERATION  PARAMETERS 

OUTR 

OUTR 

Flag  for  special  output  of  reverberation 
values  (bottom,  surface,  and  volume 
components,  total  contribution  from 
each  ping,  and  overall  total  reverbera- 
tion for  two-way  time  to  target  corre- 
sponding to  integer  ranges  in  kyd  (=0, 
no  special  printout;  =1 , prints  and  con- 
tinues to  detection  section;  =2,  prints 
and  goes  back  to  read  new  data) 

0 to  2 

0 

MO0E 

MODE 

Mag  to  indicate  sonar  operation  mode 
(=1,  near-surface  or  direct  path;  =2, 
bottom  bounce;  =3,  convergence  zone, 
surface  bounce,  or  reliable  acoustic 
path) 

1 to  3 

1 

SL 

SL 

Source  level 

dBre  ljUb 

40  to  16U 

148 

DELTH 

HBM 

Effective  horizontal  beamwidth 

deg 

0 to  360 

10 

MPR 

MPR 

Number  of  pings  for  multi-ping 

- 

1 to  9 

1 

reverberation 

TBP 

TBP 

Time  between  pings 

sec 

1 to  90 

20 

PULSI 

PULS 

Pulse  length  to  help  determine  reverbera- 

sec 

ur4  to  s 

0.1 

tion  backscattering  areas 

DETECTION  PA  R A M ETERS 

D! 

Dl 

Directivity  index  (used  only  to  deter- 

dB 

0 to  70 

0 

mine  the  amount  of  discrimination 
against  isotropic  noise  due  to  the  direc- 
tivity of  the  beam  pattern.  If  NOUT  is 
known.  Dl  need  not  be  input) 

BWS 

BWS 

Bandwidth  of  the  return  signal 

Hz 

1 to  1000 

100 

BVVR 

BWR 

Receiver  bandwidth 

Hz 

1 to  1000 

BWS 

N1N 

NIN 

Omnidirectional  noise  spectrum  level  at 

dBre  ljUb 

-80  to  0 

-45 

the  receiving-array  input 

NOUT 

NOUT 

Noise  band  level  at  the  receiving  array 

dBre  1 fib 

-10S  to  29 

-61.2 

output  or  signal  processor  input 

TVC 

LVC 

Target  closing  velocity  (target  speed  with 

knots 

-60  to  60 

0 

respect  to  the  ocean  along  the  range 
track  between  source  and  target)  used 
in  estimating  doppler  gain  of  signal  over 
reverberation. 

K'ontd) 
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Table  B-l.  (Conld) 

NIX' 

Name 

NISC 

Name 

Description 

Units 

Range 

Typical  Value 

GDOP 

GDOP 

Doppler  gain  (signal  enhancement  over 
reverberation  due  to  shift  of  reverbera- 
tion band  away  from  return  signal  band 
it  target  is  moving  relative  to  environ- 
ment) 

dB 

-2  to  300 

calculated  using  1 VC  . 
negative  value  results 
in  calculation 

15 

TGS 

TGS 

Target  strength 

dB 

0 to  40 

KL)N 

RDN 

Recognition  differential  for  ambient- 

dB 

-20  to  20 

noise-limited  detection  (corrected  by 
program  for  reverberation-limited 
detection 

0 to  I 

0 

CUP 

CUP 

1 lag  to  indicate  2-dB  clipping  of  signal  at 
processor  (=0,  no  clipping:  =1,  clipping) 

0.5  to  16 

SIGMA 

SIG 

Standard  deviation  of  signal  differential 

dB 

6 

to  be  used  with  signal  processing  model 
given  by  SP^5 

1 1 

0.001 

PI  A 

PI  A 

Probability  of  false  alarm  for  one  resolu- 

lCTUto0.1 

celi  on  some  display 

1 to  5 

5 

SP 

SP 

flag  to  denote  signal  processing  model 

= 1 . linear  correlator,  known  signal: 
=2,  linear  correlator,  constant  signal: 
=3,  quadrature  correlator,  constant- 
amplitude  signal: 

-4.  quadrature  correlator,  random- 

amplitude-and-phase  signal 
=5.  energy  detector  for  gaussian  signal 

7(> 


Table  B-2.  Data  literals. 


Data  Literal 


Resulting  Control 


RUN 

XQT 


XLQ 


P A US  l 


END 

SlOP 

IN  D-D  A I A 


Calculations  are  performed  for  the  current  values  of  the  input  parameters  and  control  is 
returned  to  SUBROUTINE  R1  ADIN  for  more  inputs. 

Same  as  RUN. 

Same  as  RUN. 

Same  as  RUN. 

Program  execution  is  immediately  terminated  and  control  is  returned  to  the  executive. 
Same  as  END. 

Same  as  END. 


Ill  ADI  R 
NUC 

NISC 

PASY 

PASN 

COMY 

COMN 


Allows  the  alphanumeric  string  on  the  next  card  to  be  used  for  headings  on  output  pages. 

Naval  Undersea  Center  input  parameter  names  are  used  to  label  the  outputs  (assumed  it 
NISC"  is  not  entered  by  the  user). 

Naval  Intelligence  Support  Center  input  parameter  names  are  used  to  label  the  outputs. 

One-way  passive  propagation  loss,  including  receiver  response  pattern  and  multipath  ray 
summations,  is  calculated  (OUTP-1  or  2 is  needed  to  display  the  values). 


Defeats  PAS'!  n it  has  been  entered. 


One-way  communications  propagation  loss,  including  transmit  and  receive  response  pat- 
terns and  multipath  ray  summation,  is  calculated  (OUTP=l  or  2 is  needed  to  displav  the 
values). 

Defeats  COMY  if  it  has  been  entered. 


COHY  Coherent  propagation  losses  are  calculated,  as  well  as  the  incoherent  losses,  for  either  the 

passive  or  communications  options.  Semicoherence  is  also  included  but  not  recommended 
because  of  the  artificiality  of  the  calculations. 

CONN  Defeats  CONY  if  it  has  been  entered. 


TTY  Reduces  output  for  demand  terminal  usage. 

NOTT  Defeats  TTY  (assumed  if  TTY  has  not  been  entered);  normal  outputs  and  page  headings 

are  printed. 
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Appendix  C 
INPUT  FORMAT 


The  LORA  program  accepts  all  inputs  through  SUBROUTINE  RLADIN. 
RFADIN  contains  only  two  read  statements.  The  input  characters  are  scanned  in  a 
manner  similar  to  NAMELIST  to  determine  if  a string  of  characters  is  an  input  name 
or  an  input  parameter  value  (a  number).  READIN  does  not  use  NAMELIST  for  the 
following  reasons: 

• It  is  important  to  keep  track  ot  the  variables  which  are  changed  from  run 
to  run  in  order  to  reduce  redundant  calculations. 

• Column  I on  the  input  record  can  now  be  used. 

• Data  literals  which  have  no  numerical  value  can  be  entered. 

The  user  enters  single-value  variables  by  typing  the  name  of  the  variable,  followed 
by  a space  or  an  equal  sign  and  then  by  the  numerical  value  in  F. !.  E.  or  D 
format.  All  input  values  are  converted  to  single-precision  Boating  point  numbers. 

There  can  be  no  imbedded  blanks  in  either  the  name  or  value,  although  there  may  be 
more  than  one  blank  separating  the  name  and  value.  A name  and  its  value  must  be  on 
the  same  line.  If  more  than  one  name-value  set  occurs  on  a line,  the  sets  are  to  be 
separated  by  a comma  or  a slash  "/'  • There  is  no  need  to  use  a comma  or  a slash 
as  a separator  between  sets  occurring  on  different  lines.  The  parameters  can  be  input 
in  any  order  or  not  input  at  ail.  as  the  user  desires. 

The  user  enters  arrays  in  a similar  fashion  to  his  entering  single-value  variables, 
by  lirst  typing  the  name,  then  a space  or  equal  sign  and"  then  the  array  values  in 
F.  I.  E.  or  D format.  Successive  array  values  are  separated  by  one  or  more  spaces.  If 
the  array  is  too  long  for  one  line,  as  many  lines  as  are  required  can  be  used  by  typing 
the  continuation  symbol  “*  ”,  that  is.  an  asterisk  and  a space,  as  the  first  characters 
on  the  successive  lines.  An  array  is  terminated  by  a comma  or  a slash  or  by  starting  a 
new  line  without  the  continuation  symbol. 

Data  literals  are  names  which  are  not  followed  by  values.  To  enter  the  name  is  to 
set  a flag  or  to  foliow  some  instructions  before  reading  the  next  input.  When  one  data 
literal  is  read,  the  line  is  automatically  terminated.  Data  literals  must  be  entered  one 
to  a line  and  be  the  last  datum  on  that  line.  Any  other  data  on  the  line  following  the 
data  literal  will  be  ignored. 
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Figure  C-l  illustrates  some  of  the  possible  ways  data  may  be  entered.  ( olumn  1 
is  indicated  by  the  arrow.  The  line  numbers  are  merely  for  reference  in  this  discus- 
sion. Line  1 in  Fig.  C-l  shows  several  examples  of  separators  and  spacings  and  the 
use  and  nonuse  of  the  equal  sign.  Line  2 contains  the  data  literal  "PASY.  which 
will  cause  a tlag  to  be  set  so  that  passive  propagation  losses  will  be  calculated.  Line  3 
is  another  data  literal.  "HFADFR.”  which  causes  the  next  card  to  be  read  as  heading 
information.  Line  5 shows  how  an  array  may  ie  entered.  A user  at  a demand  termi- 
nal may  make  a mistake.  If  so.  he  can  reenter  the  same  parameter  as  on  Line  tv  The 
latest  entry  ot  any  parameter  supersedes  any  value  it  may  have  been  assigned  by  pre- 
vious entries.  The  user  need  not  specify  the  number  of  values  in  an  array;  the  program 
automatically  counts  them.  Lines  8-10  illustrate  the  entering  ol  values  in  array 
"AHB.  employing  continuation  cards.  On  card  1 1 "I(iS~20  follows  a data  literal, 
so  it  is  ignored.  Thus.  r(iS=15.  On  Line  I 2 PFA  is  assigned  value  in  1 format  and  on 
Line  13  in  1)  format.  The  data  literal  on  Line  1 2.  "Rl  V will  cause  immediate 
termination  of  the  RFADIN  subroutine,  calculations  will  be  performed  with  the  cur- 
rent data  set,  and  control  will  be  returned  to  RFADIN  to  read  the  new  value  tor  PFA 
on  Line  13.  The  "RUN"  on  Line  14  will  again  cause  calculations  to  be  performed, 
this  time  only  in  the  detection  section,  since  PFA  is  a detection  parameter.  Line  1 s 
ends  the  program  with  "FND."  and  control  is  returned  to  the  executive. 

Col.  1 

1.  FREQ  1..  ZX  20/ZTG  300  .PHDX  5 .PHDR  .5 

2 PASY 

3 HEADER 

4 THIS  IS  A HEADING 

5 SSP  0 5000  100  5001.8  3000 

6 SSP  0 5000  100  5001.8  3000  4950  10000  5050 

7 AHB  0 5 5.0  8.5  10.  ’?  1 

8 ' 20  13  40  14 

9 * 55  ]5 

10  • 90  15  I 

11  TGS-15.  COHY.TGS  20 

12  PFA  1.E-5,  RON 

13  PFA  I D-4 

’ 14  RUN 

, 15  END  1 

1 ij;urc  C-l.  1 sample  of  some  options  in  the  input 

format.  a 

A data  set  equivalent  to  that  represented  in  Fig.  C -l  is  shown  in  f ig.  C-2. 
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Col.  1 

l 

HEADER 

THIS  IS  A HEADING 

FREQ-1  ZX-20,  ZTG-300.  TGS-15 

PFA-.00001,  PHDX-5,  PHDR--5 

PASY 

COHY 

SSP-0  5000  100  5001.8  3000  4950  10000  5050 
AHB-0  5 5 8.5  10  12  20  13  40  14  55  15  90  15 
RUN 

PFA-.0001,  RUN 
END 

Figure  C-2.  Rewritten  form  of  the  data  in  figure  C-l . 


Appendix  D 


EXAMPLES  Of  RUNS  FOR  THE  LORA  COMPUTER  PROGRAM 


I his  Appendix  contains  examples  of  inputs  and  resulting  outputs  of  the  LORA 
computer  program.  To  keep  the  report  unclassified,  no  attempt  is  made  to  make  the 
system  and  environmental  inputs  correspond  to  data  or  theoretical  values,  with  the 
exception  that  the  two  sound-speed  profiles  are  from  unclassified  measured  data.* 
These  profiles  are  depicted  in  Fig.  D-l  and  D-2  and  tabulated  in  Tables  D-!  and  D-2. 
In  the  computer  runs  the  typical  values  for  sound  speed  are  used  rather  than  the  aver- 
age values. 

1 he  transmit  and  receive  beam  patterns  are  the  internally  calculated  response 
functions  for  a continuous-line  array,  with  the  half-power  beamwidths  defining  the 
length  ol  the  line  lor  the  given  frequency.  It  is  important  to  remember  that  the 
propagation  loss  values  include  beam  response  degradations. 


Depth, m 


Sound  Speed,  ni/s 


Typical 


Average 


0 

1511.8 

1512.46 

10 

1512.0 

1512.56 

20 

1512.2 

1512.74 

30 

1512.3 

1512.85 

50 

1512.8 

1513.08 

75 

1513.5 

1513.45 

100 

1514.1 

1513.86 

125 

1514.4 

1514.1 1 

150 

1514.5 

1514.45 

200 

1514.7 

1515.44 

250 

1515. 1 

1515.92 

300 

1515.8 

1516.23 

400 

1516.8 

1516.95 

500 

1 5 1 7.9 

1517.82 

600 

1519.0 

1519.00 

700 
800 
900 
1 000 
I 1 00 


1520.2  1520,25 

1521.3  1521.47 
1522.83 
1524  41 
1526.01 


1200 
1 300 
1400 
I 500 
1750 


1527.64 
1 

1530.88 

1532.53 

1536.64 


2000 

2500 

3000 


1540  78 
1549.20 
1557.87 


Table  l)-2.  Sound  speeds  at  standard  depths  for  typical  profile  and 
average  profile.  Transition  Zone,  winter. 


Depth, m 

Sound  Speed,  m/s 

Typical 

Average 

0 

1488.0 

1490.36 

10 

1488.2 

1490.69 

20 

1488.6 

1490.87 

30 

1488.8 

1490.98 

SO 

1488.8 

1491.05 

75 

1489.0 

1491.13 

100 

1489.2 

1491.21 

125 

1489.3 

1490.95 

150 

1489.6 

1490.36 

200 

1489.5 

1488.31 

250 

1485.5 

1485.75 

300 

1479.9 

1483.23 

400 

1477.2 

1479.46 

500 

1476.6 

1477.38 

600 

1476.6 

1476.75 

700 

1476.7 

1476.81 

800 

1477.2 

1477.33 

900 

1477.8 

1477.95 

1 000 

1478.6 

1478.69 

1 100 

1479.4 

1479.64 

1200 

1480.5 

1480.67 

1300 

1481.7 

1481.82 

1400 

1483.0 

1482.94 

1500 

1484.4 

1484.20 

1750 

1487.9 

1487.55 

2000 

1491.3 

1491.12 

2500 

1498.57 

3000 

1506.52 

4000 

1523.85 
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LOR*  --  LONG  RANSt  ACTIVE  PERFORMANCE  PREDICTION  PROGRAM 

RESPONSIBLE  PERSON;  DAVID  6,  HOFFMAN 

naval  undersea  center 
code  SOM 

SAN  DIEGO,  CA  92132 

TELEPHONE  NUMBERS;  (71m)  225-2316 

AUToVON  933-2316 


• inputs,  Run  i • 


HEADER 

MEDITERRANEAN  WINTER  SOUND  SPEED  PROFILEi  NUC  TN  1006,  P,  23 

SSP  0 1 S 1 1.8  >0  1 S 1 2 • 0 20  IS12.2  30  15)  2.3  SO  1S12.B  75  1S13.5  100  ISlM.l 

• 125  1 5 1 M , M 150  15IM.5  20C  IS1M.7  250  1515. I 300  1515,8  MOO  |5!6.8  500  1517.9 

• 600  1519,0  700  1520.2  800  1521.3  900  1522.83  1000  1S2M.M1  1100  1526,01 

• 1 200  1 527. 6M  1300  1529,26  1 MOO  1530.88  ISUO  ) 532.53  1 750  1536. 6M 

• 2000  ISM0.78  2 SO j 15M9.20  3000  1557.87 

2 B M 60M0 

BEAMX  0.  B E A MR  0 

AHB  ■ 0 1 10  2 50  6 70  7 90  7.5 

2*«M0Q|  ZTG«10U0,  FR£u»2,  L«A*S,  VAI-IS,  HS»-1,  20N£S»2 
DELPX-10,  dELPR*13,  DELTm»10,  LAT=3S,  OUTR-1, 

SL  a 1 3 5 i MPR*2,  TuPaltl,  TyCal2,  BRRalo,  OUTPal,  COHY 
MODE  1 . run 


• CALCULATIONS.  tun  l . 


BEAMX  lS  THE  SlNlXl/X  RESPONSE  FUNCTION  FOR  A CONTINUOUS  LINE  ARRAY 

beamr  is  thl  sInixi/a  Response  function  for  a continuous  line  array 


THE  LETTERS  FOLLOnlNu  Tme  NUMBERS  IN  TmE  HEADING  ON  THE 
PERFORMANCE  PREDICTION  PAGE  SIGNIFY  THE  FOLUOoING, 

U USER  HAS  SPECIFIED  ThIS  INPUT  VALUE  ON  THIS  RUN 

blank  set  bY  user  or.  a previous  run 

T TYPICAL  VALUE  set  IN  ABSENCE  of  USER  SPECIFICATION 

c calculated  Internally  by  program  (user  option) 

A aRoNG  value  SPECIFIED  by  USER  (OUT  OF  RANGE  OF  ALLOW- 

ABLE values.  the  PRlNTtD  VALUE  IS  A TYPICAL  VALUE.) 
n not  an  Input,  calculateu  internally. 


SOUND  S P E E 0 PROFILE  A F T t R earth  curvature  correction  A no  CURVE  fitting 

SSP  SEGMENT  PARAMETERS 
(JASA,  VOL.  61.  no.  2,  PP.  819-638) 


LATER 

depth 

SOUND  SPEED 

number 

C F T ) 

(FT/SEc  ) 

C A • • 2 

K A 

ZA 

1 

. DO  3 

8959. 780292 

.250837  14-08 

-.11 72266-ue 

. 6098082*06 

2 

328. 1 8 1 

8967. 157258 

.286817 1408 

-.51 S07 1 2-^7 

. 6 | 25656*03 

3 

372.683 

8967.86S836 

• 2869760*08 

-.S8223SS-U8 

. 72S731 1*03 

6 

892.136 

8968. 876877 

.2865123*08 

.11811 10-08 

-. 6572257*03 

5 

1312.839 

8976.8151 17 

.2637  1 36*08 

-.9586963-10 

.2738071 *05 

6 

1591.778 

8979.68962  1 

.2752025+08 

-.6010778-10 

.6633693*05 

7 

1966.612 

8983.980388 

. 2863969*08 

.6938600-09 

-. 1 572666*08 

8 

3609.212 

SoC  7 .292068 

.2665761*08 

-.2061638-09 

.2111798*05 

9 

3820. J08 

5010.857818 

.2676960*08 

-. 1 9595 1 3-u9 

. 22u823 1 *05 

10 

8265. 857 

5018.280877 

.9853603*07 

.2637853-1 1 

-.5019662*06 

1 1 

6Q8Q.9B9 

5087.877286 

.0000000 

. DQuQCQO 

.0000000 

89 


ACTIVE  T60-..AV  LOSS,  PAGE  I 


MEDITERRANEAN  «INTER  SOuNO  SPEED  PROFILE,  NUC  TN  1006,  P.  23 


THE  LOSS  ENTRIES  INCLUDE  TrtO-WAT  PROPAGATION  LOSS 
AND  THE  TRANSMIT  AND  RECEIVE  BEAM  PATTERN  RESPONSES* 


RANGE 

MINIMUM 

BOTTOM 

DIRECT  PATH 

IKTD) 

TaO-PAV 

BOUNCE 

CONVERGENCE 

LOSS 

LOSS 

2 0 n £ LOSSES 

I 

1 3 I . 1 8 

212.02 

131.18 

2 

132.83 

211.91 

132.83 

3 

139.99 

203.76 

139.99 

9 

19S.SG 

215.79 

195.50 

S 

199.6b 

201,13 

199.65 

6 

1 S3  * 0 7 

203. 1 1 

153.07 

7 

1 5 6 . 0 8 

215.91 

1 5 6 . u 8 

8 

1 SB . 60 

202 .68 

153.60 

9 

160.63 

196.73 

163.63 

10 

162.19 

197.25 

162.19 

1 1 

163. 9U 

201.59 

163.90 

12 

1 69 . SS 

2 10.15 

169.55 

13 

165.66 

215.32 

165.66 

1 9 

166.72 

222.83 

166.72 

! S 

167.79 

216.21 

167.79 

1 * 

169.78 

209 . 1 8 

169.78 

17 

171.82 

205.09 

171.62 

1 8 

1 7 3 . 7d 

202.73 

173.70 

19 

17b.  IS 

201.99 

175.15 

20 

176.38 

2uC .65 

176.38 

21 

177.57 

20C • 9 8 

177.57 

22 

1 78.72 

2 0 C . 5 8 

178.72 

23 

177.80 

2 00*90 

177.80 

2 H 

177.99 

201.53 

177.99 

2 5 

1 77.9S 

202 . 3 1 

177.95 

26 

1 77.S9 

203.29 

177.59 

27 

177.69 

209.32 

177.69 

28 

1 66,68 

205 . 59 

168.68 

29 

i 78.S6 

206.92 

178.56 

30 

189,66 

208.98 

189.66 

31 

1 87.09 

2 1C .22 

187.09 

32 

1 86.32 

2 12.17 

186.32 

33 

1 65. 75 

219.38 

185.75 

39 

165.99 

2 16.97 

185.99 

3 S 

1 91 . 91 

22C.  1 9 

191.91 

34 

192.78 

229,99 

192.78 

37 

| 93.69 

227.26 

193.69 

38 

1 99. 97 

229 . 99 

199,97 

39 

1 95. 89 

233.31 

195.89 

90 

197.56 

237.99 

197.56 

9 1 

(99. 16 

267 . 22 

199.16 

92 

200.66 

273.91 

200 .66 

93 

202.07 

• CO 

202.07 

99 

203. 1 7 

• Oo 

203. I 7 

'9S 

209.09 

• 00 

209.09 

I 


ACTIVE  TAO--AY  LOSS,  PAGE  2 


m e d i terranEa.j  .* i„tep 

SOu^D 

SPEED  PRO!  1 LE , 

range 

MINIMUM 

BOTTOM 

DIRECT  path 

l kyd  ) 

T „ 0 - ..  A Y 

BOUNCE 

convergence 

LOSS 

LOSS 

ZOnE  losses 

96 

201.89 

•Os/ 

20 9.89 

97 

20b .72 

• U 0 

205. 72 

98 

20e.S2 

• 0 0 

206.52 

99 

207.30 

iUg 

207 . 3C 

SO 

207.72 

• w U 

207. 72 

SI 

2 u u » 9 6 

•uO 

200*96 

52 

206.5b 

• 0 0 

206.55 

53 

206.51 

• 00 

206.51 

59 

2 u 6 . 5 1 

• uO 

206.51 

SS 

206.57 

• Ci  j 

206.57 

56 

2 0 6 • 6 6 

• 00 

2 w 6 . 6 6 

57 

206.78 

• Wo 

206.78 

S8 

2U6 . 92 

*uC 

206.92 

S 9 

207 .06 

• *-  w 

207  . CB 

60 

2 0 7 , 2 d 

• o j 

207.23 

6) 

207.51 

• 0 0 

207.51 

62 

209.88 

• u j 

209.88 

63 

219.23 

• u 0 

2 19.23 

69 

215.96 

• 3 0 

2 15.96 

6 5 

2 1 o . 6 9 

• 00 

216.69 

66 

2 1 6 . 7 7 

• U w 

2 16.77 

67 

217.7b 

• 00 

2 17.75 

68 

2 2 u .09 

• SI  0 

22C . .9 

69 

22u .98 

• 0 0 

220.98 

70 

221 ,9u 

• 00 

221 .90 

7 1 

222.83 

• 0 0 

222.83 

72 

223.77 

• 00 

223.77 

73 

229. /l 

• 00 

229 . 7 1 

79 

225.66 

• 00 

225.66 

7 S 

2 2 6 . 6 2 

• 00 

226.62 

76 

227.59 

• GO 

227.59 

77 

226.56 

• 00 

228.56 

78 

229.55 

• 00 

229 . 55 

79 

237.57 

• 0 0 

237.57 

80 

23d . 6 1 

#00 

238.61 

8 1 

239.6b 

• GO 

239.65 

82 

29u . 70 

• 00 

290.70 

83 

271.99 

• 0 0 

271.99 

. 2 J 


passive  OnE-hAY  loss,  page  i 

MEDITERRANEAN  WINTER  SOUND  speed  profile,  nUC  TN  1006,  P.  23 

THE  loss  ENTRIES  include  ONE-WAY  PROPAGATION  LOSS 
and  the  receive  beam  pattern  response,  multipath  intensities 
ARE  SUMmEo.  SO  the  ONE-aAY  PASSIVE  LOSSES  ARE  NOT  IN  GENERAL 
ONE-HALF  OF  THE  TwO-WAY  ACTIVE  LOSSES. 


range 

I NCOhEREnT 

COHERENT 

SEmICOHERENT 

(KYD) 

1 - A A Y LOSS 

1-*Ay  LOSS 

1-WAY  LOSS 

l 

6S  • 53 

69.93 

69.93 

2 

6 6,36 

65.5  3 

65.52 

3 

69.93 

70.21 

70.20 

9 

72.75 

72  • 0 0 

72.80 

S 

79.79 

79.60 

79.61 

6 

76.95 

7 6.30 

76.31 

7 

77.90 

70.23 

78.19 

8 

79.08 

79.08 

79.08 

9 

8 u , 0 1 

8o  .60 

80.  SC 

10 

80.75 

82.95 

82.19 

I i 

0 1.35 

81.62 

81.56 

12 

81.99 

82.65 

82.99 

13 

82.52 

85.01 

89.90 

19 

83.06 

82.92 

82.59 

IS 

83. Sb 

82.29 

82.66 

16 

89.55 

85.96 

85.18 

17 

os  • s i 

86.77 

8 6.36 

lB 

0 6.37 

8 8 . 0 9 

87.99 

1’ 

87.09 

89.28 

88.95 

20 

87.67 

’3.29 

91.11 

21 

80.26 

8 8.65 

88.99 

22 

00.82 

89.19 

89.03 

23 

08.95 

95.26 

92.20 

29 

88.33 

87.61 

87.95 

2 S 

8 0.37 

89.39 

88.87 

24 

88.99 

88.86 

88.65 

27 

08.55 

88.91 

88.72 

28 

81.79 

80 . 82 

89.95 

29 

85.83 

89.69 

87.98 

30 

08.16 

88.20 

88.18 

31 

80.90 

97.99 

92.92 

32 

0 0.82 

85.52 

87.61 

33 

69.12 

’2.06 

90.19 

39 

09.93 

6 6.26 

88.39 

35 

91.52 

’3.01 

9 1.98 

36 

92.91 

89.50 

91.58 

37 

93.08 

’0.25 

92.33 

38 

95.33 

’3.93 

99.87 

39 

96.08 

’6.63 

96.21 

90 

97.03 

102.59 

98.15 

91 

97 . 80 

’9.18 

98.12 

92 

99.99 

’7.60 

99.60 

93 

100.73 

101.37 

100.82 

99 

|01 .38 

101.99 

101.95 

95 

I 0 | .85 

1 u2 . 9 2 

101.96 

I 


passive  one-aat  loss,  page  2 

MEDITERRANEAN  "INTER  sound  SPEED  PROFILE,  Nile  Tw  1006,  P,  23 


range 

I NcOmEKEnT 

coherent 

SEmICOHERENT 

1 K V 0 1 

l-*AY  LOSS 

i-aay  loss 

l-*AY  LOSS 

HA 

102.31 

101.02 

102.20 

H 7 

IG2.7S 

102.99 

102.76 

HB 

103.  IA 

103.12 

103. 1 A 

H 9 

103. SA 

103.HH 

103.56 

SO 

103.79 

1 oh • ss 

103.79 

Si 

96. 7B 

IOa.SH 

96.78 

52 

100.31 

102.91 

100.31 

S3 

101.31 

112.67 

101.31 

5 H 

1 0 1 .7A 

1 10. HB 

101.76 

55 

lOl.Hl 

10B.2A 

101  .HI 

SA 

102.09 

98.  AS 

102.09 

57 

102.31 

1 OH .79 

102.31 

SB 

102. HB 

1 00  • 8 H 

102. H8 

59 

|02. H9 

101.99 

102. H9 

AO 

I02.H3 

100. A9 

102. H3 

Al 

102. H9 

10H.9H 

102.  H9 

A 2 

(03.32 

101.08 

103.32 

A3 

1 OH  • A 7 

1 1 1 .70 

1 OH . 6 7 

AH 

105.13 

102. HI 

10S. 13 

AS 

10S.H9 

112.95 

105. H9 

AA 

105.81 

109. A9 

105.81 

A 7 

1 OA  • 2 3 

1 IH.77 

106.23 

A8 

1 08 . 03 

112.17 

108.03 

A 9 

1 OB  . A8 

1 0 A . H9 

108.68 

70 

109. 2A 

116.77 

109.26 

7 I 

109.79 

107.00 

109.79 

72 

110.19 

12C.  1 A 

110.19 

73 

1 10. AA 

107. 2H 

110.66 

7 ** 

111.13 

1 2C . AH 

111.13 

7 S 

1 1 1 .Al 

107. as 

111.61 

7 A 

1 12.09 

121.1S 

112.09 

77 

1 12. SB 

1 09. H8 

112.58 

78 

113.07 

1 1 A.  1 A 

113.07 

79 

| 18.H2 

1 IS. 99 

1 IB.H2 

80 

119.21 

1 18.09 

119.21 

B 1 

119.72 

1 18.65 

119.72 

82 

120.23 

1 19.03 

120.23 

B 3 

135,99 

13S.99 

1 35 ,99 

93 


KEVLhBtKAT I On  levels  at  TIMES  COKNESPONDING  TO  2-6Ar 
ACOUSTIC  TRAVEL  TIMES  BETwLEN  SOURCE  AND  TARGET. 

( 2 PINGS  --  TIME  d E T i,  E E N PINGS  IS  10. JO  SECONDS.! 
THE  SOUND  PROPAGATES  BETWEEN  SOURCE  ANU  TARGET  BT 

DIRECT  PATH 


target  t*o-aat  bottom 
RANGE  Tint  REVERB 

KTD  StC  DB 


SURFACE  VOLUME  TOTAL  PING  I PING  2 

REVERB  REVERB  REVERB  REVERB  REVERB 

Db  OB  DB  DB  OB 


.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
. 0 

1.23 
2 « R 3 
3.63 
0.03 

6*  OR 
7.2S 
0,05 
9 • 66 
10*67 

-300. 
-3R. 
-R  1 . 
■59  . 
-31  . 
■36. 
-67. 
-52. 

-91. 

• 0 

12*07 

-37. 

.0 

13.2B 

-36. 

.0 

10.08 

-37. 

.0 

15*69 

-38  . 

• 0 

16*90 

-RC 

.0 

18.10 

- R 2 • 

.0 

19.31 

-RS. 

.0 

20.51 

-R9  . 

• 0 

21.72 

-SR. 

.0 

22.92 

- 7 j • 

.0 

2R.I3 

-76. 

.0 

25.33 

-09. 

.0 

2 6 • 5 R 

-102. 

.0 

27. 7R 

-102. 

.0 

20. 9R 

-90. 

.0 

30.15 

-06  . 

.0 

3l  .35 

-82. 

.0 

32.5b 

-78. 

.0 

33*82 

-75. 

.0 

3 R . 9 5 

-73. 

. 0 

36.15 

-72. 

.0 

3 7 • 3 R 

-71  . 

.0 

30. 6R 

-71  . 

.0 

39.BR 

- 7 o . 

.0 

R 1 .UR 

-70. 

.0 

R 2 • 2 7 

-7o. 

.0 

R 3 • R 7 

-7  1 . 

• 0 

RR.60 

-71  . 

.0 

■Tma 

-71. 

.0 

-72. 

.7 
.0 
.8 
.3 
.9 
• C 
. 6 
. 1 
, R 

-300.0 
-300.0 
- 30  u • 0 
-89.  1 
-70.0 
-73.0 
-7R.8 
-78.9 
- 6 R . 0 

-36.2 
-21.5 
-22.6 
-23.9 
-2R.3 
-25.6 
-27.0 
■ 28.7 
-29.9 

-39 
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ACTIVE  PERFORMANCE  PREDICTION 

MEDITERRANEAN  AINTER  SOUND  SPEED  profile,  NUC  Tn  1006,  P,  23 


DIRECT  PATH 


OUTP 

1 .0  u 

ZX  900.0  U 

ZTG  1000. u 0 

FREU 

2.0  u 

L A A 

5.0  U 

V A I 

15.0  U 
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35.0  T 

MS 
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zbm 
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1.0  T 
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5.0  T 
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DELPR 

10.0  U 
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tav 
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5.99 
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27.0 
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ACTIVE  PERFORMANCE  PREDICTION 

MEOI  TERRAnEan  a 1 NTEK  SOUND  SPEED  PROFILE,  NUC  Tw  1006,  P.  23 
DIRECT  path 
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BE  A M X IS  THE  SIMXI/X  RESPONSE  FUNCTION  FOR  A CONTINUOUS  LINE  ARRAY 
BE  A MR  IS  THE  SInIXI/X  RESPONSE  FUNCTION  FOR  A CONTINUOUS  LINE  ARRAY 


The  LETTERS  FOLLOWING  The  NUMBERS  in  the  HEADING  on  the 
PERFORMANCE  PREDICTION  page  SIGNIFY  The  FOlLOaInG, 

F 

u user  has  specified  this  Input  value  on  this  run 

BLANK  set  BY  USER  ON  A PREVIOUS  RUN 

T TYPICAL  value  set  IN  ABSENCE  OF  USER  SPECIFICATION 

c CALCULATED  INTERNALLY  BY  PROGRAM  (USER  OPTION) 

A ARONG  value  SPECIFIED  BY  user  (OUT  OF  range  of  allo*. 

able  VALUES.  THE  PRINTED  VALUE  IS  A TYPICAL  VALUE.) 
n not  an  Input,  calculated  internally. 
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REVERBERATION  levels  at  times  corresponding  TO  2-WAY 
ACOUSTIC  TRAVEL  TIMES  BETWEEN  SOURCE  AND  TARGET, 

( 2 P1N&S  --  TIME  BETWEEN  PINGS  IS  10.00  SECONOS.I 
THE  sound  PROPAGATES  BETWEEN  SOURCE  AND  TARGET  BY 


surface  bounce  --  bounce  i 


target 

Two-way 
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reverb 

REVERB 

REVERB 
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DB 
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-50.6 
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-96.9 
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90.0 
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-SU.  1 
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-51.2 
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-61.9 
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50.71 
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5 | . 9 l 

-79.8 

-63.9 

-56. 1 

-52,1 

-55.3 

-59.8 

MM.O 
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•62.6 

-57.7 

-59,2 

-56.9 

-58.2 

MS.O 

SM.32 

-76*9 

-62.7 

-58.7 
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-57.2 

-59.9 

96.0 

55.53 

-77.9 

-69.3 

-59.7 

-56.2 

-58.9 

-60  • 2 

97.0 

56.73 

-78.9 

-65.9 

-60.6 

-55.9 

-59.9 

-57.6 

mB.O 

57. 9M 

-79.9 

-67.5 

-61.9 

-57.0 

-60.9 

-59.8 

M 9 • 0 

59.19 
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-68.5 

-62.  1 

-57.5 

-61.2 

-59.9 

S0.0 

60*39 

-8  | • 8 

-69.0 

-62.5 

-57,6 

-61.6 

-59.7 

SI  .0 

6 | . 6 1 

-83.2 

-69.5 

-52.2 

-51.9 

-52.  1 

-59.5 

S2.0 

62.75 

-89.6 

-69,9 

-57.6 

-55.2 

-57.3 

-59.3 

S3.0 

63.95 

-86.2 

-70.3 

-59.9 

-56.3 

-59. S 

-59.0 

59.0 
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-60.6 

-57,3 

-60.2 

-60.9 

5S.0 

66.36 
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-56.2 

-59.8 

-56.0 

-61.1 

56.0 

67*56 

-92.9 

-68.0 

-60.2 

-57.9 

-59.6 

-61.9 

57. 0 

68.76 

-95.9 

-67.3 

-60.7 

-58.6 

-59.9 

-69.6 

68.0 

69.97 

-99,9 

-66.7 

-63.8 

-58.8 

-59.8 

-65.5 

59.0 

7 | • 1 7 

-300.0 

-66.0 

-60.7 

-56.8 

-59.6 

—66 . 7 

60.0 
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62.0 
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63.0 
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77.16 
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-71.7 
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-65. 1 
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69.0 
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-68.5 
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-73.7 

-71,9 
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75.0 
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78.0 
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83.0 
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-300 . 0 
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-97.5 
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active  performance  FRED | CT  l ON 

MEDITERRANEAN  AJNTER  SOUND  SPEED  PROFILE,  HOC  Tn  1006,  P.  23 


SURFACE  bQONCE  --  BOUNCE  0 


DU  T P 
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Zx 
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FREO 
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s • a 

T 
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10.0 
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10.0 

hub 
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T 
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T 

lAT 

35.0 

DUMP 

.0 

7 

OUTR 
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mode 
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SL 
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10.0 

MPR 

2.0 

tbp 

IO«C 

Pulse 

, 1 O + uCT 
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C 

B 6 S 

10.0  T 

BAR 

10.0 
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T 

nqut 
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TVC 
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GOOP 
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TGS 

15.0 

T 

RON 

-2.0 

T 

clip 

•0  T 

SIGMA 

6 • J 

T 

pfa 

. 1 00- 02  T 

SP 

5. 

T 

00 

.0000 

ft 

Alfa 

. 1 265  N 

ZL 

• Ou 

N 

ZC  J 

0«  N 

PHIB  10.20 

N 

CX  *1968,  1 3 

ft 
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CL 

#00 

N 

CS  r959 ,97  N 

CB  5017.87 

ft 

TSUR 
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ft 

T AV 

59,9  N 

NOCZ 

2 

N 

I RAP 
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3 

ft 

range 

T I ME 

LOSS 

ANGLE 

SIGNAL  REVERe 

REV 

S/MN 

EXCESS 

PROB 

(KYO) 

(SECI 

1 DB  ) 

( DEG  ) 

( Db  ) 

(Db) 
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(OB) 

( DB  ) 

detect 

SEE  THE  DIRECT  PATH  OUTPUTS  (MODE»l)  FOR  FlHSWONE 
PERFORMANCE  PREuICT ions  for  this  SOUND-SPEED  PROFILE, 
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ACTIVE  performance  prediction 

MEDITERRANEAN  *INTER  SOUND  SPEED  PROFILE.  NUC  Tn  1006,  P.  23 


SURFACE  bounce  --  BOUNCE  I 


OUTP 

1.0 

Z X 

900*0 

ZTG  1000. 0 

FRED 

2.0 

L*A 

5.0 

V#I 

15.0 

Salt 

35.0  T 

HS 

5.0  C 

zbm 

6090* 

refls 

1.0  T 

2 ONE  S 

2.0 

Phox 

5.0  T 

PHDR 

5.0  T 

DElPX 

10.0 

DELPR 

10.0 

MuB 

-27.0  T 

Muv 

-99.0  T 

zc 

2000.  T 

LAT 

35.0 

DUMP 

.0  T 

OUTR 

t>0 

MODE 

3*0  U 

SL 

135.0 

delth 

10.0 

MPR 

2.0 

tbp 

10.0 

pulse 

. 10  + OOT 

01 

25.6  C 

b*s 

10.0  T 

b » R 

10.0 

N I N 

-95.0  T 

NOUT 

-60.6  C 

T VC 

12.0 

GOOP 

9.9  c 

tgs 

15.0  T 

RON 

-2.0  T 

Clip 

•0  T 

SIGMA 

6.0  T 

pfa 

. 1 00 "02  T 

SP 

5.  T 

00 

.0000  N 

alfa 

.1265  N 

ZL 

.00  N 

ZCJ 

0*  N 

PH  I B 1 

D.  20  N 

CX  9968.13  N 

CT  9973.07  N 

CL 

• 00  N 

CS  <*959.97  „ 

Cb  5097.87  N 

TSUR 

61.6  N 

Tav 

59.9  N 

NOCZ 

2 N 

I RAP 

0 N 

I AMOS 

3 N 

Range 

TIME 

LOSS 

ANGLE 

SIGNAL  RE  V E Kb 

REV 

S/Mn 

EXCESS 

PROB 

IKTD) 

(SEC) 

(OB) 

( deg  ) 

(DB) 

IDB) 

lim 

(OB) 

( DB  ) 

DETECT 

39.0 

97.0’ 

195.8 

3.36 

-95.8 

-96.9 

TES 

9.6 

9.6 

. 7679 

HO  • 0 

98.30 

197.6 

3.98 

-97.6 

-50*  1 

YES 

6 • U 

6 , 1 

.8919 

9 1 .0 

99.50 

199,2 

3.59 

-99.2 

-51.2 

YES 

5.3 

5,9 

.8105 

9 2.0 

SO. 71 

200.7 

3.71 

-50.7 

-53.9 

NO 

5.9 

5.7 

.8255 

9 3.0 

S|  .91 

202. 1 

3.83 

-52.1 

-52.  1 

YES 

3.0 

3.2 

.6993 

99.0 

53.12 

203.2 

3,95 

-53.2 

-59.2 

NO 

3.3 

3.9 

. 7328 

95.0 

59.32 

209.0 

9.05 

-59.0 

-55.3 

NO 

3.1 

3.9 

.7357 

96.0 

55.50 

209.9 

9.16 

-59.9 

-56.2 

NO 

2.7 

3,7 

. 7258 

9 7.0 

56.73 

205.7 

9.27 

-55.7 

-55.9 

NO 

1.5 

2.3 

.6951 

9 8.0 

57.99 

206.5 

9.38 

-56.5 

-57,0 

NO 

1.5 

2.7 

.6730 

99.0 

59.19 

207.3 

9.98 

-57.3 

-57.5 

NO 

. 9 

2.2 

.6916 

SO.O 

60.39 

207.7 

9.59 

-57.7 

-57.6 

no 

• 5 

1 . 9 

,6196 

SI  .0 

6)  .61 

200,  S 

3.39 

-50.5 

-51.9 

YES 

9.  1 

9.2 

. 7509 

S2.0 

62.75 

206.5 

9.85 

-56.5 

-55.2 

NO 

.5 

1 .3 

.5801 

53.0 

63.95 

206.5 

5.01 

-56.5 

-56.3 

NO 

1 • 1 

2.2 

.6365 

5 9.0 

65*16 

206.5 

5.18 

-56.5 

-57.3 

NO 

1*6 

2.9 

.6852 

SS.O 

66*36 

206.6 

5.37 

-56.6 

-59 . 8 

NO 

. 3 

1.0 

.5609 

5 6.0 

67. S6 

206.7 

5.57 

-56.7 

-57 . 9 

NO 

1.5 

2.8 

,6780 

57.0 

68.76 

206.8 

5.80 

-56.8 

-58.6 

NO 

1 . 9 

3.5 

.7135 

58.0 

69.97 

2 06.9 

6.05 

-56. 9 

-58.8 

NO 

i • B 

3.9 

.7)05 

59.0 

71.17 

207. 1 

6.32 

-57.1 

•58.8 

NO 

1,7 

3.3 

.7037 

60.0 

72.37 

207.3 

6.62 

-57.3 

-58.8 

N 0 

1.5 

3.1 

.6937 

61.0 

73.57 

207.5 

6.95 

-57.5 

-58.7 

N 0 

1.2 

2.8 

. 6797 

62.0 

79.76 

209.9 

7.26 

-59.9 

-59.9 

NO 

-.9 

.8 

, 5<*95 

63.0 

75.96 

219.2 

7.50 

-69.2 

-60,6 

NO 

-9.9 

-3.) 

.3053 

6 9.0 

77.16 

216.0 

7.70 

-66.J 

-60,7 

NO 

-6.6 

-9.8 

.2219 

65.0 

78.96 

216.7 

6.30 

-66.7 

-63,2 

NO 

-6.8 

-9.9 

.2169 

66.0 

79.66 

2)6.8 

6.72 

-66.8 

-69.5 

N 0 

-6.7 

-9.7 

. 2231 

67.0 

8Q.86 

217.7 

7.19 

-67.7 

-65. 1 

NO 

-7.6 

-5.6 

. 1782 

68.0 

81.99 

220.  1 

8.90 

-70.  1 

•66.7 

NO 

-9.8 

-7.8 

.1189 

69.0 

83.19 

221 .0 

8.57 

-71.0 

-67.9 

NO 

-10.7 

-8.7 

.0991 

70.0 

89.33 

221.9 

8.75 

-71.9 

-68.2 

NO 

-1  1.5 

-9.5 

.0793 

7 1.0 

85.52 

222.8 

8.92 

-72.8 

-69.9 

NO 

-12.9 

-10.9 

.0596 

72.0 

86.72 

223.8 

9.09 

-73.8 

-70. 1 

no 

-(3.3 

-11.3 

.0388 

73.0 

87.91 

229.7 

9.26 

-79.7 

-71.0 

no 

-19.2 

-12.2 

.0222 

79,0 

89.10 

2 25.7 

9.93 

-75.7 

-71.7 

N 0 

- 1 5 . 1 

-13.1 

.0186 

75.0 

9Q.29 

2 2 6.6 

9.60 

— 76,6 

-72.3 

NO 

-16.1 

-19,1 

.0150 

100 


ACTIVE  PErF  OHMAnCE  PREDICTION 

MEDITERRANEAN  AINTER  SOUND  SPEED  PROFILE,  NUC  Tn  1006,  P.  23 


surface  bounce 
UuTP  1.0 

; --  BOUNCE  1 

2x  960.0 

ZTG  1000. 0 

F R E 6 

2.0 

L A A 

5.0 

V*  I 

IS.O 

Salt 

35.0  T 

MS 

5 . u C 

zbm 

6090. 

REFLS 

1.0  T 

ZONES 

2.0 

PrtDX 

5. Cl  T 

PHQR 

5 • U T 

DELPX 

10.0 

DELPR 

I 0. 0 

MuB 

-27.0  T 

M U V 

-99.0  T 

zc 

2000.  T 

LAT 

35.0 

DUMP 

.0  T 

OUTR 

1 .0 

mode 

3.0  U 

SL 

135.0 

DELTH 

10.0 

MPR 

2.0 

TBP 

10.0 

pulse 

. 1G+3JT 

01 

25.6  C 

B*S 

13.0  T 

BAR 

10.0 

N I N 

-95.0  T 

NOUT 

-60.6  C 

TVC 

12.6 

GOOP 

9.9  C 

TG5 

15.0  T 

RON 

-2.0  T 

clip 

• 0 T 

SIGM* 

6 . L T 

PF  A 

. 100-02T 

SP 

5.  T 

SO 

.0000  N 

alfa 

.1265  N 

ZL 

.Cu  N 

ZC  J 

0.  N 

Pm  I B 1 

w • 20  N 

LX  8 9 6 8, 1 3 N 

Cr  8973.07  N 

CL 

.06  N 

CS  8959.97  N 

Cb  5097.87  N 

TSUR 

6 1.6  N 

Tav 

59.9  N 

NOCZ 

2 N 

I Rap 

0 N 

1 AMOS 

3 N 

Range 

TIME 

Eoss 

anule 

SIGNAL  reverb 

REV 

S/Mr. 

EXCESS 

PROB 

(KYD) 

(SEC) 

(OB  ) 

(DEG) 

( OB  I 

(08) 

L 1 M 

( 08  ) 

( 08  ) 

DETECT 

76.0 

91  .rB 

227.6 

9.77 

-77.6 

-72.9 

NO 

- 1 7 , u 

-15.0 

.0113 

77.0 

92.6? 

2 2 d . 6 

9,98 

-78.6 

-73.5 

i.O 

- 1 8 , u 

-16.0 

.0076 

78.0 

93.86 

229.5 

10.11 

-79.5 

-78.1 

N 0 

- 1 9,  u 

-17.0 

.0039 

79. C 

95.19 

237.6 

9.61 

-87.6 

-78.3 

N 0 

- 2 7 » u 

-25.0 

.0000 

BO  • 0 

96.38 

238.6 

9.79 

-88.6 

-79.6 

nO 

-28.0 

-26.0 

• OCOO 

o 1 • 0 

97.57 

239.6 

9.98 

-89.6 

-86.5 

nO 

-29.0 

-27.0 

.0000 

02. 0 

98.76 

280.7 

10.16 

-90.7 

-87.6 

NO 

-30.1 

-28.  1 

.0000 

0 3.0 

100*00 

272.0 

• 10.07 

• 122.J 

-97.5 

N 0 

-61.8 

-59.8 

.0000 
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INPUTS,  run  3 


a sspRo  ^8«!oArj,^;i'<jE50s;i|j2.n3DiPKOFlLE? N0C  TN  io°6* p*  ^ 

• loo  jl/I.’f  loo  gaa  IT*5  250  ^*s"s5300,'?j79.J%oJeJ;®  ‘“°  ,H89«* 

• 1200  M80.S  Uoo  H9I.7  U0  2 | ’o°.  5r  ?7<  * ' **  ° 0 ,<,;8-6  «I00  1<476 

i-.J8m}r8,S7  los-  15b6-52  ««!  7i*s.ii0-  ,7b0 

2a  20,  2 T 0 SCO 
PMdX«0,  PhDKmj 
00  TP  0 

m u o £ ■ i ,Run 


• calculations,  nun  3 * 


8EAMX  IS  THE  SlNUI/X  RESPONSE  FUNCTION 

bEAMR  IS  the  sin.xi/x  response  function 


For  a CONTINUOUS  line 
for  A CONTINUOUS  line 


array 

array 


PERFORMANCE  PREo  ^^^'pAOE^sfoNIFY^THE^FOLEOiiyNG* 


u 

blank 

T 

c 

A 


t5mcIl“S^E^F*” VALUE  0N  TH,S 
calColate^i^nternally *b y ^progr ^ US£k  ^fEcification 

a R ONU  VALUE  SPEC  IF icn  fl  PROGRAM  (USER  OPTION) 

a^e  mlSeSs:ec^;c?k;:tse"  ,°f  °p  aglo*. 

Rot  an  input,  calculated  internal!!  * K*L  v*lue-> 


SOUND  SPEED  PROFILE  AFTER  EARTH  CURVATURE  CORRECTION  AND  CURVE  FITTING 


SSP  SEGMENT  PARAMETERS 


LAYER 

depth 

SOUND  SPEED 

(JASA,  VOl 

. 41.  NO.  2,  PP 

. 419-438) 

number  <ftj 

1FT/SEC  1 

C A * • 2 

ka 

ZA 

1 

• 000 

4881. 824251 

.2214913*08 

.4066880-09 

-. 1317770*05 

2 

6S • S66 

4883.67561 7 

. 2385609*06 

-.  1383460-06 

• 1074963*03 

3 

9 2*310 

4864. 191610 

.2386202*08 

-.  1573519-07 

.2258650*03 

<1 

1 6 H • 198 

4884.730599 

.2385417*08 

.3503356-08 

-. 1 129618*03 

S 

910*060 

4886.4155)7 

.2388801*08 

-.734)937-08 

.6599851*03 

6 

995.890 

4887.052483 

.2388441*08 

-.3073871-07 

.5350405*03 

MAX 

535*090 

4887. 1679002 

.2388441*08 

-.3073871-07 

.5350405*03 

a 

596*801 

4887. 15751 | 

,3582575*08 

-.5880864-12 

-.9214970*06 

9 

625*626 

4887.018249 

.2388296*08 

-.2389725-06 

*6241 136*03 

10 

652*989 

4886,533793 

.2388232*08 

-.2760177-06 

. 6279857*03 

1 1 

680*292 

4886. 1 1QI77 

.2399920*08 

-.3723707-07 

. 290677e ■ 03 

1 2 

964.329 

4858.022883 

.2355957*08 

. 3506278-06 

• 103456  1 *r,« 

13 

984.849 

4855,924476 

.2357178*08 

.8723283-06 

.1004841*04 

|4 

998.522 

485S. 162433 

.3530054*08 

-.  1369661-09 

-.59271 39*05 

1 s 

1273.013 

4847.81 6046 

.2348977*08 

.6197013-07 

. 1 362057*04 

1 6 

1 3 2 7 * 7u 6 

4846,802134 

.2347903*08 

• 853  1 808-08 

. 1577096*04 

l7 

149( .709 

4845.666728 

.2347286*08 

. 1632508-08 

• 1937840*04 

M 1 N 

1937*840 

4844.6794312 

.2347286*08 

. 1632508-08 

• 1 937840*04 

|9 

2952*988 

4848.959949 

.2333293*08 

. 303d l 65-09 

-.2462254*04 

20 

328) .020 

485 1 ,682678 

.2348170*08 

. 1 273930-08 

. 1 9008 13*04 

21 

3609*212 

4854.825910 

.2334022*08 

.4814 184-09 

-.8843195*03 

22 

992) *878 

4870 . 859697 

.225598 1*08 

. 1485947-09 

- . 1-326C  1 2*05 

23 

5673.810 

4881 .516253 

,2595228*08 

-.  | 1 7574  1 -09 

.3320 168*05 

29 

6562*554 

4894 , 20:663 

.2252955*08 

. 1243471 -09 

-. 1530363*05 

25 

9844.916 

4944 , 954444 

.2169439*08 

.8038726-10 

-.2761407*05 

26 

1 1 7 2o  « j 3 9 

4977.532418 

.2959689*08 

-.9563049-10 

. 5682859*05 

27 

1 36 1 2 • u84 

Sul  1 . 1 651 79 

.3792102*08 

-.5496851-10 

. 1099431*06 

28 

17379.253 

5 j77 , 355440 

• GuQOOOG 

.0000000 

• 0003000 
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KEVERBERAT  1 0 i'J  LEVELS  AT  TIMES  CORRESPONDING  TO  2-,YAY 
ACOUSTIC  TRAVEL  TIMES  BETWEEN  SOURCE  AND  TARGET. 

( 2 PINGS  — TIME  BETaEEN  pings  IS  lc»00  seconds.) 
THE  SOUNO  PROPAGATES  BETWEEN  SOURCE  ANO  TARGET  BY 


Surface  Channel 


target 

T A 0 - u A Y 

bottom 

SURFACE 

VOLUME 

total 

PING  1 

PING  2 

range 

T I ME 

REVERb 

REVERB 

REVERB 

reverb 

REVERB 

REVERB 

KYD 

SEC 

DB 

DB 

Ob 

06 

OB 

08 

1.0 

1.23 

-300.0 

-5.0 

.7 

1.7 

1 . 7 

-30.5 

2.0 

2 . M 6 

-300.0 

-6.3 

-8.7 

-5.5 

-5.5 

-32. M 

3.0 

3.69 

-300.0 

-6.5 

-IM.3 

-5.8 

-5.8 

-3M.0 

H . 0 

M . 9 2 

- 300 ■ 0 

-25.3 

-1  8.  M 

-16.6 

-16,7 

-35. M 

5.0 

6>15 

- 3 uu  . 0 

-27.6 

-21.7 

-19.7 

-19,8 

-35.9 

6.0 

7*37 

-57.2 

-33.8 

-2M.M 

-22.6 

-22.9 

-3M.7 

7.0 

6 * 60 

-55.1 

-MO.  9 

-26.8 

-25,3 

-25,6 

-37.0 

e.o 

9.63 

-57. M 

-M6.2 

-29.3 

-27,8 

-28.0 

-M  1 .7 

9.0 

1 1 . 'j  6 

-58.6 

-M  8 , 5 

-31.6 

-30.0 

-30.2 

-MM.  1 

10.0 

12.29 

-7M.2 

-M7  . 1 

-33.8 

-31.9 

-32.1 

- M 5 . 5 

1 1 .0 

13.52 

-53.7 

-M5 . 6 

-35.8 

-33.6 

-33.8 

-M  6 . 6 

12.0 

1 M • 7 S 

-50. M 

-M3, a 

-37.7 

-35,0 

-35.2 

- M 7 . 6 

1 3.0 

IS. 98 

-52.8 

-MO. 7 

-39.5 

-35.6 

-35.9 

-M8.7 

1 MO 

17*21 

- 6 0 . 0 

-37.0 

-M  1 .2 

-3M.9 

-35.0 

-M9  • 7 

<5.0 

1 8 • MM 

— 8 u . 7 

"37.8 

-M2. 8 

-36.0 

-36.  1 

-50.8 

16.0 

19.66 

-82.6 

-M7.2 

-MM.M 

-MO, 9 

-Ml  .2 

-51.9 

1 7.0 

20*69 

-72. M 

-62 . 3 

- M 6 . 0 

-M3  . M 

-M3. 9 

-53.0 

ia.o 

22.12 

-65*1 

-70.3 

-M7.6 

-MM  . 7 

-MS.  3 

-53.9 

l 9 , 0 

23.35 

-61*5 

-70.0 

-M8. 7 

-M  5 . 6 

-M6.  M 

-5M.7 

20.0 

2m.  58 

-59.9 

-69.3 

-M9.9 

-M6.8 

-M7.5 

-55.  1 

21.0 

2S.B1 

-59. M 

-68.5 

-51.0 

-M7.6 

-MS. 5 

-5M.8 

22.0 

27. OM 

-59.9 

-67.0 

-52.2 

-M8 . 9 

- M 9 . 6 

-57.0 

23.0 

28.27 

-60*6 

-65.9 

-53. M 

-50.  1 

-50.7 

-59.3 

2 M • 0 

29.50 

-62.  1 

-65.  1 

-6M.5 

-51,3 

-51.8 

-61  ,0 

25.0 

30*73 

-63*6 

-6M.0 

-55.7 

-52. M 

-52.8 

-62. M 

26.0 

3 I • 9S 

-65.  M 

-63.1 

-56.8 

-53.  M 

-53.8 

-63.7 

27.0 

33*18 

-67.5 

-61.7 

-57.9 

-5M.2 

-5M.6 

-6M.8 

28.0 

3 M • M 1 

-69 . 8 

-60.0 

-59.0 

-5  M , 7 

-55,  l 

-65.9 

29.0 

3 S * 6 M 

-72. M 

-57.9 

-60.2 

-5  M , 6 

-5M.9 

-66,9 

30*0 

36  *67 

-75.5 

-60.2 

-61.3 

-56.3 

-56.6 

-67.9 

31.0 

38*10 

-79. M 

-65. 1 

-62. M 

-"5  8.6 

-59.0 

-68.9 

32.0 

39.33 

- 8 M . 6 

-69.9 

-63.5 

-60.3 

-60.8 

-69,8 

33.0 

Mo  .56 

-92.  1 

-7M.7 

-6M.6 

-6  1,7 

-62.3 

-70.8 

3 M • 0 

M|  .79 

-89.9 

-67.8 

-65.7 

-63.0 

— 63,6 

-71.7 

35.0 

M 3 • u2 

-89.3 

-66.9 

-66.8 

- 6 M , 0 

- 6 M , 7 

-72.7 

36.0 

MM  • 25 

-88.7 

-69.9 

-67.9 

-65. 1 

-65.7 

-73.6 

37.0 

M5.M7 

-88.0 

-70.8 

-68.9 

-66.  1 

-66.7 

-7M.6 

36.0 

M 6 • 7 0 

-88 . 6 

-71.8 

-70.0 

-67.1 

— 67,8 

-75.5 

39.0 

M7.93 

-9  1.0 

-72.6 

-71.1 

-68.1 

-68.7 

-76.5 

MO  • 0 

M 9 . 1 6 

-93.  1 

-73.5 

-72.  1 

-69.0 

-69.7 

-77. M 

Ml  .0 

50-39 

-95.2 

“ 7 M . 3 

-73.2 

-70.0 

-70.7 

-78.  M 

M 2 • 0 

Si  .62 

-98.0 

-75.0 

-7M.2 

-70.9 

-71.6 

-79. M 

M 3 • 0 

52*65 

-106.0 

-7S.9 

-75.3 

-28.7 

-72.5 

-28.7 

mm.o 

5m  « 08 

-112.0 

-76.7 

-76.3 

-33. M 

-73.5 

-33.  M 

M 5 . 0 

5S  • 3 1 

-118.9 

-77.5 

-77.3 

-33.5 

-7M.M 

• 33.5 

104 


-38.8 

1.2 

3.2 
6 • l 
8.8 
8.1 
8.0 
0.8 

8.3 
-80.8 
-82.8 
-85.6 

7.8 

0.  I 

1.9 

3.3 

8.7 

5.9 

7.1 

8.1 

9.0 

7.7 
-73.5 
-75.7 
-80.7 
-82.7 
-88.6 
-87.0 
-88.8 
-90.7 

-102.9 

-103.7 

-108.6 

-105.5 

-107.8 


-78.8 
-79 
-80 
-81 
-82 
-83 
-29 
-89 
-52 
-55  - 
-57.8 
-60  • 1 
-62.5 
-68.8 
-67. 1 
-69.3 
-71.3 
-73.  1 
-78.7 
-76.0 
-77.0 
-8  3 “ 


103.5 
108.8 

105.2 
106.1 

108.6 
109.8 

110.3 
111.2 
112.1 

113.1 

118.1 

115.2 

116.3 


-73.8 
-76.0 
-81.3 
-83.2 
-85.2 
-87.  I 
-89.0 
-91.0 

109.8 

110.8 

111.9 

112.9 
300.0 


SURFACE  CHANNEL 


TARGET 

T60-AAY 

BOTTOM 

SURFACE 

VOLUME 

TOTAL 

P I NG  1 

PING  2 

RANGE 

TIME 

reverb 

REVERB 

REVERB 

REVERB 

REVERB 

REVERB 

KYD 

SEC 

DB 

DB 

DB 

DB 

DB 

OB 

REVERBERATION  LEVELS  at  TIMES  CORRESPONDING  TO  2-rtA v 
ACOUSTIC  TRAVEL  TIMES  BETWEEN  SOURCE  AND  TARGET, 

( 2 PINGS  --  TIME  BETWEEN  PJNGS  IS  10.00  SECONDS.) 
THE  SOUND  PROPAGATES  BETWEEN  SOURCE  AND  TARGET  BY 


125 
125 
121 
117 
1 18 
1 I 1 
108 
106 
108 
102 
101 
100 
100 
-99 
-99 
-99 
-99 
-99 
-99 
-99 
100 
100 
100 
101 
101 
102 
102 
103 
1 08 
108 

105 

106 
106 

107 

108 


56.58 

57.76 

58.99 


78*66 
79.89 
81.12 
82.35 
83.57 
89.80 
86.03 
87.26 
88*89 
69.72 
90.95 
92*18 
93*81 
99.68 
9S>  86 
97.09 
98.32 


67*60 
68.83 
5 


active  performance  predict  ion 

..ESTERN  N,  ATLANTIC  i.INjER  SOUND  SPEED  PHOTlLE.  NUC  TU  lOu6,  P.  9 7 


surface  Cha,»i 
dutp  .0  u 

.EL 
Z A 

2C  • Z 0 

ZTG  3LL.U  U 

f Ktw 

2.3 

L A A 

5.0 

V A I 

15.0 

salt 

3 5 . 3 T 

MS 

5 • w C 

ZbM  17372.  u 

REFL5 

1.0  T 

ZONES 

2*0 

Phgx 

• J U 

PhDK 

• w U 

DEtPX 

10.0 

delpr 

10.0 

mob 

-27. C T 

Muv 

9 9 . L T 

ZC  2000.  T 

LAT 

35.3 

Dump 

.0  T 

OU  TR 

I *0 

mode 

1.3  U 

SL  135.L 

delth 

10.3 

MPR 

2.0 

TBP 

10*0 

pulse 

, 1 0 + C J T 

01  25.6  C 

b«S 

13.0  T 

BnR 

10.0 

N I N 

— 9 S . 0 1 

MOOT  - 

6 C • 6 C 

TVC  12.  j 

GDOP 

9.9  C 

TGS 

15.0  T 

RON 

-2.0  T 

Clip 

•0  T 

SIGMa 

6*0  T 

pfa 

. 100-02T 

SP 

5.  T 

00 

. 0 1 07  N 

Alfa  . 

1 283  N 

ZL  992. 13  N 

ZC  J 

6 15  2>  N 

P M I 6 15.93  N 

CX  9682.39  N 

CT  9885. S3  n 

CL  9 a 6 7 . 19  N 

CS  9661.8?  N 

CB  5077.36  N 

TSUR 

9 9.0  N 

T A V 

36. H N 

NOCZ 

0 N 

1 Rap 

0 N 

I AMOS 

1 N 

range 

T 1 ME 

LOSS 

angle 

signal 

REVERB 

RE  ? 

S/MH 

EXCESS 

PROB 

IKTD) 

(SEC  ) 

( Ob  ) 

(OECi) 

( DB  ) 

( 06  ) 

L I M 

(DB) 

( OB  ) 

DETECT 

1 • 0 

1.23 

12C.6 

• GG 

29.9 

1.7 

TES 

32.  1 

32.  1 

1 .0300 

2.0 

2.96 

133,2 

• CG 

16.8 

-5.5 

YES 

26.7 

26.7 

1 ,0300 

3.0 

3.6? 

190.6 

• w)0 

9.2 

-5.6 

YES 

l ? , 9 

19.9 

1 .0000 

9.0 

9 . 92 

196,9 

• Ou 

3.6 

-16,6 

YES 

2 9 • 7 

29.7 

1 .0000 

5.0 

6 . 1 S 

(53,9 

• GO 

-.9 

-19.7 

YES 

23.2 

23.2 

1 .0000 

6.0 

7.37 

159.6 

• GO 

-9.8 

-22.6 

YES 

22.3 

22.3 

1 .0000 

7.0 

8.60 

1S6.3 

• wG 

-8.3 

-25.3 

YES 

21.9 

21.9 

I .0000 

8.0 

9.83 

16  1,6 

• wO 

-11.6 

-27.8 

YES 

20.7 

20,7 

1 ,0300 

9.0 

l 1 .0* 

169.5 

• oC 

-19.5 

-30.3 

YES 

1 ? . 9 

19.9 

1 .0000 

10.0 

12.2? 

167.3 

• GG 

-17.3 

-3  1.9 

YES 

1?.  1 

1 9.  1 

1 . OuOO 

11.0 

13.52 

1 6?  . 8 

• GO 

-19.8 

— 3 3 * 6 

YES 

18,2 

16,2 

1 .0000 

12.0 

19. 7S 

1 72.3 

• JO 

-22.3 

-35.0 

YES 

17.1 

17.1 

. 9966 

1 3.C 

15.98 

179,6 

• GG 

— 29.6 

-35.6 

YES 

1 5 . 5 

15.5 

.9903 

19.0 

17.21 

176.8 

• GG 

— 2 6.6 

-39.9 

YES 

12.5 

12.5 

. 9789 

15.0 

18.99 

178.? 

.00 

-28.9 

-36.3 

YES 

1 1 . 5 

11.5 

.9652 

16.0 

19.66 

160.9 

• GO 

-30. 9 

-90*9 

YES 

19,3 

19.3 

.9857 

17.0 

20.8? 

182.6 

• 00 

-32.8 

-93.9 

YES 

19.8 

19.8 

.9879 

18.0 

22.12 

189.2 

• GO 

-39.2 

-99.7 

YES 

19.7 

19.7 

. 9872 

19.0 

23.3b 

165.7 

• GO 

-35.7 

-95.6 

YES 

19.2 

19.2 

.9655 

20.0 

29.56 

187.1 

• 00 

-37.  1 

-96.8 

YES 

13.7 

13.7 

.9839 

2 l .0 

25.81 

1 66.5 

• 00 

-38.5 

— 97,6 

YES 

1 3 . u 

13.0 

. 98U9 

22. 0 

27.09 

169.8 

• 0 0 

-39.8 

-98.9 

YES 

1 2 , 7 

12.7 

. 9799 

23.0 

28.27 

191.2 

• 00 

-91.2 

-50.  1 

YES 

12.9 

12.9 

. 97B7 

2 9.0 

2 9.5m 

192.5 

• 00 

-92.5 

-51.3 

YES 

1 2 • u 

1 2.  1 

. 9773 

2 5.0 

30.73 

193,9 

• Gw 

-93.9 

-52. 9 

YES 

11.9 

11.7 

.9690 

26.0 

3|  .95 

195.2 

• 0 c 

-95.2 

-53.9 

N 0 

1 0 . 6 

11.2 

,9589 

27.0 

33.10 

196.5 

• GO 

-96.5 

-59.2 

NO 

ID.  1 

10.6 

.9999 

28.0 

39.9  1 

197.7 

• 00 

-97.7 

-59.7 

N 0 

?.  1 

9.7 

.9299 

29.0 

35.69 

199.0 

• GO 

— 9 9 . j 

-59.6 

NO 

7.6 

8 . 9 

.899  1 

30.0 

36.87 

203.3 

*uO 

-50.3 

-56.3 

?„0 

7.9 

6.9 

,8956 

3]  .0 

38.10 

231.5 

• 0 G 

-5  1.5 

-58.6 

NO 

7.1 

8.7 

.9018 

32.0 

39.33 

2:2.6 

• 00 

-52. 6 

-60*3 

NO 

6.9 

8.3 

. 8?  1 9 

33.0 

90.56 

2:9.: 

• 00 

-59,3 

— 61.7 

nO 

5, 6 

7.5 

.8739 

39.0 

9 1.7? 

*05.2 

• <jCi 

-55.2 

- 6 3 • J 

NO 

9.6 

6 . 6 

.8526 

35.0 

93.02 

2 C 6 , 9 

• 0 J 

“56,9 

■69.!, 

nO 

3,5 

5.5 

.8157 

3 6.0 

99.25 

2m  7.7 

• GO 

-57.7 

-65.  1 

■*0 

2.9 

9.9 

. 7620 

37.0 

95.97 

*08 , 9 

• GO 

-58.9 

-66.1 

>0 

1.3 

3.3 

. 707  1 

ACTIVE  PERFORMANCE  Prediction 

aESTERN  N.  ATLANTIC  A I N T £ R SOUND  SPEED  PROFILE.  NOC  Tn  1006,  P.  97 


Surface  c h a 
OuTP  «G  u 

nel 

Za 

20.  Z U 

2 T G 

3GC.0  D 

FREU 

2.0 

L*  A 

5.0 

V A I 

15.0 

Salt 

35.0  T 

mS 

5.0  C 

ZBM  17372*  U 

kefls 

1.0  T 

ZONES 

2.0 

PhOX 

•0  u 

PhOR 

.0  0 

UElpX 

10.0 

DELPR 

10.0 

Mua 

-27.0  T 

MuV 

99.j  T 

zc 

2000.  T 

lat 

35.0 

DUMP 

.0  T 

OUTR 

1 *0 

MODE 

10  U 

SL 

135.0 

delth 

10.0 

MPR 

2.0 

TaP 

10. 0 

pulse 

. 10*OOT 

D I 

25.6  C 

6 6 S 

lu.O  T 

BnR 

10.0 

N I N 

•95.0  T 

nout  - 

60.6  C 

TVC 

12.0 

gdop 

9.9  c 

TGS 

IS.O  T 

RON 

-2.0  T 

clip 

•C  T 

SIGMA 

6.0  7 

hf  a 

. I00-U2T 

SP 

5.  T 

00 

. C l 0 7 N 

Alfa  . 

1283  N 

ZL  992.13  N 

ZC  J 

6 1 5 2 • N 

PHIS  15.93  N 

cx  9882.39  N 

CT  9885.53  N 

CL  9687.19  N 

CS  9881.89  N 

CB  5077.36  N 

TSOR 

99.0  N 

TAV 

36.9  N 

NOC  Z 

u N 

1 R AP 

0 N 

1 AMOS 

1 N 

range 

T I ME 

LOSS 

ANGLE 

SIGNAL  REVERb 

REV 

s/mn 

EXCESS 

PRUB 

(KYD) 

(SEC) 

( 06  ) 

( 0 E 6 ) 

( DB  ) 

(DB) 

L 1 M 

(08) 

( OB  ) 

DETECT 

38.0 

96.70 

2 1 j . 1 

.00 

-60. 1 

-67. 1 

NO 

• 2 

2.2 

.6912 

39.0 

97.93 

211.3 

.00 

-61.3 

-68.  1 

NO 

-.9 

1 . 1 

.5698 

9 0.0 

99  . | 6 

2(2.9 

.30 

-62.9 

- 6 9 . u 

NO 

-2.0 

-.0 

.9981 

Ml  .0 

SO. 39 

213.6 

• 00 

— 63.6 

-70.0 

N 0 

-3.2 

-1.2 

.9259 

9 2.0 

51.62 

219.8 

.00 

-69.8 

-70 . 9 

N 0 

-9.3 

-2.3 

. 3536 

*13.0 

52.  BS 

2|e.u 

• 00 

-66.0 

-26.7 

YES 

-32.9 

-32.9 

.0000 

99.0 

59.06 

217.1 

.00 

— 67.| 

-33.9 

YES 

-29.3 

-29.3 

.0000 

*15.0 

55.31 

218.3 

.00 

-68.3 

-33.5 

YES 

-30.9 

-3U  • 9 

. 0000 

96.0 

56.59 

2 19.5 

« u C 

-69.5 

-38.8 

YES 

-26.3 

-26.3 

.0000 

9 7.0 

57.76 

220.6 

. oc 

-70*6 

-91.2 

YES 

-25.  1 

-25.1 

• GOOD 

9 8.0 

58.99 

22  1.8 

. 30 

-71.8 

-93.2 

YES 

-29.3 

-29.3 

.0000 

99.0 

60.22 

222.9 

• Co 

-72.9 

-96. 1 

YES 

-22.8 

-22.8 

.0000 

SO  • 0 

6 1.95 

229. 1 

.00 

-79.1 

-98.9 

YES 

-21.8 

-21.8 

.0000 

S|  .0 

62.68 

225.2 

• 00 

-75.2 

-26.  1 

YES 

-99,7 

-99.7 

.0000 

S2.0 

63.91 

226.3 

.00 

— 76.3 

-28.9 

YES 

-93.9 

-93.9 

.0000 

S 3 • 0 

65.19 

227.5 

.00 

-77.5 

-30.9 

YES 

-92.7 

-92.7 

.0000 

59.0 

6 6*37 

228.6 

.00 

-78.6 

-38.3 

YES 

-35.9 

-35.9 

.0000 

ss.o 

67.60 

229,7 

.00 

-79.7 

-90.8 

YES 

-39.6 

-39.6 

.0000 

S 6 * 0 

68.83 

230.8 

.00 

-80.8 

-92.8 

YES 

-33.8 

-33.8 

• 0000 

S7.G 

70.05 

232.0 

.00 

— 8 2 . u 

-95.6 

YES 

-32.3 

-32.3 

• 0000 

SS.O 

7|  .28 

233. 1 

.00 

-83.1 

-97.8 

YES 

-31.9 

-31.9 

.0000 

S 9 • 0 

72.51 

239.2 

.00 

-89.2 

-50. 1 

YES 

-30.6 

-30.6 

.0000 

60.0 

73.79 

235.3 

.00 

•85.3 

-51.9 

YES 

-30.9 

-30.2 

.0000 

61.0 

79.97 

236.9 

• oc 

• 86.9 

-53.3 

N 0 

-30.5 

-30.  1 

.0000 

62.0 

76.20 

237.5 

.00 

-87.5 

-59.7 

NO 

-30.7 

-30.1 

.0000 

63.0 

77.93 

236.6 

.00 

-88.6 

-55.9 

NO 

-31.1 

-30.2 

.0000 

69.0 

78.66 

239.7 

.00 

-89.7 

-57. 1 

NO 

-31.7 

-30.5 

.0000 

65.0 

79.89 

290 . 8 

.00 

-90.8 

-58.  1 

NO 

-32.9 

-30.9 

.0000 

6&.C 

81.12 

291.9 

.00 

-91.9 

-59.0 

NO 

-33.  1 

-31.5 

.0000 

67.0 

82.35 

293.0 

• 00 

-93.0 

-67.7 

NO 

-32.7 

-30.7 

.0000 

68.  C 

83.57 

299  . 1 

• OG 

-99.1 

-73.5 

NO 

-33.6 

-31.6 

.0000 

69.0 

89.80 

295.2 

.00 

-95.2 

-75.7 

NO 

-39.6 

-32.6 

.0000 

70.0 

86.03 

296.3 

• OG 

-96.3 

-80.7 

N 0 

-35,7 

-33.7 

.0000 

7 | .0 

87-26 

297.9 

.00 

-97.9 

-82.7 

NO 

-36.6 

-39.8 

.0000 

72.0 

68.99 

296.5 

.00 

-98.5 

-89.6 

N 0 

-37.9 

-35.9 

.0000 

73.0 

89.72 

299.6 

.00 

-99.6 

-67.9 

N 0 

-39. U 

-37.0 

.0000 

79.0 

9u  . 95 

250.7 

.00  — 

1 GO  . 7 

-68.8 

NO 

-90.  D 

-38.0 

.0000 

KL‘ 
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active  performance  PREDICTION 

AESTERN  N,  ATLANTIC  AINTER  SOUND  SPEED  PROFILE*  NUC  TN  |006,  P,  97 


n 


surface  channel 
UuTP  *0  U 2x 

20.0  U 

ZTG 

300.0  U 

FREH1 

2.0 

LA  A 

5.0 

V A 1 

15*0 

Salt 

35.0  T 

HS 

5.0  C 

ZBM  17372.  U 

refls 

1.0  T 

Z ONES 

2.0 

PhDX 

• 0 U 

PHDR 

.0  u 

DELPX 

10.0 

delpr 

10.0 

MuB 

-27.0  T 

M u V - 

99.0  T 

zc 

2000.  T 

LAT 

35.0 

DUMP 

.0  T 

OUTR 

1 >0 

MODE 

1*0  U 

SL 

135.0 

dEltm 

10.0 

MPR 

2.0 

TBP 

10.0 

Rulse 

. 10*OOT 

0 I 

25.6  C 

BAS 

10.0  T 

BaR 

10.0 

MN 

-95.0  T 

nout  - 

60 . 6 C 

TVC 

12.0 

GDOP 

9.9  C 

tgs 

15.0  T 

RON 

-2.0  T 

Clip 

•0  T 

SIGMA 

6.0  T 

pfa 

. I00-02T 

SP 

5.  T 

00 

.0107  N 

Alfa  . 

1283  N 

ZL  992.13  N 

zcj 

6152*  N 

PH1B  15.93  N 

CX  9882.39  N 

CT  9885. S3  N 

CL  9887.19  N 

CS  9881.89  N 

CB  5077.36  N 

TSUR 

99.0  N 

Tav 

36.9  N 

NOCZ 

0 N 

IRAP 

0 N 

I AMOS 

1 N 

RANGE 

T I ME 

loss 

ANGLE 

SIGNAL  reverb 

REV 

S/MN 

EXCESS 

PROB 

(KYO) 

(SEC) 

(DB  ) 

(DEG) 

(DB  ) 

(DB) 

LI" 

(DB) 

(DB) 

OETECT 

75*0 

92.  IB 

251.8 

.00 

-101.8 

-90*7 

NO 

-91  . t 

-39.  | 

.0000 

76*0 

93.9) 

252.8 

.00 

-102.8 

-102.9 

NO 

-92.2 

-9U.2 

.0000 

77.0 

99.69 

253.9 

.00 

-103.9 

-103.7 

N 0 

-93.3 

-91.3 

.0000 

78.0 

95.86 

255.0 

.00 

-105.0 

-109.6 

NO 

-99.9 

-92.9 

,0000 

79.0 

97.09 

256.  1 

.00 

-106. 1 

-105.5 

NO 

-95,9 

-93.9 

.0000 

80.0 

98.32 

257 . 1 

.00 

-107.1 

-107.8 

NO 

-96.5 

-99. S 

.0000 

i 
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• INPUTS,  RUN  <1  # 


PMDX-RO,  PhDR»HU 
M0DE»2,Run 


• calculations,  run  r • 


beamx  is  t rl  sinixi/x  Response  function  for  a continuous  line  arrat 

BE  A MR  IS  TMt  S l N ( A I / X RESPONSE  FUNCTION  FOR  A CONTINUOUS  LINE  ARRAY 


the  letters  following  the  numbers  in  the  heading  on  the 
PERFORMANCE  prediction  PAGE  SIGNIFY  The  FOLLOWING, 


u user  has  specified  this  input  value  on  this  run 
blank  set  oY  UBER  On  a previous  RUN 

T TYPICAL  value  set  IN  ABSENCE  OF  USER  SPECIFICATION 

C CALCULATED  INTERNALLY  BY  PROGRAM  (USER  OPTION) 

ti  AMONG  VALUE  SPECIFIED  BY  USER  (OUT  OF  RANGE  OF  ALLOW- 

ABLE values,  the  printed  value  is  a TYPICAL  VALUE,) 
n not  an  Input,  calculated  internally. 
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REVERBERATION  levels  AT  times  CORRESPONDING  TO  2-rAY 
ACOUSTIC  TRAVEL  TIMES  BETWEEN  SOURCE  AND  TARGET, 

( 2 PINGS  --  TIME  oET„£EN  PINGS  IS  IC.uL  SECONDS.) 
the  sound  propagates  between  source  anu  target  by 


bottom  bounce  --  BOUNCE  | 


TARGET 

T*0-wAY 

BOTTOM 

surface 

VOLUME 

TOTAL 

PING  1 

PING  2 

RANGE 

TIME 

REVERB 

REVERB 

REVERB 

REVERB 

REVERB 

REVERB 

KYD 

SEC 

DB 

OB 

DB 

DB 

DB 

DB 

tf 

1.0 

IR.  1 1 

-52. R 

-R  9 . R 

-85.2 

-31.7 

— M 7 . 6 

-31.9 

2.0 

1 R . 30 

-62. 1 

-R7.3 

-83.6 

-32.7 

-R7 , 1 

-32.9 

f 

3.0 

IR.S2 

-71 .8 

-RR.9 

-83.0 

-33.7 

-MR. 9 

-3R.0 

R • 0 

1M.82 

-77.5 

-R2.5 

-82. 1 

-3R.8 

-M2. 5 

-35.6 

S.O 

15.31 

-76.6 

-R3.7 

-80.3 

-37.1 

-M3. 7 

-38.2 

1 

6.0 

15*80 

-76.1 

-RR.9 

-78.2 

-39,3 

-RR.8 

-R0.8 

7.0 

16.36 

-73.0 

-R6.B 

-80.  C 

-R2.G 

-M6.8 

-R3.7 

8.0 

17.07 

-R9.8 

-53.6 

-82.  1 

-RR.9 

- M 8 , 3 

-R7.5 

9.0 

17-77 

-R7.0 

-60. R 

-8R.5 

-R6.7 

-R6.8 

— 6 1,2 

10.0 

18.52 

-R6 . 7 

-6R  • 8 

-75.8 

— R 6 . 7 

-M6 . 7 

-73.6 

1 1 .0 

19.36 

- R 7 . 6 

-52.2 

-62. R 

-R6.2 

- R 6 . 2 

-66.6 

I 2 . 0 

20.21 

-R9.B 

-39.7 

-R8.9 

-38.8 

-38.8 

-58.3 

13.0 

2 1 . u 8 

-52.  R 

-28.0 

-R6.3 

-27.9 

-26,0 

-R9.6 

I R • 0 

22.03 

-55.8 

-29.1 

-R6.5 

-28.9 

-29.0 

-R5 . 7 

1S.0 

22*97 

-59.7 

-3U  • 1 

-R7.  1 

-29.9 

-30.0 

-R6.R 

16.0 

23*93 

-63.8 

-31.2 

-50 . 1 

-31.0 

-31.1 

-R7.0 

17.0 

2 M • 9R 

-68.5 

-36.3 

- 5 R . R 

-35.9 

-36.2 

-R  7 . 7 

1 8.0 

25.95 

-73. R 

-R  1 . 7 

-59. R 

-R0.7 

-R  1 . 6 

-R8. 3 

19.0 

26*97 

-78. R 

-R7.0 

- 6 B . 7 

-R5.7 

-M  7 . 0 

-51.8 

1 

20.0 

28  ■ u3 

-79.3 

-66.0 

-7R.9 

-55.0 

-66.9 

-55.3 

21.0 

29.09 

-69.5 

-93.2 

-96.2 

-5B.5 

-69.5 

-58.6 

22.0 

30.27 

-57.7 

-97.6 

-10R.3 

-56.5 

-57.7 

-62.6 

23.0 

3)  .23 

-R8. 1 

-96.8 

- 1 G R • 8 

-R  8 . 1 

-R  8 , 1 

-68.2 

2R.0 

32.32 

-R5 . 9 

-95.9 

-99.2 

-R5. 9 

— M 5 . 9 

-73.0 

25. 0 

3 3 * R 2 

— R 6 . 7 

-83.9 

-80.5 

-R6 , 7 

- M 6 , 7 

-75. 1 

26.0 

3m.  52 

-R7.R 

-67.8 

-80.0 

- R 7 4 R 

- M 7 . R 

-77.2 

27.0 

35*  6R 

-R8 . 2 

-67.6 

-79.2 

-R8. 1 

-R8.  1 

-79.6 

28.0 

36.75 

-51.2 

-67. R 

-78.8 

-51.0 

-5  1.0 

-82. 1 

l 

29.0 

37.67 

-55.1 

-67.2 

-79.7 

-5R.8 

-5R.8 

-8R.7 

30.0 

39*01 

-59.0 

-68.6 

-80.7 

-58. 6 

-58.6 

-87.5 

3 1.0 

R 0 • 1 R 

-63.0 

-70.  I 

-82.0 

-62.2 

-62.2 

-90. 1 

32.0 

R 1 .26 

-72.0 

-71.5 

-83.7 

-68.5 

— 68,6 

•86.  1 

> 

33.0 

R 2 • R 2 

-90>3 

-73.6 

-85.6 

-72.9 

-73.2 

- 8 R . 7 

3R.0 

R 3 • 5 7 

-108.6 

-75.7 

-87.5 

-7R.9 

-75. R 

— 8R  • 8 

35,0 

RR.  7 1 

-126.9 

-77.9 

-89.8 

-76.8 

-77.6 

-8R.  7 

36.0 

R 5 • 8 7 

-1R5.3 

-80. R 

-92.2 

-78,8 

-80.  1 

-8R.7 

1 « 

37.0 

R7.U2 

-137.3 

-83.0 

— 9 M . 7 

-80.6 

-82.7 

-8R.7 

38.0 

R8. 18 

-123.6 

-85.7 

-97.5 

-82.6 

-65. R 

-85.8 

39.0 

R 9 • 3R 

-109.9 

-86.7 

-100.3 

- 8 R « 5 

-88.  R 

-86.8 

if 

MO.O 

50.50 

-96.2 

-91.7 

-103,3 

-85.8 

-90.2 

-87.8 

Ml  .0 

5).66 

-85.0 

-9M.9 

-106. 5 

• 83.2 

-6R.6 

-88.8 

•? 

4 

R 2 • 0 

52.83 

-8R  . 9 

-98. R 

-109.7 

-28. R 

-8R.7 

-28  . R 

M 3 • 0 

53.99 

-8R.8 

-101.8 

-113.2 

-3R.7 

- 8 R . 7 

- 3 R . 7 

mr.o 

55.16 

-8R  . 7 

-105.5 

-1  16. R 

-31.2 

- 8 R . 7 

-31.2 

M 5 • 0 

56.33 

-8R.6 

-109.1 

-119.6 

-38.5 

-8R.6 

-38.5 

110 


1 


! 


REVERHLrAT I ON  LEVELS  at  TIMES  CORRESPONDING  TO  2-*aY 
ACOOSTlC  TRAVEL  times  6ET6EEN  SOURCE  and  target. 

( 2 RINGS  --  TIME  b E T « E E N RINGS  IS  1C. GO  SECOnDS.I 
the  SOUnO  PROPAGATES  BETWEEN  SOURCE  ANO  TAH6ET  bY 


BOTTOM  bounce  -«  BOUNCE  I 


T ARGET 

T.VO-..AY 

BOTTOM 

SURE  ACE 

VOLUME 

TOTAL 

PING  | 

PING  2 

RANGE 

TIME 

REVERb 

REVERB 

REVERB 

reverb 

REVERb 

REVERb 

KYD 

SEC 

OB 

Ob 

Db 

06 

Ob 

DB 

*6.0 

S 7 * S 1 

•65.2 

-112.9 

-122.1 

-90.9 

-b5 . 2 

-*0.9 

*7.0 

5 6 • 6 d 

• 86*2 

-115.2 

-129.0 

-92.9 

-66.2 

-92.9 

*6.0 

59*06 

• 87.2 

■117.1 

-125.6 

-95.6 

-87.2 

-95.6 

*9.0 

6 1 . 03 

•8B  . 3 

-1  18.7 

-126.8 

-97,9 

-68.3 

-97.9 

50.0 

62.21 

-89.3 

-120.3 

-126.  C 

-99.8 

-89,3 

-99.8 

Si  .0 

63.36 

-90*8 

-121.3 

-33.0 

-33.  u 

-33.0 

-51.6 

S2.0 

6*  • 56 

-92.9 

-81.8 

-36.0 

-3S.9 

-36.0 

-53.0 

53.0 

65.7* 

-93*9 

-80  * - 

-36.9 

-36.8 

-36.9 

-59.3 

5 9.0 

66.92 

• 95.9 

-78.9 

-39.8 

-39.7 

-39.8 

-55.5 

SS.O 

68*09 

-97.1 

-79.0 

-*1.9 

© 

• 

T 

• 

-91.9 

-56.6 

56. 0 

69.27 

-99.0 

-79.1 

-* 9 . 2 

— 99,0 

-99.2 

-57.7 

57.0 

70  * *5 

-100. B 

-82.5 

-96.8 

-96.5 

-96.8 

-58.6 

58.0 

7 | . 63 

-102.7 

-89.5 

-*8.9 

-98.5 

-98.9 

-59.5 

59.0 

72.01 

-10  9.6 

-86.3 

— 50.8 

-5U.8 

-5u.  8 

-72.9 

60.0 

73.99 

-106.8 

-88. 1 

-52.3 

-52.3 

-52.3 

-79.6 

61.0 

75.17 

-10  6.9 

-89,6 

-53.7 

-53.6 

-53.7 

-76.7 

62.0 

76*35 

-111*1 

-91.7 

-5*.9 

-59,9 

-59.9 

-8  1.6 

active  PERFORMANCE  PREDICTION 

aESTERN  N.  ATLANTIC  aINTER  SOUND  SPEED  PROFILE.  NUC  Tn  1006,  P,  97 


BOTTOM  BOUNCE  --  BOUNCE  l 


UUTP 

• 0 

Zx 

2C  • C 

ZTG 

300.0 

FREQ 

2.0 

LA  A 

5.0 

V 6 I 

15.0 

Salt 

35. C T 

hS 

S.u  C 

Z BM  1 7372* 

KEFLS 

1.0  T 

ZONES 

2.0 

PhOX 

90*  C U 

PHDR 

90. w 0 

dElp* 

10.0 

OELPR 

10.0 

MuB 

-27.0  T 

M U V 

99.0  T 

zc 

20  00  • T 

LAT 

35.0 

DUMP 

.0  T 

UUTR 

1«0 

mode 

2*0  U 

5L 

135.0 

DELTH 

10.0 

MPR 

2.0 

TBP 

10.0 

pulse 

. |0*00T 

D 1 

25.6  C 

BmS 

10.0  T 

BAR 

10.0 

N I N 

-95.0  T 

N 0 u t - 

6 C . 6 C 

TVC 

12.0 

gdop 

9.9  C 

TGS 

15.0  T 

RDN 

-2.0  T 

clip 

•0  T 

SIGMA 

6 . u T 

PFA 

. 100-02T 

SP 

5.  T 

00 

.0107  N 

alfa  . 

1277  N 

ZL  992.13  N 

zc  0 

6152*  N 

PHIS  15.93  N 

CX  9882,39  N 

CT  9685.58  N 

CL  9887.19  N 

CS  9881.89  N 

CB  5077.36  I, 

TSUR 

99.0  N 

Tav 

36.9  N 

NOCZ 

0 N 

I RAP 

0 N 

I AMOS 

I N 

RANGE 

T 1 ME 

LOSS 

ANGLE 

signal  reverb 

REV 

SVMn 

EXCESS 

PROB 

1 K Y 0 ) 

( SEC  i 

( oa ) 

(DEG) 

( OB  ) 

( Db  ) 

L IM 

(OB) 

( DB  I 

DETECT 

1 >0 

19.11 

223,7 

85,20 

-73.7 

-31.7 

YES 

— 37,6 

-37.6 

.0000 

2.0 

19.30 

225.1 

80,93 

-75.1 

-32.7 

YES 

-37.9 

-37,9 

. 0000 

3.0 

19.52 

07.6 

-75.67 

— 87,6 

-33.7 

YES 

-99.5 

-99,5 

.0000 

9.0 

19.82 

2 1 9 , D 

71.16 

-69.& 

-39.8 

YES 

-29.6 

-29.8 

.0000 

s.o 

15.31 

216,8 

6 7,06 

-68.8 

-37.  1 

yes 

-27.3 

-27.3 

.0000 

6*0 

15*80 

269, 1 

62.97 

-119.1 

-39.3 

YES 

-70.9 

-70.9 

.0000 

7.0 

16*36 

213.0 

59.15 

-63 . 0 

-92.0 

YES 

-16.7 

-16.7 

.0099 

8*0 

17*07 

208,9 

55.89 

-58.9 

-99.9 

YES 

-9.9 

-9.9 

.0711 

9*0 

17.77 

2 2 3,0 

52,63 

-73. C 

-96.7 

YES 

-22.9 

-22.9 

.0000 

10*0 

18.52 

213.9 

99.62 

-63.9 

-96.7 

yes 

-13,3 

-13.3 

.0182 

11*0 

19.36 

197,0 

97.12 

-97,0 

-96.2 

YES 

3.3 

3.3 

. 7039 

1 2.C 

20.21 

169,1 

99.61 

-39.1 

-38.8 

YES 

9.U 

9.0 

.7911 

1 3.C 

21.08 

165.7 

92.28 

-35.7 

-27.9 

YES 

-3.3 

-3.3 

. 2990 

19.0 

22.03 

185.1 

90.37 

-35.1 

-28.9 

YES 

-1.8 

— 1,8 

.3863 

1S.0 

22.97 

186.2 

38,97 

-36.2 

-29.9 

YES 

-1.8 

-1.8 

.3893 

16*0 

23.93 

188,8 

36.69 

-38.8 

-31.0 

YES 

-3 . 9 

-3.9 

.2903 

17*0 

29.99 

192,9 

35.19 

-92,9 

-35.9 

YES 

-2.1 

-2.1 

. 3686 

18*0 

25.95 

197.9 

33.79 

-97.9 

-90*  7 

YES 

-2.3 

-2.3 

. 3531 

19.0 

26.97 

209.3 

32.30 

-59.3 

-95.7 

YES 

-9.5 

-9.5 

.2333 

20.0 

28. 03 

212.2 

31,17 

-62.2 

-55.  a 

NO 

-5.2 

-9.5 

.2326 

21.0 

29.09 

229.2 

30.06 

-79.2 

-58.5 

NO 

- 1 5 , 5 

-19,0 

.0159 

22.0 

30.27 

298,2 

29.39 

-98.2 

-56.5 

no 

-90.5 

-39.9 

.0000 

23.0 

31.J3 

292,  1 

28.01 

-92.  1 

-98.1 

YES 

-90.2 

-90.2 

.0000 

29.0 

32.32 

229.6 

27.15 

-79.6 

-95.9 

YES 

-29.7 

-29.7 

.0000 

2 S • 0 

33.92 

229.3 

26.26 

-79.3 

-96.7 

YES 

-23.7 

-23.7 

.0000 

26.0 

39.52 

222.0 

25.50 

-72. C 

-97.9 

YES 

-2  0.7 

-20.  7 

.0000 

27.0 

35.69 

22  1.2 

29.82 

-71.2 

-98. 1 

YES 

-19.3 

-19.3 

.0000 

28.0 

36.75 

221.2 

29.19 

-71.2 

-5  1.0 

YES 

-16,8 

-16.8 

.0098 

29.0 

37.87 

221.8 

23.99 

-71.8 

-59.8 

no 

-19.9 

-19.2 

.0199 

30.0 

39.01 

222.8 

22.95 

-72.8 

-58.6 

no 

-1*1. 2 

— 12,6 

.0208 

31.0 

90.19 

229.2 

22.90 

-79.2 

-62.2 

NO 

-19.5 

— 12.6 

.0208 

32.0 

9|  .28 

225.9 

21.87 

-75.9 

-68.5 

nO 

-15.5 

-13.5 

.0179 

33.0 

92.92 

227,7 

21.93 

-77.7 

-72.9 

no 

-17.2 

-15,2 

.0109 

39.0 

93.57 

229.8 

20.99 

-79.8 

-79.9 

NO 

-19.2 

-17.2 

.0030 

35  . u 

“9.71 

232.1 

20.56 

-82.1 

-76.8 

NO 

-21.5 

-19.5 

.0000 

36.0 

95.87 

239.5 

20.20 

-89.5 

•78.8 

NO 

-23.9 

-21.9 

.0000 

37.0 

97.02 

237.1 

19.a3 

-87.1 

-80*6 

nO 

-26.9 

-29.9 

.0000 

ACTIVE  PERFORMANCE  prediction 

AESTERN  N,  ATLANTIC  AINTER  SOUND  SPEED  PROFILE.  NUC  TN  |006.  P.  H7 


IL. 


BOTTOM  BOUNCE  --  BOUNCE  | 


OuTP 

• 0 

ZX  20.0 

ZTG 

300.0 

freb 

2.0 

LAA 

5.0 

V*  I 

15.0 

Salt  35. o t 

MS 

5.0  C 

ZBM  17372. 

REFLS 

1.0  T 

ZONES 

2.0 

PhDX  ho.o  u 

PHDR 

HO.O  U 

delpx 

10.0 

DELPR 

10.0 

MUB 

-27.0  T 

Muv  -H’.o  T 

ZC 

2000.  T 

lat 

35.0 

DUMP 

.0  T 

OUTR 

I .0 

mode  2.o  u 

SL 

135.0 

dElth 

10.0 

MPR 

2.0 

tbp 

10.0 

pulse  .io+oot 

01 

25.6  C 

BAS 

10.0  T 

BAR 

10.0 

NIN 

-M5.0  T 

NOUT  -60*6  C 

TVC 

12.0 

GDOP 

H.H  C 

TGS 

15.0  T 

RON 

-2.0  T 

clip  «o  t 

SIGMA 

6.0  T 

pfa 

. I00-02T 

SP 

5.  T 

60 

. 0 | 0 7 N 

Alfa  .1277  n 

ZL  H 9 2 , 1 3 N 

ZCJ 

6152.  N 

PHIB  15.93  N 

CX  H882.39  N 

CT  H 885  t 53  N 

CL  H887.IH  N 

CS  H 8 8 1 . 8’  N 

CB  5077.36  N 

TSUR 

H 9 . 0 N 

TaV  36. H N 

NOCZ 

0 N 

|Rap 

0 N 

I AMOS 

1 N 

RANGE 

TIME 

LOSS  ANGLE 

signal  reverb 

REV 

s/mn 

EXCESS 

PROB 

( K T 0 ) 

(SEC) 

(OB)  IDEG) 

( Ob  ) 

( D8  ) 

LIM 

(OB) 

(DB) 

detect 

36. 0 

HB.  | B 

23’.’  1 ’ . H 9 

-89.9 

-82.6 

NO 

-2’.2 

-27.2 

.0000 

39.0 

H 9 . 3 H 

2H2 , 7 1 ’ . 1 9 

-92.7 

-8H.5 

NO 

-32.1 

•30.  1 

.0000 

<40.0 

50.50 

2h5 . 9 18,89 

-95. 9 

-85.8 

NO 

-3S.2 

-33.2 

.0000 

HI  .0 

5 | . 68 

2h’.2  18.62 

-”.2 

-83.2 

NO 

-38,5 

-36.5 

.0000 

H 2 • 0 

82.83 

252.6  18.37 

-102.6 

-28. H 

YES 

-69.8 

-69,8 

.0000 

H 3 . 0 

53. 9» 

256.6  18,12 

-106,6 

-3H.7 

Y t S 

-67, H 

-67. H 

.0000 

HH.O 

55.16 

260.6  17.91 

-110.6 

-31.2 

YES 

-7S.0 

-75.0 

.0000 

H 5 . 0 

56.33 

265.3  17.70 

-115.3 

-38.5 

YES 

-72,  H 

-72. H 

.0000 

M6.0 

57.51 

268.3  -17.53 

-118.3 

-HO.  9 

YES 

-73.  1 

-73.  | 

.0000 

H7.0 

58.68 

269.6  -17.35 

-119.6 

- H 2 . 9 

YES 

-72.5 

-72,5 

.0000 

H8.0 

59.88 

270.’  -17,19 

-120.9 

-H5.6 

YES 

-71.2 

-71.2 

.0000 

H9.0 

61.03 

272.2  -17.03 

-122.2 

-H7 . 9 

YES 

-70.  H 

-70. H 

.0000 

SO.O 

82.21 

273.5  -16.90 

-123.5 

-H9  • 8 

YES 

-70.  1 

-70.  1 

.0000 

8 1.0 

63.38 

2 7 H , B -16.77 

-12H.8 

-33.0 

YES 

-87. H 

-87. H 

.0000 

82.0 

6m  .56 

276.2  -16,65 

-126.2 

-35.9 

YES 

-BS.  9 

•85.9 

.0000 

83.0 

65. 7H 

277,6  -16,55 

-127.6 

-36.8 

yes 

-86,  H 

-86. H 

.0000 

8 H • 0 

66.92 

279.0  -16, H5 

-129. C 

-39.7 

YES 

- 8 H , 9 

-8  H , 9 

. 0000 

85.0 

68.0’ 

280. H -16,37 

-130. M 

-H  1 . 8 

YES 

- 8 H . H 

-8H.H 

.0000 

86 .0 

69.27 

281.9  -16.29 

-131.9 

-HH.O 

YES 

•83,8 

-83.8 

.0000 

57.0 

70.  M5 

283.5  -16.22 

-133.5 

-H6 . 5 

YES 

-83.0 

-83.0 

.0000 

56.0 

71.63 

285.1  -16.16 

-135.1 

-H8.5 

YES 

-82.8 

-82.8 

.0000 

59.0 

72.81 

285.9  16,10 

-135. y 

-50.8 

YES 

-81.8 

-81.8 

.0000 

60.0 

73.9’ 

286,7  16,06 

-136.7 

-52.3 

YES 

-8  1 • H 

-81.2 

.0000 

6 1.0 

75.  |7 

287.8  16.02 

-137.8 

•53.6 

NO 

-81.7 

-81.3 

.0000 

62.0 

76.35 

289.5  15.99 

-139.5 

-5 H . 9 

NO 

-82.5 

-81.8 

.0000 

• inputs,  run  s • 


PMDX»S  > PH0K«5 
MOOE»3  i «UN 


CALCULATIONS,  nun  s 


ueamx  is  the  siNtxi/x  response  function  foh  a continuous  line  array 
beamr  is  the  sinixi/x  response  function  for  a continuous  line  array 


the  letters  following  the  numbers  in  the  heading  on  the 

PERFORMANCE  PREDICTION  PAGE  SIGNIFY  the  following, 


u 

blank 

T 

c 


user  has  specified  this  input  value  on  this  nun 
set  OY  USER  ON  a PREVIOUS  RUN 

TYPICAL  VALUE  set  IN  ABSENCE  OF  USER  SPECIFICATION 
CALCULATED  INTERNALLY  BY  PROGRAM  (USER  OPTION) 

*KONG  value  SPECIFIED  BY  USER  (OUT  OF  RANGE  OF  ALLOW- 
ABLE VALUES.  THE  PRINTED  VALUE  IS  A TYPICAL  VALUE.) 

not  an  Input,  calculated  internally. 


i 


REVER8ERAT J UN  LEVELS  AT  T J ME  5 CORRESPOND  J NG  TO  2-AAV 
ACOUSTIC  TRAVEL  TIMES  BETWEEN  SOURCE  AND  TARGET, 

( 2 PINGS  --  TIME  BETWEEN  PINGS  15  IO«U0  SECONDS.) 
the  sound  propagates  between  source  and  target  ar 


CONVERGENCE  ZONE  --  20N£  | 


target 

UO-«AY 

bottom 

SURFACE 

VOLUME 

TOTAL 

PING  1 

PING  2 

range 

T I ME 

REVERB 

REVERB 

REVERB 

reverb 

REVERB 

REVERB 

KYD 

SEC 

DO 

Ob 

08 

DB 

UB 

DB 

Si  .0 

6 3 • u 6 

-88.2 

-89.2 

-3  1.7 

-31.7 

-31.7 

-50.8 

S2.0 

6m. 29 

-88.9 

-35 . 6 

-35. M 

-32.9 

-32.5 

-52.9 

S3.0 

65.S5 

-89.7 

-91.9 

-36.2 

-35.0 

-35.  1 

-53.9 

5 9.0 

66.72 

-90.3 

-MM  . 2 

-39. M 

-38.1 

-38.2 

-55.2 

5S.0 

67.93 

-91.2 

— M 7 . 6 

-Ml  . 7 

-90.6 

-MO.  7 

-56.9 

5 6 • 0 

69.19 

-92.3 

-So. 5 

-M3. 8 

-92.8 

-92. 9 

-57.5 

57.0 

7o.  MM 

-93.5 

-55.3 

- M 6 . 7 

-95. 9 

— M 6 . 2 

-58 . 5 

58.0 

7 1.67 

-9M  , 6 

-58.5 

— M 9 . 0 

-98.2 

-MB. 5 

-59.5 

59. 0 

72.90 

-95.9 

-61.9 

-50.9 

-50.5 

-50.6 

-72.3 

60.0 

7m.  12 

-97.9 

-69.2 

-52.5 

-52.2 

-52.2 

-79. S 

61.0 

75. 35 

-98.8 

-67.0 

-53.9 

-53.6 

-53.7 

-79.3 

62.0 

76.88 

-I0u.3 

-69.8 

-55.2 

-55.0 

-SS.0 

-81.7 

63.0 

77.01 

-101  .9 

-72.7 

-56. M 

-56.3 

-56.3 

-83.7 

68.0 

7 9 • U 3 

-103.7 

-75.6 

-57.5 

-57.9 

-57.9 

— 86.0 

6S.0 

80.26 

-105.9 

-78.3 

-58. M 

-58.9 

-58.9 

-88.0 

6 6.0 

8 1 . M9 

-107.2 

-81.0 

-59.  M 

-69.3 

-59.3 

-89.9 

67.0 

82.72 

-109.2 

-83.5 

-72.2 

-71.9 

-7  1.9 

-10  1.3 

68.0 

83. 9M 

-111.3 

— 86*0 

-79.5 

-79.2 

-79.2 

- 1 1 1 .0 

69.0 

85*17 

-113.9 

-88.5 

-76.7 

-76.9 

-76.9 

-112.1 
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REVERBERAT ION  bEVELS  AT  TIMES  CORRESPONDING  TO  2-«AY 
ACOUSTIC  TRAVEL  TIMES  BETWEEN  SOURCE  A11D  TARGET. 

( 2 PIN&S  --  TIME  d E T a £ £ N PINGS  IS  IC.JC  SECONDS.) 
The  SOUND  propagates  between  SOURCE  AND  TARGET  BY 


CONVERGENCE  2 ONE  --  ZONE  2 


TARGET 

T * 0 - »i  A Y 

bottom 

surface 

VOLUME 

TOTAL 

PING  1 

PING  2 

range 

TIME 

kl VERb 

REVERB 

REVERB 

REVERB 

REVERB 

REVERB 

K Y 0 

SEC 

0 o 

06 

Db 

DB 

DB 

DB 

103.0 

1 2 7 • 3 V 

- 3 u J • o 

• 3 C b . C 

-62.8 

-62.6 

-62.8 

-75.2 

109.0 

128.62 

- 3GG  • 0 

— 6 6 . b 

-67.1 

-69,2 

-69.5 

-76.3 

10S.0 

1 29 • o2 

-3uo  • 0 

-72.3 

-67.9 

-65.9 

-66.1 

-78.3 

106.0 

1 3 1 .09 

• 3 u w • 0 

-79.2 

-69.3 

-67.8 

— 66,1 

-79.6 

107.0 

132.2b 

-36o»u 

-75. B 

- 7 v • 8 

-69.3 

- 6 9 , 6 

-81.1 

108. u 

133.97 

- 3uu  • u 

-77.3 

-72.2 

-70.  7 

-7  1 . u 

-82.3 

109.0 

1 3 9 . OB 

• 30 u • 0 

-79.9 

-73.5 

-72.  1 

-72.5 

-63.9 

1 10. 0 

1 35.09 

-3Cw  • 0 

-81.0 

-79.6 

-73,9 

-73.7 

-89.5 

I 1 1 .0 

137.10 

-300*0 

-82.6 

-75.7 

-79,6 

-79.9 

-85.6 

112.0 

138*30 

-300*0 

-89 . 2 

-76.7 

-75.7 

-76. U 

-66.5 

113.0 

139*66 

-300*0 

— 86.5 

-78.6 

-77.5 

-78.0 

-87 . 9 

*18.0 

190*09 

- 30u • u 

-89.2 

-80.  1 

-79.1 

-79.6 

-88.3 

1 1 5 . 0 

192.12 

•30u  • Ci 

-91.0 

-81.9 

-80.3 

- 6 U , 9 

-89.2 

116.0 

193.39 

- 3Gu  • 0 

-92 . 6 

-82.5 

-81.5 

- B 2 , 1 

-90.0 

117.0 

199.57 

-300*  C 

-99.2 

-83.7 

-62 . 6 

-83.3 

-90.6 

118.0 

195.80 

— 3 0 \j  • 0 

-95.8 

-89.8 

-83.7 

-89 . 9 

-91.6 

119.0 

197. u3 

- 3 u u • 0 

-97.9 

-B5.8 

-89 . 7 

•85.5 

-92.3 

120.0 

198.25 

- 300  • u 

-96.9 

-86 . 7 

-85 . 6 

-86.9 

-93.1 

12  1.0 

199.98 

— 3 0 y • >J 

- 1 Ob.  3 

-87.5 

-86.9 

-67. 3 

-93.7 

122.0 

l So. 7 1 

- 3 0 o • o 

-10  1.9 

-88.9 

-67.3 

-88.2 

-99.9 

123.0 

151.99 

-30  J • o 

-103.9 

-89.2 

-88.1 

-89.1 

-95.0 

129.0 

153.1b 

- 30 y • 0 

-1 05. C 

-90  • C 

- fa  B • 9 

-B  9.9 

-95.7 

1 2S.0 

159.09 

•3uo»o 

— 1 o 6 . 5 

-90 . 8 

-90.5 

-90.7 

-103.9 

126.0 

155.62 

- 30w  • 0 

-108.1 

-91.6 

-91.9 

-91.5 

- 1 07 . 2 

127.0 

1 5 6 • o S 

- 3 0 - • 0 

- 1 C 9 . 7 

-92.3 

-92.1 

-92.2 

- 1 CB  . 5 

128.0 

1 58 . b 7 

- 3G  J • G 

-111.2 

-93.0 

-92.9 

-93. u 

- 1 09  . B 

129.0 

1 59 • 3b 

-3iiu*w 

-112.7 

-93.7 

-9J. 6 

-93.6 

- 1 11.0 

130.0 

1 6Q . S3 

- 3 0 u • 0 

-119.2 

-99.3 

-99 . 2 

-99.3 

-112.3 

131.0 

1 6 J .76 

— 3 0 o • 0 

- 1 1 5. 6 

-95.0 

-99,9 

-99.9 

-113.5 

132.0 

162.96 

-300*0 

-117.0 

-95.6 

-95.5 

-95. 6 

-119.6 

133.0 

169.21 

-300*^ 

- 1 18.9 

-103.9 

- 1 U 3 • 1 

- 1 03. 3 

-116.1 

139.0 

165.99 

- 3 0 u • 0 

- 1 19.9 

-107. 3 

-10  6.9 

- 1 «7  • 1 

- 1 22  . B 

1 35.0 

166*67 

- 3 0 o • w 

-121.3 

—106.5 

-106.2 

-106. 3 

- 1 29  . b 

136.0 

167*09 

-300*0 

-122.6 

-109.8 

-109.5 

- 1 b 9 , 6 

-125.2 

137.0 

169*12 

-30w  • 0 

-129.0 

- 1 1 1 . 1 

-l  IU.  7 

-110.8 

-127.0 

138.0 

1 70.35 

- 30u  • 0 

-125.9 

-1)2.3 

-112.0 

-1)2.1 

-128.2 

139.0 

171.57 

- 30u  • o 

-126.8 

- 1 13.6 

-113.3 

-113.9 

-129.7 

190.0 

172*60 

- 3 G u * 0 

- 1 2b  . 1 

-119.8 

-119.6 

-119.6 

-300.0 

1 9 1 • 0 

1 7 9 . b 3 

-30o*C 

-129.5 

-116.1 

-115.9 

-115.9 

-300.0 

f 
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ACTIVE  PERFORMANCE  prediction 

WESTERN  N,  ATLANTIC  i.INTER  SOUND  SPEED  PROFILE*  NUC  TN  1U06,  P,  97 


CONVERGENCE 
OUTP  *0 

ZONE  - - 

Zx 

ZONE  1 
20.  0 

ZTG 

300.0 

FREQ 

2.0 

La  A 

5.0 

V A I 

1 S • 0 

salt 

35. C T 

hS 

5.0  C 

ZBM  17372. 

REFLS 

1.0  T 

ZONES 

2.0 

PhDx 

5.0  U 

PHDR 

5.0  U 

DELP* 

10.0 

OElPR 

10.0 

MuB 

-27.0  T 

71 U V 

89.3  T 

zc 

2000.  T 

l*t 

35.0 

DUMP 

.0  T 

UUTR 

I .0 

node 

3.0  U 

SL 

135.0 

delth 

10.0 

MPR 

2.0 

TBP 

10.0 

Rulse 

• 1 0 + 0 u T 

01 

25.6  C 

BPS 

10.0  T 

bwR 

10.0 

N I N 

-95.0  T 

NOUT  - 

6 0,6  C 

TVC 

12,8 

gdop 

9.9  C 

TGS 

15.0  T 

RON 

-2.0  T 

Clip 

•C  T 

S I GMA 

6.0  T 

pfa 

. 10C-02T 

SP 

5.  T 

GO 

.0107  n 

Alfa  . 

1 2 7 7 N 

ZL  892.13  N 

ZC  J 

6152.  N 

PM  I 8 15.93  N 

CX  9082.39  N 

CT  8885. S3  N 

CL  9887.18  N 

cs  9 8 ( 

6 1,89  N 

CB  5077.36  N 

TSUR 

9 9.3  « 

T A V 

36.8  N 

NOCZ 

C N 

[Rap 

U N 

I AMOS 

1 N 

range 

T I m£ 

LOSS 

ANGLE 

SIGNAL  reverb 

rev 

S/MN 

EXCESS 

PROB 

l K Y D ) 

( SEC  ) 

103) 

( deo  i 

( DB  ) 

I DB  ) 

LIN 

(08) 

( DB  ) 

DETECT 

5 1.0 

63.08 

180.1 

3 .98 

-33.  1 

-31  . 7 

TES 

6.0 

6.0 

. 8385 

52.0 

68.2V 

1 86 . 6 

7.78 

— 36.6 

-32.9 

YES 

• 2 

• 2 

.5158 

S3 . 0 

6S  . SS 

18  9.8 

8.76 

-39.8 

-35. U 

YES 

-.9 

- • 9 

. 9796 

50  .0 

66.72 

195.1 

9.70 

-95.  1 

-38.1 

YES 

— 2,6 

-2.6 

.3338 

SS  • 0 

67.93 

198,9 

10.52 

-98.9 

-90*6 

YES 

-8.1 

-9.  1 

.2572 

S 6 • 0 

69  . 1 R 

2 3 3,9 

11.30 

-53.0 

-92.8 

YES 

-6,9 

-6.8 

.1512 

B7  • Q 

7 0,99 

207.2 

3.50 

-57.2 

-95.9 

YES 

-7.2 

-7.2 

.13  16 

sa  .0 

71.67 

209,6 

3.83 

-5  V.  6 

-98.2 

YES 

-7.6 

-7.6 

.1229 

S9.C 

72.90 

211.8 

3.35 

-61.8 

-50.5 

YES 

-7.9 

-7.8 

.1186 

6 0 » Q 

7 8 . 1 2 

2 1 3 , 8 

3.27 

-63,8 

-52.2 

YES 

-8.6 

-8.5 

. 1090 

91*0 

7 S . 3 S 

2 15,6 

3.19 

— 65,6 

-53 . 6 

NO 

-9.5 

-9.  I 

,0896 

62. 0 

7 6.58 

2 1 7.0 

3.  1 1 

— 67,0 

-55.0 

no 

-10,9 

-9.7 

.0760 

63.0 

77.81 

2 l 9 , 0 

3.03 

— 69*0 

-56,3 

no 

-11.3 

-10.3 

.0612 

68.0 

79.03 

220.5 

2.95 

-70.5 

-57.9 

NO 

- 1 2 , 0 

-11.1 

.0989 

6S  . 0 

BO. 26 

222. L 

2.87 

-72 . 3 

-58.9 

NO 

-13.9 

-11.9 

.0255 

66.0 

6 1.09 

223.8 

2.79 

-73.8 

-59.3 

nO 

-19.5 

-12.8 

.0199 

67.0 

82.72 

228 . 8 

2.72 

-79.8 

-7  1.9 

N 0 

-19.2 

-12.2 

.022  1 

68.0 

83.98 

226,0 

2 .69 

-76.3 

-79.2 

N 0 

- 1 8.5 

-13.5 

.0173 

69 . C 

85.17 

2 27.3 

2.56 

-77.3 

-76.9 

NO 

-16.7 

-19.7 

.0127 

117 


ACTIVE  PERFORMANCE  prediction 

WESTERN  N,  ATLANTIC  WINTER  SOUND  SPEED  PROFILE.  NUC  TN  1006,  P.  97 


CONVERGENCE  ZONE  --  ZONE  2 


OUTP 

• 0 

Zx 

20.0 

ZTG 

300  • 0 

FREW 

2.0 

L PA 

S.O 

V A I 

)S>0 

salt 

35.0  T 

hs 

5.0  C 

ZBM  17372* 

REFLS 

1 .0  T 

ZONES 

2*0 

PrtDX 

5.0  U 

PHDR 

5.0  U 

DELPX 

10.0 

delpr 

10.0 

MoB 

-27.0  T 

N U V 

99.0  T 

ZC 

2000.  T 

LAT 

35.0 

DUMP 

.0  T 

OUTR 

1 >0 

mode 

3.0  U 

SL 

135.0 

delth 

10. 0 

MPR 

2.0 

TBP 

10.0 

pulse 

. 10.00T 

D 1 

25.6  C 

8 a S 

10.0  T 

BAR 

10.0 

N | N 

-NS . 0 T 

nout  - 

60.6  C 

1 V c 

12.0 

GDOP 

9.9  C 

TGS 

15.0  T 

RON 

-2*0  T 

Clip 

•0  T 

S I GMA 

6,0  T 

PFA 

. 1 00-02T 

SP 

5.  T 

<>0 

.0107  N 

Alfa  . 

1 277  N 

ZL  992.13  N 

ZCJ 

6152*  N 

PHI  B 1 

5.93  N 

CX  9862,39  N 

CT  ‘,885.53  N 

CL  9887,19  N 

CS  9881.89  N 

CB  5077.36  N 

TSUR 

**9.0  '« 

Tav 

36.9  N 

NOCZ 

0 N 

1RAP 

0 N 

1 AMOS 

1 N 

range 

T 1 Mt 

LOSS 

ANGLE 

SIGNAL  REVERb 

REV 

S/MN 

EXCESS 

PROB 

IWD) 

( SEC  1 

1 OB  1 

(DEG) 

( Db  i 

(DB  ) 

L 1 M 

(DBI 

(DB) 

DETECT 

103.0 

127. 3» 

213.7 

9.35 

-63.7 

-62.6 

NO 

-9.U 

-2,0 

. 3706 

1 Or  . 0 

1 2 8 . 6 2 

2 1 9 . 1 

3.97 

-69.1 

-69.2 

NO 

-».  I 

-7.  1 

.1395 

1 US . 0 

129.82 

222.6 

7.69 

-72.6 

-65.9 

NO 

— 1 2 , 9 

-10.9 

.0603 

1 U6  • 0 

|3l  .0*1 

22i  . 9 

8 .93 

-79.9 

-67.6 

NO 

-n.6 

-12.6 

.0208 

IU7.0 

1 32. 26 

22  7 . 1 

8.99 

-77.1 

-69 . 3 

NO 

- 1 6 . 7 

-n.7 

.0128 

1 08  • 0 

1 33.97 

229.  7 

9.92 

-79.7 

-70.7 

NO 

-n.2 

-17.2 

.0029 

1 09 . 0 

1 3r . 6» 

2 32  . 1 

9.86 

-62.  | 

-72.1 

NO 

-2  1 • 6 

— 19,6 

.0000 

1 10.0 

1 35. 89 

239.2 

10.28 

-89,2 

-73.9 

N 0 

-23.7 

-21.7 

.0000 

1 1 1 .0 

137.(0 

236.3 

10.70 

-86.3 

-79.6 

NO 

-25.7 

-23.7 

.0000 

112.0 

1 3b. 30 

236.9 

11.10 

-88.9 

-75.7 

no 

-27.6 

-25.8 

.0000 

113.0 

139.66 

2 9 v . 6 

3,62 

- 9 0 . 6 

-77.5 

no 

-30.  U 

-28.0 

.0000 

119.0 

no.  87 

292. 1 

3.59 

-92.1 

-79.1 

NO 

-31.5 

-29. S 

. 0000 

1 1S.0 

112.  | 2 

29  3 . 6 

3.55 

-93.6 

-80*3 

NO 

-33.0 

-31.0 

.0000 

1)6.0 

no. 39 

2 9 9,9 

3 , 5 | 

-99.9 

-8  1.5 

NO 

-39.3 

-32.3 

.ouoo 

117.0 

199.57 

296.2 

3.97 

-96.2 

-82.6 

N 0 

-35,6 

-33.6 

.0000 

I 1 B • 0 

ns. BO 

297.5 

3.93 

-97.5 

-63.7 

NO 

-36.9 

-39.9 

. 0000 

119.0 

197.03 

298,7 

3 . 90 

-98.7 

-89.7 

N 0 

-38.1 

-36.1 

.0000 

120.0 

198. 2S 

299,6 

3,36 

-99.8 

-6S.6 

NO 

- 3 V . 2 

-37.2 

. 0000 

12  1.0 

199.96 

250.9 

3.32 

— 1 0 U • 9 

-86.9 

NO 

-90.3 

-3  8.3 

. 0000 

122.0 

1 SO. 7 1 

252.0 

3.28 

- 1 02.  0 

-87.3 

NO 

-91.9 

-39.9 

.0000 

123.0 

1S| .99 

253.1 

3.25 

-103.  I 

-68.1 

NO 

-92,9 

-90.9 

. 0000 

129.0 

1 S3.  1 6 

259.1 

3.21 

- 1 DR  . 1 

-68.9 

NO 

-93,9 

-91.9 

.0000 

125.0 

159.39 

255.  ) 

3.17 

-105.  1 

-90*5 

NO 

-99,9 

-92.9 

• OUOO 

126.0 

1 S 5 . 6 2 

256 . 0 

3.13 

-106.  „ 

-91.9 

NO 

-95. 9 

-93,9 

.0000 

127.0 

1 &6 . 8S 

256.9 

3 . U 9 

-106.9 

-92.  1 

NO 

-96.3 

-99.3 

. ouoo 

128.0 

1 Sb . 07 

257.6 

3 , u 6 

-107.6 

-92.9 

NO 

-97.2 

-95.2 

• OOOQ 

129.0 

lS9.3u 

256.7 

3.02 

-108.7 

-93.6 

no 

-96.1 

-96.  | 

■ OOOO 

130.0 

160.53 

259.6 

2.98 

-10  9.6 

-99.2 

NO 

-99.0 

-97.0 

. 0000 

13  1.0 

161.76 

260.5 

2.99 

-1  ID. 5 

-99.9 

NO 

-99,8 

-97.8 

. OGGC 

132.0 

162.96 

26  1 . 3 

2.90 

-111.3 

-95.5 

NO 

-50.7 

-98.7 

. 0000 

133.0 

169.21 

262.  1 

2.67 

-112.1 

- 1 9 3 . 1 

NO 

-51.5 

-99  . s 

. 0000 

139.0 

165.99 
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